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Abstract 

The perturbed-stationary-states (PSS) method, modified to 
include electron translation factor (ETF) effects, has been 
used to calculate detailed cross sections for charge 
transfer and direct excitation in H*-H(1s), He?*-H(1s) and 
H*-He*(1s) collisions at low to intermediate energies (E< 10 
KeV/amu). The close-coupled equations are solved using a 
classical trajectory description of nuclear motion, and 
electron translation factors appropriate to molecular states 
have been constructed using the method of Switching 
functions, after the formulation of Delos and Thorson(12). 
Switching functions for molecular states of H,* and HeH?* 
Systems have been determined by the "optimization" method 
developed in earlier work from this laboratory (29). 

Our cross sections and other calculated results are 
compared with other theoretical studies and with recent 
experimental findings. A brief summary of the relevant 
conclusions is as follows: 

1) Ht-H(1s) collisions 

(a) Our differential and total charge transfer 
cross sections are in good agreement with 
earlier calculations and with experimental 
values. 

(b) The H(2p) charge transfer cross section is 
in good agreement with the molecular. state 
calculations of Crothers and Hughes (48) and 


with recent experimental values. Discrepancies 
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2) 


between our results and those of Refs.(48) may 
be traced directly to differences between the 
ETF descriptions used here and in Refs.(48). Our 
results appear to be in somewhat closer 
agreement with the experimental measurements. 
(c) Our H(2s) charge transfer and direct 
excitation cross sections differ substantially 
from those found by Crothers and Hughes(48) and 
our results are in excellent agreement with the 
most recent experimental data. 

(d) We have studied the convergence of cross 
sections as a function of increasing basis size 
and the excitations to higher (Rydberg) levels. 
We find that a small but non-negligible portion 
of the flux always escapes to the highest levels 
available in the basis set and conclude that the 
most important ionization mechanism in these 
collisions involves a "ladder-climbing" process 
of sequential excitations involving such states. 
He?*-H(1s) collisions 

(a) Our total charge transfer cross sections are 
in good agreement with those calculated by 
Winter and Hatton(73) (see also Hatton, Lane and 
Winter«/4)) and with. recent experimental 
measurements, but disagree strongly with the 
values reported by Vaaben and Taulbjerg(72). As 
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3 ) 


corrections greatly improves convergence of the 
results as function of basis size, and also that 
the results seem to converge on _ slightly 
different limits. 

(b) Our detailed excitation/charge transfer 
cross sections for individual states differ 
slightly from those of Winter et _al(73), though 
both are in good agreement with experimental 
values. 

H*-He*(1s) collisions 

(a)Our total charge transfer cross sections are 
14%-32% smaller than those of Winter, Hatton and 
Lane(/5) the ditterencesvappear to result from 
differences in ETF descriptions used. 

(b) Even larger discrepancies (as much as a 
factor of two in one case) are found between our 


(sia lwer )TeErossa Sect lONSTOMOriPundividual “state 


excitations, and those given by Winter et al; 


these differences persist as basis size is 
increased and are traceable toin the. ETF 


descriptions used. 
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PREFACE 


When I handed the manuscript of my thesis to Prof. 
Thorson, he corrected and extensively revised my Japanese 
English. Therefore not a single Japanese phrase in my 
Original manuscript survived in the revised thesis. This 
preface is supposed to be my last territory where no one is 
allowed to trespass, so that at least a part of my thesis 


contains my original word. 


inehduky o1b975 forl srecemvedraihetter wfrom \ProficeThorson 
after I asked him about the possibility of my study with him 
at the Univ. of Alberta. He kindly provided me with a lot of 
information about graduate study in Chemistry in this 
University land "Ckwas ‘so excited after i read the letter. A 
few days later, I received a telegram from Dr. S. ° Davis 
confirming my acceptance as graduate student at the Univ. of 
Alberta. When I received it, the telegram man said to me "a 
couple of words in the telegram were lost somewhere between 
Vancouver and Tokyo. We are now trying to find out what they 
are". I was worried that the missing word was "not" in the 
sentence "for not accepting you". But fortunately, it was 
not "not". Then I rushed into the bookstore and looked at 
the map of Canada to find out exactly where Edmonton is. 
When I told my father about my plan to study in Canada, he 


Said "It must take a long time to finish it and maybe I 
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won't see you in this world, but good luck". I am pleased to 
Say that he is doing fine now. 

On Christmas Eve, 1975, I flew from Tokyo to Edmonton 
with more than a hundred of Vietnamese orphans who seemed to 
me to be one or two years of age. It was my first shock of 
seeing what is happening in the world. I arrived in Edmonton 
at 5:30 PM and the temperature was -7 C, which wasS not’ so 
bad. But I soon realized it was just a short break. I stayed 
on a huge "queen"-size bed at the MacDonald Hotel during the 
first night in Edmonton without any blink of sleep, thinking 
of my study, life, etc. (I thought the MacDonald Hotel must 
be a cheap place for a poor graduate student, since it was 
suggested Ilestayotherevby’) thesChemeeDept., officer. yActually, 
it was not a cheap place.) Next morning, on Christmas day, I 
met Prof. and Mrs. Thorson for the first time at his office 
and thought, "Gee! Rte 1S “NOt stair to work with him as a 
graduate student, because of such a huge difference in 
physical size and mass". But it was just too pessimistic, 
because I am still surviving. As a matter of fact, I really 
enjoyed working with him, and learned not only his attitude 
towards science, but also his attitude towards life which I 
had never experienced before. In the last five years, I have 
met many nice friends inside and outside the university 
community, and they are among the most valuable benefits 


that I obtained in Canada. 
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I would like to dedicate this thesis to my teacher, 
colleagues, friends and parents. The work in this thesis was 
not done by myself alone, but in collaboration with many 
people. 

Prof. Thorson guided me in the proper direction from 
kindergarden to the PhD as a Chemistry professor and an 
English teacher, and from Lake Louise to Sentinel pass as a 
mountain guide. 

Drs. S. Knudson and J. Delos taught me a lot about physics. 
SeieChor, JIsvaDavis, ade Rankinivand SenBartontcreated an 
enjoyable atmosphere in this laboratory. The members in the 
Theoretical Div. (Drs. Birss, Huzinaga, Fraga, Clarke, Mrs. 
J. Charter and colleagues) supported me in many ways. And my 
friends, Oh! what nice people. 

Without those people, I could not have completed my work. I 
cannot thank them enough. 

Following the North American tradition, I also would like to 
thank my wife, Keiko, and my daughter, Kana for substantial 
moral support and sharing life with me, and of course, would 
like to thank to our parents for giving me many kinds of 
Support. 

I should mention special thanks to Mrs. Thorson who has 
considered us so warm-heartedly that Kana feels as if she 


were her "auntie" in Edmonton. We really appreciate her. 
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I. INTRODUCTION 


A. PURPOSE 
We have investigated fundamental processes occurring in 
Slow ion-atom collisions. The systems we have studied are of 
the simplest 3-body type (one electron and two nuclei), but 
serve aS adequate prototypes for studies on the _ basic 
aspects of the general ion-atom scattering problem with 
which we are concerned in this thesis. Although these 
primitive systems are not alwaysS convenient from an 
experimental standpoint, they are most useful for 
theoretical purposes because 
(1). The electronic wavefunction for the one-electron, 
two-nucleus system is known exactly so that we can avoid 
uncertainties about the accuracy of the electronic basis 
functions; 
(2). The fundamental dynamical effects which - cause 
physical changes ate. the System can be clearly 
represented without any interference of effects arising 
from multi-electron system phenomena (e.g., 
autoionization, configuration mixing effects at curve 


crossings,etc.) 


A major objective of this thesis is..to determine . the 
degree and importance of the electron translation factor 


effects on calculated cross sections for charge exchange and 
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direct excitation processes in one-electron systems, in 
particular, the H,* and HeH’?* systems. Only recently has it 
been recognized that these effects are of importance, and 
that their adequate treatment may require quite elaborate 
descriptions of electron translation factors. A second 
objective of this work is to examine the overall validity of 
the close-coupling method as it is normally employed, which 
neglects loss of flux due to ionization. We present some 
evidence that a small but not negligible portion of the 
Scattering flux in these collisions would correspond in 
reality to ionizing events, and we can make some estimate of 
its magnitude. 

To study these collisions, we used the impact parameter 
or classical trajectory multistate perturbed stationary 
State (PSS) method, corrected for electron translation 
effects by following the recent formulation of Delos and 
Phorson Gir, 127) 3 This "method @treats e*thehenuciears* motion 
classically, and the electronic problem is solved by the 
close-coupling method ; quantum mechanically, using 
molecular states as basis states. 

We have mainly applied the theory to investigate charge 
exchange and direct excitation for the symmetric H*-H(1s) 
System and the asymmetric He’*-H(1s) and He*(1s)-H* 
collision systems. Reliable experimental measurements have 
been performed recently and extensive molecular state 
calculations have also been reported. These serve as useful 


comparisons with the present theory. 
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B. THEORETICAL BACKGROUND 

In a low to intermediate energy collision, where the 
relative velocity of the nuclei is smaller than the orbital 
electron speed, the colliding system may be regarded as 
forming a quaSi-molecule. If this is the case, we might 
Suppose , aS Born and Oppenheimer did for their treatment of 
molecules(1), that the electron system can adjust smoothly 
to the slowly changing nuclear position. It would appear 
that the best basis in which to expand the scattering wave 
function would be the set of Born-Oppenheimer(1) (BO) (fixed 
nuclei) adiabatic molecular electronic states, i.e., the 


expansion 


where the BO states (P:R) are the eigenfunctions of the 


electronic Hamiltonian at each (fixed) nuclear position R, 


>> >> >> 
° ° = Hs 
Hy (eR) o, (ee) epee Ar ) 


and depend parametrically upon R. 
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Substitution of this expansion into the complete 
(time-independent) Schroedinger equation leads to coupled 


equations for the components x (2) of a vector Xx ; 


=) ; > a abe Ey (R) 1-3 
(2u) “[-ihlV, + P(R)] x(R) + €(R)X(R) me 


where £ is the diagonal matrix with elements €,(R) . In 
this formulation, all couplings (and electronic transitions) 
result from the nonadiabatic coupling matrix P(R). This 
matrix has both angular and radial coupling components; the 


radial component, for example, is defined as 


> = > 
=< -ihV > I-4 
PL fR) k| ih sa 


and is supposed to represent the effects of deformation of 
the electronic wavefunctions due to changes in the 
internuclear distance R. This formulation of Scattering 
problems for slow collisions is called the 
perturbed-stationary-states (P.S.S..) method, and 
anvestigations uSingsite« gow back «to.s1953,.0With.»work by 
Bates,Massey and Stewart(3). The method has been widely used 
in studies of elastic and inelastic collisions at thermal 
energies. However, it has been known for a long time that it 
suffers from some fundamental difficulties, which become 


quite evident when the method is applied to practical 
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problems involving charge transfer or direct electronic 
excitations. Firstly, it can be shown generally that the 
asymptotic boundary conditions appropriate to a scattering 
problem are not properly satisfied by the molecular. state 
expansion given in Eq.(I-1). Secondly, the coupling matrix 
elements Pan (®) (both radial and angular components) entering 
into the coupled equations(I-3) are found to depend markedly 
on the reference origin chosen for electron coordinates 
(i.e., they do not exhibit proper translational invariance), 
and furthermore they show a quite unphysical behavior; for 
example, some matrix elements do not tend correctly to zero, 
but to constant asymptotic values, as R- 9c. Also, some 
matrix elements.,havey- unreasonably wlarge: magnitudes at 
intermediate nuclear separations, and may be expected to 
give cross sections which are much too large. 

“All the above defects in the PSS method originate in 
the fact that the Born-Oppenheimer expansion(I1-1) entirely 
ignores the translational motion of an electron due to 
motion of a nucleus to which it is attached. Suppose we 
consider any reference origin (for example, the centre of 
mass), with respect to which a nucleus 1S moving: an 
electron bound to that nucleus is carried along with it, and 
a factor describing that translation for the electron is a 
necessary part of the electronic wavefunction. But no such 
factor appears in the Born-Oppenheimer expansion. 

Numerous attempts have been made to reformulate PSS 


theory... correctly..by including such electron translation 
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factors CET)? In the pioneering work of Bates and 
McCarroll(4), this is achieved in a manner appropriate to an 
atomic state (i.e., a single-centre orbital): for atomic 
orbitals associated with nucleus A, a factor describing 
translation with A is attached, and correspondingly a factor 
describing translation with B for orbitals centred on 
nucleus B. If a molecular state is asymptotically correlated 
tO8an “atomielstate  4on  *nucleus®* A*°C(or ‘on (B), "Bates and 
McCarroll assign to it the corresponding atomic translation 
factor. (Some additional complications are involved for 
molecular states of symmetric systems like H.* which have gq 
Or u symmetry and do not correlate in a one-to-one’ fashion 
with unique atomic states). Most attempts to correct the 
defects of PSS theory have been based on the Bates-McCarroll 
scheme or minor formal modifications of it. 

The Bates-McCarroll(4) description of ETF's ensures 
that asymptotic scattering boundary conditions and 
translational invariance requirements are satisfied; a 
formally correct theory is then obtained and defects such as 
the incorrect asymptotic behavior of coupling matrix 
elements are removed. However, physically, in the region of 
small to intermediate internuclear distances where the 
system forms a quaSimolecule, the electron moves under the 
influence of the molecular field: it belongs completely to 
neither nucleus but is partially attached to both. Therefore 
a proper ETF for molecular states should have molecular or 


two-centre character. 
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Thorson and Delos (11,12) have shown that ETF's 
appropriate for molecular states can be constructed using 
the idea of switching functions introduced originally by 
Schneiderman and Russek(9); the resulting formulation meets 
all requirements of translational invariance and scattering 
boundary conditions, and _ leads to coupled equations 
identical in form to those of PSS theory (Eq.(I-3)), except 
that the coupling matrix P is corrected for the effects of 
inclusion of electron translation factors. This scheme 
effectively recovers the same asymptotic description of 
ETF'S as the method of Bates and McCarroll does for 
Single-centre states, but the flexibility of the switching 
Function Balso. *permiter ETEAS to Pdescrabest theptwo-centre 
character of a molecular state at finite internuclear 
distances. Initially, Thorson and Delos assumed that a 
common switching function may be used for all electronic 
States of a given system, but more generally it may be 
expected that the most appropriate switching function will 
be different for each distinct molecular bound state. More 
recently, Delos and Thorson(12) have given a formulation of 
slow collision theory, based on this more general 
assumption, within the framework of a classical trajectory 
theory. 

However, these formulations do not fully specify the 
choice of molecular state switching functions to be used, 
and as a matter of fact their selection is an important 


practical problem since the nonadiabatic coupling matrix 
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elements for many electronic transitions are very sensitive 
to the switching functions used. 

AS @iag fvests.pane cofwlthissefthesrsp  itwsits eshownsthat 
several independent lines of thought can be used to make a 
consistent and definite choice of the proper molecular state 
Switching functions for the systems considered. 

The remaining part of the thesis then proceeds to a 
Gest maok otheheconurectederPSSvemethod by applying it to 
calculations of cross sections for the following collision 
processes: 

a. Symmetric System 
H*+ H(1s) ——> H*+ H(nl) (Direct excitation) 
H*+ H(1s) ——> H(nl) + H* (Charge exchange) 
b. Asymmetric System 
He?*+ H(1s) —— > He?*+ H(nl) (Direct excitation) 
He?*+ H(1s) ——»> He*(nl) + H* (Charge exchange) 
c. Asymmetric System 
He*(1s) + H* ——> He*(nl) + H* (Direct excitation) 


He*(1s) + H* ——> He?* + H(nl)(Charge exchange) 


Our results are compared with other recent molecular 
state calculations and with experimental measurements. 

In the symmetric H*-H(1s) collision system, I have made 
some studies with basis sets involving a large number of 
molecular states (up to 16 ungerade states), in an effort to 


test the limits of the close-coupling method in this 
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context. Close-coupling methods at present use only bound 
State basis functions and give no account of flux loss due 
to 1onizing events. The results of these studies suggest 
that a small but not completely negligible portion of the 
collision probability, even at the lower collision energies, 
will in reality appear as ionization by a sequence of 
"ladder-climbing" excitations, and we can form a= rough 


estimate of the magnitudes involved. 


C. HISTORICAL SURVEY 

Unresolved problems) still remain fOr} tithere st udyn7 of 
ion-atom collisions in the one-electron prototype systems 
despite their simplicity. In all energy regions there are 
some major discrepancies among theories, among different 
experimental measurements, and between theory and 
experiment. These systems have been extensively studied; 
specifically more than fifty theoretical and more than ten 
experimental investigations have been reported for _ the 
H*-H(1s) collision system and only slightly fewer for the 
He?*-H(1s) system. Here we will not attempt a survey of all 
the work reported, but will concentrate most attention on 


Studies of low to intermediate energy collisions. 
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1. Proton-Hydrogen Atom Collisions 

Numerous theoretical studies have been reported which 
are most likely to be valid at energies above the 
intermediate region(E210 Kev), and these have obtained 


reasonably accurate cross sections for higher energies. 


Methods used include Born and distorted-wave 
approximations(50),Glauber's method(51,52), the eikonal 


approximation(54), impulse approximation(53), Faddeev's 
method(55), and close-coupling treatments involving 
pseudo-states(56,57,58), as well as a variety of hybrids of 
these(59). However, all these methods fail increasingly at 
lower collision energies, giving wrong positions for cross 
section maxima and in many cases qualitatively quite 
incorrect behavior. The remedy for these difficulties is the 
use of molecular state descriptions which should be more 
appropriate below 10 Kev. 

A pioneering study was done by Bates and Williams(40) 
using the straight-line impact parameter approximation and a 
three-state PSS close-coupling theory (corrected by 
Bates-McCarroll ETF's). Their calculations showed that at 
low collision velocities the most important inelastic 
process is the population of H(2p) via 2pQ -2pt% rotational 
coupling. (Their work has been of major importance in 
understanding primary excitation events in more complex 
lon-atom systems--especially in respect to inner-shell 


vacancy production processes). 
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Later, impact parameter multichannel PSS-type 
calculations were done by Rosenthal(45), Chidichimo-Frank 
and Piacentini(46), and Schinke and Kruger(47-a). All these 
Studies neglected ETF effects as the uncorrected PSS method 
does and therefore the results lack the required Galilean 
invariance with respect to the choice of reference origin 
for electron coordinates; in addition, radial couplings were 
entirely neglected in these studies. 

A more sophisticated molecular state calculation has 
been reported recently by Crothers and Hughes(48). Their 
treatment of ETF effects is based on a modification of the 
Bates-McCarroll scheme: they construct linear combination 
fragments of g and u molecular states which correspond 
asymptotically to one-centre atomic states, and with each of 
these fragments they associate the corresponding 
Bates-McCarroll ETF; finally they recombine the modified 
fragments to form g and u molecular basis states. However, 
this approach has the same general defects as_ the 
Bates-McCarroll scheme because the "atomic fragments" are 
not really one-centre atomic functions at finite 
internuclear separations but also have a two-centre 
molecular character. To offset the resulting deficiencies, 
Crothers and Hughes introduce a variable parameter )(R) in 
the ETF, and determine its optimal value at each R by the 
Euler-Lagrange variational principle(21). [It can be shown 
that their procedure is equivalent (in low-velocity limit) 


to the use of a crude switching function(33).] 
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The results obtained by Crothers and Hughes show much 
closer agreement with experiment than do the earlier 
Studies, which completely ignore ETF effects. However, the 
Gefects in their somewhat ad hoc description of ETF effects 
do have significant effects on their calculated cross 
sections in certain cases (as I show in this thesis): in 
addition, they incorrectly used Hermitian averages of 
non-Hermitian couplings which appear in the theory, although 
it can be shown that the non-Hermitian character is 
necessary to preserve unitarity, and in some cases this may 
lead to incorrect values for individual state cross 
sections. 

One of the purposes in my thesis investigation was to 
clear up ambiguities in the 1 KevV-10 KeV region resulting 
from these defects, by using a rigorous theory with more 
detailed molecular state switching functions to describe ETF 
effects, including an attempt to demonstrate convergence of 
the molecular state close-coupling results as a function of 
basis size. 

(Some studies have been carried out for energies below 
1 KeV, but these have been primarily concerned with the 
Gifferential cross sections for charge exchange and Lyman- 

or excitation and its relation to the classical 
trajectories used; below 1 KeV it can be shown that effects 
of deviation from a straight-line (impact parameter) path 
are important (see, for example, papers by Knudson and 


Thorson(42), Bates and Sprevak(43), and Schinke(47-b)). 
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Until recently, experimental data ons Ze -4) ebiaise 
excitation cross sections at low to intermediate energies 
have not been available; early measurements concentrated on 
the resonant charge exchange cross section (Fite, Smith and 
Stebbings(80) and McClure(63)). Lately, Reyeveae Morgan, 
Geddes and Gilbody(66) measured direct and charge-exchange 
excitation to H(2p) and H(2s) atomic levels in the range 
2-30 KeV. In addition to the above, H(2s) production cross 
sections have been measured by Bayfield(64) (from 3 Kev), 
Hill; Geddes and Gilbody(67) (from 1.9 KeV), and most 
recently, Morgan, Stone and Mayo(68) (from 2 KeV). Hill et 
al give somewhat larger values than those of Morgan et al 
below 3 Kev, but above this energy are in excellent 
agreement with them. On the other hand, Bayfield's much 
older H(2s) measurements disagree with these, showing a 
pronounced minimum near 3-5 KeV, though above 10 KeV his 
data agree with those of other authors. 

This major disagreement among the experiments is 
resolved by the rigorous theoretical study given here; our 
results are in excellent agreement with the more recent 


measurements by Morgan et al . 


2. He?*-H(1s) and H*-He*(1s) Collisions 
Again most theoretical studies have been done for the 
high energy region (E2100KeV) where atomic state methods are 


more appropriate; among the methods widely used are the 
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elose-couphing ‘method(76/77,78) 9 and unitarized distorted 
wave approximation (UDWA)(79). As expected, close agreement 
between experiments and the theoretical results is found at 
these high energies. 

Below 100 KeV ,however, the situation becomes quite 
different. Attempts to reproduce the experimental behavior 
using atomic state basis schemes fail miserably and 
different theories give widely divergent predictions. It 
therefore seems clear that a molecular state approach is 
essential for this energy region. 

The earliest molecular state treatment for He?*-H(1s) 
collisions was reported by Piacentini and Salin(70). Their 
3-state work was much more fully extended in a very thorough 
Study by Winter and Lane(71) using up to 22 basis states. 
Both these studies took the reference origin for electron 
coordinates on the proton but did not take ETF effects into 
account (hence they lack proper translational invariance). 
Since asymptotic boundary conditions for all states 
H(nl)-He?* are satisfied by this choice of reference origin, 
the total charge exchange cross section may in principle be 
calculated cciromitthe uidirecit mcross sections by using 
probability conservation; but individual state cross 
sections will not be correctly predicted, especially for 
charge exchange excitations. Even for the total cross 
sections, relatively poor agreement with experimental 
results is found when only a few basis states are used; 


agreement is improved if the basis set is extended greatly, 
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but convergence as function of basis size is slow. 

Recently, thorough calculations using PSS_ theory 
corrected with Bates-McCarroll (one-centre) ETF's have been 
done by Winter and Hatton(73) and Hatton,Winter and 
Lane(74). They examined convergence of cross sections with 
respect to basis size and calculated some individual 
He*(n=2) state cross sections, as well as the total charge 
exchange. Their results show very markedly improved 
convergence as function of basis size and are in generally 
good agreement with measurements on total and He*(2s) 
transfer cross sections. However, as pointed out previously, 
while some molecular states of HeH?* are nearly "atomic He*" 
States at most internuclear distances (e.g.,1s®@ ,2pTy ), 
others, such as 2p® _ , which correlates to H(1s), are truly 
molecular states with substantial two-centre character at 
the small R-values where coupling occurs, and 
Bates-McCarroll ETF's may not be adequate for these cases; 
molecular-state ETF's should be more appropriate. 

[Vaaben and Taulbjerg(72) have reported a molecular 
state calculation of total charge exchange cross sections in 
He?*-H(1s) collisions using a molecular state ETF based on a 
common iswitching function for all states. Their results are 
in very substantial disagreement with all the results cited 
above, particularly at the lower energies studied (their 
cross: section: at 1 KeV is only 50 % of the value reported by 
Winter et al); this is very surprising, in view of the fact 


that at lower energies the cross section is dominated by the 
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ZpDO Zp and’ “2p0-—"-s00 ‘couplings’ and’ fs relatively 
unaffected by ETF effects. In this thesis, I have done 
calculations uSing the switching function of Vaaben and 
Taulbyetcga andy rindy=s contrary “towytherr "reported values* 
values in good agreement with Winter et al at 1 Kev (see 
SectionV1I.B.2)]. 

Winter, Hatton ~and ° Lane(75)~ have carried out 
calculations of direct excitation and charge exchange in 
H*-He*(1s) collisions, again using Bates-McCarroll ETF's; 
only charge exchange into H(1s) is included. (They also 
checked the probability for capture to H(n=2) levels in a 
22-state PSS calculation by the method used by Winter and 
Lane(71). The same questions may be raised as_ before 
regarding the validity of Bates-McCarroll ETF's (versus a 
molecular ETF description based on switching functions)). 

In this thesis I have performed calculations for the 
HeH*?* collision systems using molecular States, with 
molecular state ETF's based on different switching functions 
for each molecular state. Convergence of the results as 
function “or *basis “size FP as® ‘well’ as®*behavior®*of both 
individual state and total cross sections, is studied. It is 
shown that individual state cross sections especially are 
sensitive to molecular effects on ETF's. Surprisingly, even 
for thes quatitatively “He*-like” 1s0- state, molecular 
character has a big influence on coupling matrix elements 
and this effect is particularly large for cross sections in 


H*-He*(1s) collisions (1s@ initial state). 
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Experimental meaSurements for He?*-H(1s) collisions 
have been done by Fite, Smith, and Stebbings(80) and Shah, 
and Gilbody(83), Nutt, McCullough, Brady, Shah, and 
Gilbody(84) and Bayfield and Khayrallah(82). All _ the 
measurements except those of Bayfield and Khayrallah are in 
reasonably good agreement above 10 KeV. (Bayfield et al give 
Somewhat larger values in this energy region). Below 8 KeV, 
only two experimental measurements are available (Fite et 
1(80) & Gilbody et al(83,84)); however, data of Fite et al 
lie about an order of magnitude above the more recent values 
of Gilbody et ai. For H*-He* (1s) collisions, the 
experimental values have been determined from the measured 
He*?* formation by Peart, Grey and Dolder(85) and Angel, 
Sewell, Dunn, and Gilbody(86,87). Dolder et al measured the 
He** formation cross section from c.m. energy 3 KeV_ and 
higher and their values tie in well to those of Gilbody et 
al (who measured only at higher energies 260 KevV.). These 
values for the He?* formation cross section include both 
ionization and charge exchange to both ground state H(1s) 
and excited states H(nl) and this fact should be taken into 
account when one compares the theoretical results with the 


measured values. 
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II. GENERAL THEORY 


The original perturbed stationary state(PSS) method and its 
serious deficiencies are described. The corrected PSS theory 
incorporating electron translation factors (ETF's), 
following Delos and Thorson(12) (see also Thorson and Delos 
11), is summarized. The resulting coupled equations (in the 
classical trajectory formulation) are presented in a_= form 
adequate for low energy collisions. Both differential and 


total cross section mfiormulae are shown and a brief 


discussion of the coherent effects is given. 


A. THE PSS THEORY 

In an atomic collision, if the relative velocity a haan 
of the collision is smaller than the orbital velocity of a 
bound electron oi 4, ce. << V. , the electron, might 
be supposed to adjust smoothly to changes in the nuclear 
position. In such a case, the molecular adiabatic states 
should form the best expansion basis to represent the 
electronic state of the colliding system: fixed nuclei 
Born-Oppenheimer (B-O) adiabatic states are chosen as the 
basis set and the nuclear motion itself may be regarded as 
the perturbation which causes electronic transition between 


states. This formulation was named the perturbed stationary 
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State (PSS) theory by Mott and Massey(2), who originally 
developed it. 

The conventional form of the PSS method can be 
developed either in a classical trajectory form or in a 


fully quantum mechanical form. 


1. Quantum Mechanical Form 
In a quantum mechanical version, the total scattering 


wave function can be expanded in B-O wavefunctions ahah ; 


yo= DX p(B) [n> II-1l 


where x, (FR) describes the nuclear motion associated with 
internal state n. The electronic state |n» isimeathe 
eigenfunction of the electronic Hamiltonian, at fixed 


nuclear configuration R, 


>> =) 
7R) n> = -c_ (R) ime ae 
fe edise ie See aieay 


where Heg iS the sum of potential energy v(r,R) and 
electronic kinetic energy Te . 
Both the eigenvector |n> and its eigenvalue depend 


Pparametrically on R. 
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The total wavefunction satisfies the time-independent 


Schroedinger equation for the entire system, 


2 aS) 
+ = ~ 
ea LEE EDAD eel, Ti=3 


Coupled equations for the components Pm) Of (an eabStract 
VeCtOrs \(R) ‘are © obtained by multiplying Eq. (11-3) from the 


left by bra €m| and have the form (see Ref.(11)) 
(ou) ah, + B(R)) x (®) + €(R)X(®) = EX(R) —II-4 


Since eGR) is a diagonal matrix in this adiabatic 
representation, all the couplings between different states 
are. caused by the, matrix Ps P(R) is called the nonadiabatic 
coupling matrix; it is responsible for all the electronic 


transitions and is defined as 


Pp (R) = ~ihcm|V | n> ti-5 
mn R 


2. Classical Trajectory Form 

Alternatively, in the classical trajectory version, one 
can assume that the nuclei A and B move along a classical 
path, specified by the internuclear vector R(t), and that 


the wavefunction for the electron satisfies a time-dependent 


&f& 


dy 
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Schroedinger equation, 


Bi te) eee) = iK3/att (r;t) t1-6 


If the total wavefunction tT is expanded in the set of B-O 


basis kets |n> , 


T = Dag t) [n> Se 


coupled equations for the components a, (t) of an abstract 
Vector a, Which result from substitution of Eq.(11-7):. into 


Eq. (11-6) and multiplying again by m | from the left, are 


> ae as 
tha = [e + VePla II-8 


where 


V = aR/dat 


The classical trajectory version of the theory may be 
derived as an approximation from the more rigorous quantum 
mechanical version, under certain conditions, which will be 


discussed in the following section. 
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3. The Validity of the Classical Trajectory Version of the 
Theory 

The problem of the derivation of the classical 
trajectory version of scattering theory from the more 
rigorous quantum-mechanical version, and the corresponding 
reduction of the coupled second-order equations (11-4) to 
first-order time-dependent classical trajectory equations in 
Bq. (r?=879; has been widely discussed in the 
literature(17,18,19) . The reduction depends on the use of a 
semiclassical approximation to describe the nuclear motion. 
A simple illustration of the methods involved is provided if 


we set 


y(R) = exp [inV-R/fi]a(R) hese 


ne ‘ : Sy AP 
where V is a constant velocity of magnitude (2E/u) 


(impact parameter approximation.). Applying Eq.(II-39) to 
Bo. -—4)) and aisorting: waccordang: = tomipowerstmof cumini owe 


immediately obtain, 


= = 
[-inveV + © + VeP]a(R) = Olu”) II-10 


Eq.(II-8) is recovered if the RHS of Eq.(II-10) is neglected 
and a(R) is interpreted as a(t) according to the classical 


connection R(t). More generally, the plane-wave factor {iin 
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Eq.(11-9) may be replaced by a suitable semiclassical 
approximation describing nuclear motion, in which V=dR/dt 
meed° not be constant(18). Conditions under which a 
semiclassical reduction to the classical trajectory version 
of the theory is valid vary somewhat according to the 
methods used, but always require that(18): 
(1).The momentum Ak = [2u(E-e (R)) 177? must be 
large,i.e. de Broglie wavelength for the state |n> must 
be small compared to distances over which € (rR) changes 
appreciably (this is the usual JWKB condition) 


(2). Differences ett San peeks , i1.€., momentum and 
n 


m 
energy transfer in a transition must be a small fraction 


of the total momentum and energy in collision. 


In most collisions of atoms at KeV or higher impact 
energies (and even lower fOr heavier atoms),such 
semiclassical approximation to the nuclear motion is 
certainly adequate, and in this study the classical 
trajectory form of the coupled equations is adopted for 
actual calculations. A separate question arises, whether a 
Straight-line (impact parameter) approximation is valid or 
other more sophisticated trajectories (Coulomb trajectory or 
other trajectories) are needed, and this point will be 
Ovecussed in-a later section(11.C.4). However, even for such 
high energies that the semiclassical approximation is valid 


and even that a straight-line trajectory assumption is 
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valid, a collision can still be slow in the sense that the 
collision speed is less than the orbital velocity of bound 
electrons, so that a molecular state description of the 


electron system 1S appropriate. 


B. DEFECTS OF PSS THEORY AND ORIGIN OF ELECTRON TRANSLATION 
FACTORS 
In both quantum mechanical and classical trajectory 
formulations of the PSS theory, transitions arise from the 
Same non-adiabatic coupling matrix,P, which is supposed to 
describe changes in the adiabatic state as the nuclear 
position alters. However, closer investigation reveals a 
fundamental ambiguity in the meaning of P: the derivative in 
Eq.(II-6) is taken while keeping the electron coordinate 
fixed with respect to its reference origin, but this means 
that -P depends upon the reference coordinates chosen to 
represent the system. Major defects of the PSS approximation 
arise from this situation. We can explain the problem as 
follows) 129 : 
(1). In the usual treatment of molecular systems, the 
centre-of-mass of the nuclei (CMN) is taken as reference 
Origan fori@they electron ‘system? Ate large internuclear 
separations, however, the physically appropriate origin 
for an electron coordinate is not the CMN but the atomic 


centre to which the electron is bound: the electron 
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moves with the nucleus to which it is attached. 

(2) i P is assumed to represent effects of polarization, 
GYVStoOrr1on change-or-character, etc, in the: molecular 
wavefunctions as they adjust to the changing molecular 
fields, but if the reference origin for the electron is 
taken (for example) as the CMN, then a contribution to P 
also arises from simple displacement of the basis 
functions with the moving nuclei to which they are 
attached. This introduces spurious couplings into the 
EaMabul xe scence melt 6S ObViOUS thatethe: part jof P which 
merely represents displacement of basis functions with 
respect to electronic reference origin is not a real 
coupling which could lead to electronic transitions; the 
"displacement part" must somehow be removed from P, to 
remove these fictitious couplings. 

(3). Since P is origin-dependent, the PSS theory as it 
stands does not meet physical invariance requirements. 
An electron bound to a particular nucleus has 
translational motion, with respect to some reference 
origin such as CMN, due to the nuclear motion with 
respect to that origin. Transformation theory requires 
that .when a state which is at rest in one reference 
frame is described in a new frame, with respect to which 
it is uniformly moving with velocity a then the 
wavefunction (as seen in the new frame) must be 
Subjected to a unitary transformation to maintain 


Galilean invariance: 
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v (r/R) = Ud (r,R) II-1l 


where U is a unitary transformation operator with the 


form 


U = exp[(im/h) (Wer - 5V°t)] ragealy 


This operator is here conventionally referred to as an 
electron “translationz=factor (ETF): it describes the 
momentum and kinetic energy of the electron, associated 
with its translation along with a moving nucleus to 
which it is bound. Because PSS theory ignores’ such 
factors, ie does not meet Galilean invariance 
requirements, and in particular the individual molecular 
States which make up the PSS expansion cannot meet 
asymptotic boundary conditions for scattering and are 
not in one-to-one correspondence with correct initial 


and final scattering states. 


The defects of PSS theory which result from ignoring 
these considerations are not purely formal ones, but may 


lead to grossly incorrect physical predictions in some 
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cases. The nonadiabatic coupling matrix,P, has the following 
unphysical properties: 
(1). For certain discrete-discrete transitions, P tends 
to a finite value as R-9c00o (e.g., H,2* 1803 —2sg5 
coupling Figure(II-1)), and the magnitude in the small R 
region 1s unreasonably large (e.g., H,* IsQZ ——3dgq 
coupling Figure(II-2))(see also Ref.(38)). 
(2). For ionizing transitions, the P-matrix elements are 
physically unrealistic: each initial bound electronic 
state 1S coupled to 30-40 partial waves in the 
electronic continuum and the couplings extend out to 
extremely large distances (750 a.u.)(see Rankin(29b), 


Ranken and Thorson \(29a)); 


It might be argued formally that, in spite of these 
defects, since the molecular states {|n>} (including the 
continuum) form a complete set, the resulting coupled 
equations and hence PSS theory are still completely rigorous 
as a scheme for solving scattering problems. This would be 
true, provided that: 

(1). a complete expansion basis were actually used 

(including continuum states); and 

(2). the resulting wave function were projected onto a 


rigorously correct set of final scattering states. 


However, such a procedure is completely impractical, and 
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Isqj—> 2so, RADIAL COUPLING 


COUPLING ( ALU.) 


(P +A) Matrix 


Figure(II-1) Nonadiabatic coupling matrix, P (PSS matrix) 
and ETF (see text) corrected matrix, (P + A) for 1s0, -2s05 


coupling matrix element. 
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Isog>3do, RADIAL COUPLING 


P- Matrix 


0.184 


COUPLING ( ALU) 
0.062 


(P+A) Matrix 


0.016 


Figure(II-2) Nonadiabatic coupling matrix, P (PSS matrix) 
and ETF (see text) corrected matrix, (P + A) for 1sq -3dq5 


coupling matrix element. 
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when the expansion basis used is truncated to a finite set, 
the defects of PSS theory can lead to serious errors. 
Defects of stiesP mmatrix jall eresult aefirom ) the: + :spur ious 
"displacement" terms described earlyer.acTo srobtaimora 
physically correct theory, it will be necessary to correct 
the nonadiabatic coupling matrix by removing such terms. In 
the next section, we shall see that when a description of 
electrons translation factors «sis, ‘properly included) in: the 
theory, the most amportant. result »:is, that .the«-matrix P 
appearing in the coupled equations is in fact replaced by a 


corrected coupling matrix, PtA. 


Bates and McCarroll(4) were the first to recognize the 
serious defects of the PSS method due to its neglect of 
electron translation effects, and they developed a 
modification of the theory which at least removes’ the 
asymptotic defects. In their formulation, a plane-wave ETF 
of the form (II-12) is associated with each molecular state. 
Assuming that the molecular state |ny is asymptotically 
correlated in unique fashion to a particular atomic state on 
nucleus J, they associate with state n a factor (11-12), in 
which the velocity V; which appears is the velocity of 
nucleus J with respect to (say) the CMN. This guarantees 
that the basis states in the expansion are in one-to-one 
correspondence with the correct atomic scattering states, 
and removes such defects as the spurious finite values as R 


—+>o for some P matrix elements, and much of the very 
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long-range behaviour of the ionizing-transition matrix 
elements. (In the special case of the symmetric molecular 
systems (for example,H.* ), where molecular states have 
gerade (g) or ungerade (u) symmetry and do not correspond 
uniquely to atomic states on a single atom, a more 
complicated procedure is used by Bates and McCarroll(4): 
they combine pairs of g and wu molecular states to form 
fragments which do correspond asymptotically to states on 
unique atoms ,J, and then associate an ETF of form (II-12) 
to each such fragment.) 

In spite of these major improvements, the formulation 
of Bates and McCarroll(4) still contains assumptions which 
are physically incorrect, because they: assign to an electron 
in a molecular state an ETF which describes transport 
associated with a particular atomic centre. At the finite 
internuclear distances where real couplings OCGUr, a 
molecular electron may be partially attached to both centres 
(i.e., a molecular state may have two-centre character). For 
Such situations, a proper ETF for a molecular state must 
itself have molecular, two-centre character. AS we will see 
in the next section, this can be done using the device of 


Switching functions. 


he 


ci oie had Sievert etka vend f 


elucainn, of v2 snes dh ae 


yr sesec2 SAS a. On ? fone. ty 
—wogasixos idf ob mE 7 Taam 
- pa ceate't 
i 8 ee eee.’ = ba 8 No 
2 | i. OM bi LE kf i is 
eel 
+ eves raelusetom(y baw’ ih ey ‘exis 


i eedere of yitetiyoagedae  wipaicesitae hie | 
-fety n2ei. 3 : {os a8 a3ef mre wae bpsaihs see : 
| fe - 


Ltetwmioid. ads BS a aN dh sot ear obs 3 te si kis v | 
” ; 
fiw grersghuqes anes tno athe. on Etqnatidete: “t 


ne dogiles me oF ips Bat Yond seat Uy 
sadastd aetlasdsb, geide yee ne iad 
sagen 
ative? #2 3A sapsiina atone satvala veg 5 Biv 8 ai 
Ae 7 
5 vere A . epi) f, bea Teor lator: shsraceet 4 7 
Serined died or hedsagine. at igitian’ oo te: poauate a 
10% (ae 7Ob elo gagoe spies, ovate ys: rapstom. 
{cum AIH12 tuirGalon S 258, so 3Hapt & etna 
see (ddu one 2A assosted 3 (dss ov atu sa bom 
$0 es ivea S07 grifay srolk ait nap ae Prt: 1388 


: : ano bese | 


” 
ae 
“< a F _ 
s 
. as + : PT 7 
; i : biel +P 
/ ’ 5 
4 ) 


35 


C. THE CORRECTED PSS THEORY (MODIFIED By ELECTRON 
TRANSLATION FACTORS (ETF) ) 

We now present an account of a modified PSS _ theory 
which is corrected to include electron translation factors 
METS) “and their effects. The ‘formulation meets all 
requirements of Galilean invariance and asymptotic boundary 
conditions. The -method used ~ is-— that First presented 
rigorously by Thorson and Delos(11), who showed that ETF's 
appropriate for molecular states can be constructed using 
Switching functions. In effect, these describe an "electron 
translation” which is a function of the electron's local 
behavior: alternatively they can be viewed as defining a 
local reference origin for an electron adapted to its 
"degree of attachment" to one or the other of the nuclei. 
Asymptotically, provided the switching Bune tion meets 
certain boundary conditions, a description effectively 
identical to that of Bates & McCarroll is recovered as a 
limiting form. 

The original formulation of Thorson and Delos(11) was 
based on the assumption that a single switching function may 
be used for all electronic states of a given system. 
However, study of the shapes of electronic wavefunctions 
reveals clearly that some states show substantial two-centre 
character at finite internuclear separations, while other 
groups of states have almost entirely one-centre character 
over nearly the whole range of internuclear separations. 


This as well as other evidence strongly supports the 


; 
1 


aise: vieuorts aime oe — ae 


©. - = ; - ; 4 - 


worries te cylkKdOM) ae 


‘ . wh! i 
wiaett 223 Osnsifogm Go hull 


amy | treabete Te sco Ss Me 


a) hy 
(1) aekett mam ong gailowt,. 2 
solau Harwictteags 0 bead baht, 684 
theta” “eg atti toeab aeany $hntXe! ot “ans gt, iS 
fuool a pontoels « or io Ho old oF et ake, be 
ty initeh sé Bewaty 6A es F aan tdseitooaneah 
Hoi an Bapqeba!, Howes ity fa es ip Het 9s mee a ny 


Y dnp. 2 r ee | sire! ge ' 


= 
ar | 


> oy) Wak tut 
sae A Lat AM F 7 iJ a eos of 7 a on ae a apie 


seen  cotaorwt ~ petting si ‘act Reser, ‘gies ! 
vleovigotiie nos dgaae yen A at Eins igi be ai 
n 28 Seeeveuss~. at oa s tose ta ros “8: 
asy (Ts bento ons aun 9° it rmtnrt tanieine 
you qo sdoned pnia@rsive nisAPe 4 Jas cotaguuées om 
Te Iays er & to Peey piderdoele ite. 4% : 
edelgomrinavew tinorchata” to. saaade eit> 26 hus a 
eranes~ ast lebtnetanee voz eeipze sue ted? Yio 
wise ed tity , nl I6> KEg2 rexbsunteyai Vent 8. 


_SetaeetS «6attieo*uno | uiasiine Jecwis sl aege | 
ane tYaragqes aol see want Wo. aptea ‘shady oe 


ny tT 
Le 

VF 

7 


g * ; 1 


36 


assumption that different switching functions for each 
molecular state are more appropriate. A formulation based on 
this more general assumption has been given by Delos and 
Thorson(12) in the framework of classical trajectory 


collision theory. 


This theory 1S reviewed in slightly modified form for 
practical use bearing on the processes eccurring in 
(symmetric) HoH (1s) and (asymmetric) He?*-H(1s) and 


H*-He’(is) collisions. 


1. Coupled Equations for Slow Collsions 

It is assumed that the nuclei A and B move along a 
classical path, specified by the internuclear vector RCL). 
and that the wavefunction for the electron satisfies a 


time-dependent Schroedinger equation 


> > => > 
iho /otT (r;t) = Ho (erR(t))T (xe) ti 16 


The total electronic wavefunction Y 1s expanded in the ETF 


modified set in the form 


T= ae a (OE ¢. 11-14 


where db are the basis states describing internal motion of 
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the electron and F, are associated electron translation 
factors (ETF) describing the electron's transport momentum 
ano kineticvenergysas it 1S carried alongmin state “do. with 
one nucleus or the other. If the basis & were a set of 
atomic states, then the appropriate F, 's would be 
Bates-McCarroll factors of the form (I1I-12), but since 

d, are molecular states (for a PSS expansion) we shall 
Gefine the ETF's F,in a somewhat different way. 

For convenience, we take the geometric centre of the 
system (GC), rather than CMN, as the reference origin for 
electron coordinate r. AS Delos and Thorson(12) point out, 
the physical description is in any case invariant to choice 
of reference origin, if Galilean invariance requirements are 
met correctly for that choice; furthermore, since the 
velocity of the CMN with respect to the geometric centre is 
simply -tXR »where ~ ® (Ma —Ms) /(Ha + Me) , it is easy 
to construct ETF's with respect to CMN as reference origin, 
if these were needed, by the rules of transformation theory 
eSchmid (27) )-. 


With GC as reference origin, we define the ETF's as 
8 


Ee. = exp [ (im/h) (w +r - i reeds) A II-15 


— . ° 
where w, is defined as 
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and 
<> 
Vere R = aR/at 
Wa is effectively a local transport velocity for an 


electron in molecular state n, and it is specified by the 
Switching function f,» (#,R). Since the electron in a bound 
State is attached asymptotically to one of the nuclei, we 


have the asymptotic constraints 


: EARS _- > a be a oes 
LiMo Rh) = ly, ir eed = |r = =R| ranite A eae 
n B 2 
Ro 
p > =e Ie 5 : 
i ad a We + 5R| finite 


The atomic-type ETF's of Bates & McCarroll,(II-12), are then 
effectively recovered as R-— ,ocoo , since a "B-atom" state 
projects out a velocity W =+ 5 V, and an "A-atom" state a 
velocity W =- > V, in Eqmilil=16): 

The idea of using a switching function to achieve’ such 
a molecular ETF description was first introduced by 
Schneiderman and Russek(9), who also considered the 
constraints which must be met by such a function. Thorson 
and Delos(11) carried out the first rigorous formulation of 


a scattering theory using the device. 
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Apart from the asymptotic constraints given above, the 
foumenot the switching function ££. (t,R) is “not fully 
specified. The determination of appropriate switching 
functions for molecular states is an important practical 
problem, because the corrected nonadiabatic coupling matrix 
elements which result are in some cases extremely sensitive 
to the choice used for fm. Detailed methods for determining 
Switching functions are discussed in Chapter III. 

In this formulation, the expansion basis states , 
are again just the Born-Oppenheimer adiabatic eigenfunctions 
For the electronic Hamiltonian, ,which obey Eq.(I1I-2). To 
obtain the coupled equations, we substitute Eq.(II-14) into 
Schroedinger equation Eq.(II-13), multiply on the left by Fe 


k 
, and close with the bra Sal. with the result 


* * 
) <m|F_F |n>iha/at Bek ) tn Ea ~ ihV-V_) [n> 


+ <m|F_{ [H ]) - ifid/ot F } [n> a 
m n n 


F 
ein 


II-18 


These equations are equivalent to the Schroedinger 
equation(II-13) and, within the approximations in the 
classical trajectory theory, they are exact. After 


simplifying by neglecting small terms related to aV/dt and 
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> 


> 
V Wf , the coupled equations have the form (in the 


abstract vector notation) 


> > > 
s(v)#ha/dta = [h(V) + Ve (P(V) + A(v))la II-19 
where 
* 
S(V)= «mir le n> II-200 
mn INn 
* 
h(v) = <m|F_F_h|n> II-20b 
i AY) | II-20 
= < -j > - Ce 
dla <m|F OF wi oe n 
* 
ae hee: > II- 20d 
Ee (im/h)<m|F FL Sj |n ° 
ere eye II-20e 
n rag 4) 


For low and intermediate energy collisions, which are of 
particular interest in our present studies, we can apply a 
Substantial simplification to the above coupled equations . 
Expanding the ETF's in power series of velocity,V, and 
neglecting all terms of order V* and higher, we can obtain 
Simpler coupled equations by retaining terms to first-order 


in V; and to second-order in V, as required. 
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Coupled equations to first-order in velocity 


When terms of order (Vv?) are also neglected, the 


coupled equations have the simple form 


tha/at a= (€ + Ve(B + A))a II-2\ 
where 

oar <m|{h|n> Smn II- 226 
eo 
p= <m| (-iV_) [n> it>22b 
“mn R 

= (ih/m)<m|[h,2_] |n> II-22c 
mela: n 


In these equations ¢ and P are familiar terms 
which arise in the PSS theory in which ETF's are 
neglected. P is the nonadiabatic PSS coupling term. 
The new coupling matrix A arises from the presence 
of the ETF in expansion(I1I-14); it is the correction 
matrix which identifies and removes’ from P that 
portion which represents just simple "displacement" 
of basis states with nuclei to which they are 
attached. Eqs.(II-21) have just the same form as the 
OFiginaly eqbations¥!(11-B8)avof ePSSntheorysi but? the 


nonadiabatic coupling matrix has been corrected to 


remove the fictitious couplings due to displacement. 
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Inclusion of second-order terms in velocity 


When second-order terms are also included, the 


coupled equations are 


=> 


tha/at a = POS ear 
/ eee dv (P + Aja, 


+2, LamAngcx| 2 )-2 0100 ah, + amy tH 8 1) > 
< y <k| [V> (S475) ) |3><3|V> (2K, + (im/h) ie al |n>} 


Heh 


2 >2 >> 20> > 2 
(m/8) [(£)-1)V POY 22) (VEY) 


Ao) ray cae eae ee SAO aes 
r)(£ VeV_ £ + 20-¥, Ila, II- 23 


These equations contain the same first-order terms 
as (II-21), but also include the very messy and 
complicated second-order terms shown. These terms 
become increaSingly important as_ the collision 
velocity becomes comparable with the orbital speed 
of a bound electron, but aS a rule we have_ found 
that for energies below about 5 KeV/amu their 
effects in the systems studied are negligibly small. 
Their importance, however, also varies somewhat from 
system to system; if one considers couplings 
involving states which have strongly "molecular" 
character out to relatively large internuclear 


separation, differences between switching functions 
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for different states have a relatively much larger 
effect; for example, in the symmetric H*-H(1s) 
collision system, this effect is much more 
pronounced than for He?*-H(1s) or H*-He*(1s), where 
many of the states are effectively He*(nl) atomic 
states at relatively small internuclear separations. 
As a result, in H*-H(1s) collisions , second-order 
terms have a noticeable effect already at collision 
energies 4-5 KevV, while for the HeH?* system they 


remain quite small even up to 20 Kev. 


These are the required forms for the coupled equations 
for practical applications to slow and intermediate energy 
collisions. They can be solved numerically. Because of the 
neglected terms in (dV/dt) and the higher order terms in ye 
Bont 21) and ~ Eq. (ll-—23) may not exactly conserve 
probability in general, but this failure of unitarity is 
small, of the order of the neglected terms, and we will 


prove this in the next section. 


2. Unitarity 

As pointed out earlier, the coupled equations (II-19) 
are equivalent to the time-dependent Schroedinger equation 
(II-13) from which they were obtained, and, within the 
classical trajectory description they are exact. Since the 


Schroedinger equation conserves probability, it follows that 
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the quantity 


* 4 
<T|T> = a (ta, (t) <m|F_F_|n> ear yn 
m,n mn 


is conserved, i.e., 


* e ox 
thayat<T|T> = ) fas. (thal - [-iha Isa, 
m,n 


+ ae hoy oe Syeiavii= 10 II- 24b 
m mn nh 


If (as we have done) the exact equations (11-19) are 
approximated by neglecting certain terms, then there will be 
a corresponding failure of the conservation of probability 
Poe ha. Gli=24b)% 

To derive Eqs.(II-21) or (11-23), we have (1) neglected 
terms in (dV/dt) and Ose and (2) after expanding ETF's 
in powers of the velocity, we have kept terms in the coupled 
equations only up through first-order in velocity 
(Eq. (11-21)) or second-order in: velocity (Eq.(II-23)) and 
have discarded higher terms. Correspondingly, it can be 
shown that errors in the conservation of probability are of 
the order of the discarded terms. Thus, for the first-order 
eduations (11-21). it .can be shown explicitly that 
Eq.(II1-24b) is satisfied to within errors of order V?(and, 


of course, dV/dt and V-Vi£. ). It 1s, illuminating st oma cave 
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this proof explicitly (Green(14)). Using Eqs.(II-21) in 


eqs. (II-24b), we obtain 


, ki * > * 
thd/at<T|T> = a [V+ (P-P) + Ve (a-A*)Ja 


re | = 
+ Va (OS his hs) (3) a 


* ee D 
- crete mde VeV_S TOMY §))) ai II-2$ 


Using the relations 
* -1 -1 
A-A =-(im/h)[os h- hs 6] TI-26 


we find 


* * 2 
tha/at<t | T> = a [V+ (P-P Nae hv -V_S tC (Ve) arene ca ot 


and since it is easily shown that 


* = * 
-ihVs =P - P II-28 


it follows that Eq.(1I-24b) is satisfied to within errors of 
order V2. This proof is relevant to a particular point in 


the practical calculations. Since the switching functions 
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are different for different molecular states, it follows 
that Mer is not equal to Awn , and hence the coupling 
matrices are not Hermitian. Some workers who have included 
treatment of ETF effects in calculations have incorrectly 
assumed that these distinct forward and reverse couplings 
should be replaced by their Hermitian averages 
(e.g.,Crothers and Hughes(30)). However, as Eq.(11-24) 
shows, what is to be conserved is not ( aren) eel buteael 

<T|T> ) given by (11-24a) (to errors O(V?)), and as we 
have just shown using Eqs.(II-27), the non-Hermitian 
character of the matrix elements is just what is needed to 
achieve this. Similarly, one can show that if terms up to 
order V? are included as in Eqs.(II-24), then errors in 
probability conservation in Eq.(11I-24b) will be of order Vv? 
(as before, there are also errors’ of order (aV/dt) and 
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3. Components of P- and A-Matrices in the Rotating Molecular 


Frame 

In the foregoing developments we have _ shown the 
nonadiabatic coupling matrices (P+A) as vectors in 
three-dimensional space: 
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Now, it is well known that the Born-Oppenheimer eigenstates 
[n> depend parametrically only upon the earecuee of R, and 
not upon its direction, and it might therefore appear that 
the only non-vanishing component of P would be the radial 
one. This is not correct, because the gradient in Eq.(1I-29) 
is to be taken keeping the electron coordinate r fixed in a 
non-rotating reference frame; but the Born-Oppenheimer 
eigenstates depend only upon R, provided that the electron 
coordinates are described inj) they rotating) molecular 
reference frame. Let TH=(epyy 2) be the coordinates of the 
electron in the lab. frame, while r'=(x',y',z') denotes the 


Same point in the molecular frame. These are related 


By eos 6) 
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where @ @$ are the polar coordinates of the vector R. 


From this it can be shown that 
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where fis, ate ee , are the components of the electronic 
Orbital angular momentum operator in the molecular frame ( 
with respect to GC as origin). It follows that P has angular 
as well as radial components. 

inva classical trajectory description of a collision, 

the vector R(t) moves ina plane, which we may take to be 
$ =0. The velocity V=dR/dt has components: Ve =R, 

V. = RO . Corresponding to these components 


the vector P has components 


Rier.: > II- 3a 
eee -ith<k| (8/8R) preg |n 
and 
=] 
PE = -R <k |Z, [n> : II- 33b 


The corresponding components of the ETF correction matrix A 


are: 
aR = (im/ah) (e, - €)<k|£,2" [n> II- 34a 
' II- 4b 
AG =ten/ehrts, -ee,,)<klgae 12 


In the developments which follow, we shall understand 
Y(x,y,z) to refer to rotating molecular frame electron 
coordinates, i.e., we shall drop the prime Superscript in 
the above equations. In the molecular frame coordinates, not 
only the eigenstates n> ,but also the switching functions 
f(t;R) depend parametrically only upon R; therefore it 


follows from the above equations that both radial and 
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angular COuUDLING  samatrix elementsealhile3 3) nbLi234) 4 are 

functions only of R. The Born-Oppenheimer eigenstates are 
Z “~ 

also eigenfunctions of Lg , the component of orbital angular 


momentum on the molecular axis: 


> 
i n m “i | n> II- 35 


; R 
It follows that radial couplings (ps4 A ee may connect only 
States with the same m, while angular couplings connect only 


States for which es ua lee 


4. Choice of Trajectory 

In the derivation of a classical trajectory formulation 
of a collision problem from the fully quantum-mechanical 
description, the definition of the classical trajectory R(t) 
to be used is not completely arbitrary; it must be 
determined in some way which is self-consistent with the 
assumptions used in the derivation. In most 
caseS---particularly where the JWKB approximation is the 
basis for the reduction---the trajectory R(t) is defined by 
defining some suitable average momentum (in terms of the 
classical momenta associated with different states, |n> ). 

A great deal of discussion in connection with 
semiclassical reductions has concentrated on the problem of 


determining the best average trajectory. While it can be 
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shown in many cases that some choices appear to give 
somewhat more accurate results than others, this question 
remains still open. Delos and Thorson(18) state that if it 
makes very much difference which average trajectory must be 
used, then the validity of the semiclassical reduction 
itself may be questionable. 

For the collision energies considered in this study, 
the validity of the semiclassical reduction and therefore of 
the classical trajectory approach is not in question. For 
example, studies in the H*-H(1is) collision problem by 
Knudson and Thorson(42) show that the classical trajectory 
equations are able to reproduce fully quantum-mechanical 
calculations for c.m. collision energies as low as 30-50 ev, 
which is more than an order of magnitude below the lowest 
energies considered in the present studies. 

However, Knudson and Thorson(42) also showed that the 
classical trajectory equations must be solved using a 
classical trajectory which iS appropriate to the actual 
behavior of the system. In particular, for the case of 
strong rotational coupling in the united atom limit (R—> 0) 
between 2p0>5 and 2PM molecular states of the H’*-H(1s) 
system, it is important to use a classical trajectory which 
correctly describes how the system passes through the strong 
coupling region near R—»O. In earlier studies of atomic 
collision processes, it was a common practice to use the 
classical trajectory equations with the arbitrary assumption 


of a straight-line, constant-velocity trajectory: 


> > > 
R(t) = 9p + vt II-36 
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(this is often called the impact-parameter approximation). 
If the collision energy is high enough, this assumption can 
always be justified in the limiting case. However, Knudson 
and Thorson(42) found that as R-—> 0 the distortion of 
trajectory from a straight line due to Coulomb repulsion of 
the nuclei leads to major differences between the impact 
parameter results and the correct behavior at lower 
energies. Correct results are obtained for E2100 ev, if 
Coulomb trajectories are used instead. Additional 
modifications of trajectory due to electronic potential 
curve details were found to be much less important.Since the 
united-atom rotational coupling is known to be ae =™main 
excitation mechanism in many types of atomic collision 
problems (including cases we study here), the question of 
accuracy of impact parameter (straight-line) trajectories, 
versus Coulombic trajectories, 1S perhaps important at the 
lower energies considered (<1 Kev). 

We have tested the validity of the straight-line 
trajectory by comparing results with those found for 
CoulombicYtragectories in yH*-H(1s) collisions at 1. Kev 
(projectile energy). As shown by Knudson and Thorson, the 
maauneetfect. of ithe Coulombic: trajectory 1s to shift. the 
position of the peak of maximum probability to smaller 
impact parameters and create a new peak at very small impact 
parameter,’ which contributes. avery small effect to the 
cross sections at higher energies. In our studies at 1 Kev 


we found that the total deviation between excitation and 
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charge exchange cross sections computed with straight-line 
and with Coulombic trajectories is at most 5 %. Needless to 
Say, as the collision energy increases, the results will 
coincide still more closely. 

For the HeH?* collision system, Winter and Lane(71) 
made a careful study of differences between straight-line 
and Coulombic trajectory cross sections at 1 KeV, and found 
they were at most 8 %. 

Therefore we concluded that the straight-line 
trajectory could be used safely in these studies and we have 
Carried out our calculations on this assumption because of 
its Simplicity. For this case, the radial and angular 


velocities can be written, 
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where f is the impact parameter and V is the (constant) 


relative velocity. 


D. CROSS SECTION FORMULAS 
1. Total and Detailed Cross Section Formulas 
The coupled equations, Eqs.(11-21) and (11-23) should 


be solved to obtain the amplitude and probability as t-—»+@ 
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. These equations are solved numerically with the following 
initial condition: if the index "1" designates the initial 
state (1S05 and 2pqQ> for H,* system, and 2paq for the 
He?*-H(1s) , 1s0 for He*(1s)-H* system), the initial 


conditions for the coupled equation are 
a,(- 99) = V7 Cee Ou tH ss II-38a 
bor gHeH*" II-38b 


and the final state amplitudes af+o,f ) for given energy 
and impact parameter are computed. The probability of 


excitation to molecular state k is defined as 
P(E oO = la (+ o;0)|° II-39 
k k : 


The corresponding integrated cross section for molecular 


State k is conventionally given 


eo 
Qy = 2x [oas P(E,p) II-40 
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However, tO compute probabilities and cross sections for 
excitation to specific atomic States(j), atomic state 
amplitudes byt (+0> 5 PF ) must first be formed by an 
appropriate coherent addition of molecular state amplitudes 


before using Eq.(II-40),i.e., 


Bois sores as) Tee eantn 0;0)exp[-(1/h) Ce, (t le tase ] 


II-41 


where Y is an arbitrary phase, and €. is the atomic state 


energy. 


The coherences considered in Eq.(II-41) are of two 
types(a,b) for the H,* case and one type(b) for the HeH?* 
case. (In the asymmetric system, the MO correlates in 
one-to-one fashion to a corresponding AO.) 

a. g and u Coherences 

The direct and charge exchange amplitudes for 
H,* system are formed by addition and/or subtraction 
of the amplitude for matching pairs of g and u 
molecular states. For the incident channel where 
final amplitudes are significant for both g and u 
components, the well known resonant charge exchange 
oscillations result(e.g. resonant H(1s) charge 
exchange involves both the initial states 1s and 


2P05 states). For the excited channels , on _ the 
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other hand, amplitudes for excitation of g molecular 
States are smaller below 5 KeV, with the result that 
the direct and charge exchange cross sections are 
nearly the same(3p% and 4f£G for u-State and 2s 


5 and 3d9 for g-state). 


b. Hybridization Coherences 


In a hydrogenic ion or atom where atomic levels 
like 2s and 2p are essentially degenerate, the 
molecular states are asymptotically correlated not 
to atomic eigenstates but to the hybrid (Stark 
field)is'States.(@For HH, , ‘the pairs of molecular 
states (2505 » 3005 ands A3od> rAtOs ) 
COLrelatewstOn OG) Lande tue sp  tybrids” ef “H(n=2),, 
respectively. For n=2 level in HeH?*, the molecular 
States (2se , 3007 ) correlate to the sp-hybrid 
of He*(n=2) atomic states. These mixings must. be 
unscrambled to obtain 2S and 2p atomic amplitudes, 
and cross sections Lor H(2s) ,H(2p) and 


He*(2s),He*(2p) individual atomic states. 


2. Differential Cross Section Formula 
A formula for the differential cross section is given 
by McCarroll and Salin(20) within the frame of the Eikonal 


approximation, 
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where p is the reduced mass, Vv is the collision velocity, 
6 is the scattering angle , Jm(x) is the Bessel function 
of order m, and mis the difference between the magnetic 
quantum numbers in the initial and final molecular states. 
The total exchange probability may be therefore calculated 


from the expression 


eee meee 
ibe 
ih, Ey Bits saan ac 11-43 


=x P, (oS (ey +00 (8) ) 


which describes the ratio of differential cross sections, 
and d,c denote direct excitation and charge-exchange 
excitation, respectively. Eq.(I1-42) has been derived under 
the assumption that the forward scattering which occurs at 
the angle 6 -0 is the dominant process in the scattering. 
[McCarroll et al used this assumption and expanded the terms 
which involved the angle, 0 , and then retained the first 
term of the expansion only.] Therefore this formula is valid 


only for small angle scattering. 
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III. SWITCHING FUNCTION AND THE COUPLING MATRIX 


A. INTRODUCTION AND SUMMARY OF OTHER WORK 

In this work, we have used the switching function, 
f£(r;R), as a device to represent the translational motion of 
an electron in a molecular electronic state, as a local 
function of its position. In a molecular state, it is not 
possible to Savea thet the elect romwaisSe attached» to aa 
particular nucleus (as would be true for an atomic. state); 
in effect, the switching function shows the local degree of 
attachment to either nucleus. Asymptotically (R-— oo), we 
can Say that a bound state electron must be attached to a 
Gefinite nucleus, and therefore it follows that there are 


asymptotic constraints 


lim [£(%;R)] = +1 ¥n = |f-R/2| finite III-la 
Roo E 

lim [£(F;R)] = -1 r= \r+R/2| finite et tT—1h 
Roo A , 


In, the, united. atom limit.R.=7 0, 2t.can be ashowna that 
the bound electron should translate with the centre of 
charge of the system (rather than the geometric centre). For 
an asymmetric system, the switching function (which 
describes translation with respect to the geometric centre) 
must therefore tend to a finite mean value, not to zero (as 
has sometimes been assumed in earlier work). This has_ been 


pointed out earlier by T. A. Green(32-a) and more recently 
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has been emphasized by Taulbjerg and Vaaben(28) in a_ recent 
paper on switching functions. 

Apart from these two types of limiting constraints, 
however, the specific form of the switching function is not 
defined. 

The formal requirements of a scattering theory are met 
using ETF’s based- on an arbitrary choice of switching 
function, provided it meets the asymptotic constraints 
(IIlI-1a,b). The results of a scattering calculation should 
thus be invariant to the choice actually made for f, 
assuming a complete basis set is used. However, practically, 
the basis sets used for close-coupling calculations are far 
from complete, and one then finds that the choice made _ for 
the switching function directly affects the coupling matrix 
elements and through them the resulting cross sections. 

This problem has great theoretical importance because 
the corrected nonadiabatic couplings for many transitions 
are in fact extremely sensitive to the form and parameters 
used for the switching function, and a precise determination 
See swicching functions formmolecular States of HH. “and Hen” 
waS an important investigation topic for this research. 

Var1ous © "authors (21-33) have” proposed” methods” for 
treating ETF descriptions and, in particular, fOr 
constructing switching functions. A good discussion of 
earlier aspects of the ETF description problem has_ been 


given by Riera and Salin(34). 
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Apart from the methods we have used ifs this work, which 
are described more fully below, most attempts to. find 
Criteria for determining switching functions have taken the 
Euler-Lagrange variation method(21a) as the formal starting- 


point for obtaining such criteria. 


A first application of the Euler-Lagrange method to the 
determination of ETF's was given by Riley and Green(21b). 
However, in that paper they considered only 
Single-centre-type ETF's, i.e., they considered only ETF's 
based on R-dependent parameters and did not treat the more 
general problem of ETF's based on switching functions. 
Recently, a much more elaborate formulation, including 
Switching functions aS a possible device, has been presented 
byey THO AvesGreen(13)% ‘funmekthasintormulation) ve" optimum" 
Switching functions may be determined variationally by 
solving a set of coupled differential equations. However, 
these equations are extremely complicated and their solution 
seemed to us to be quite impractical for purposes of our 
present work. 

More limited application of the Euler-Lagrange method 
to ETF description in the case of a single ‘molecular state 
leads to equations closely resembling the equation of 
continuity for the probability density. Some authors have 


attempted to use such equations, or equations closely 


‘We thank Dr. Green for showing us these equations in 
advance of publication. 
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related to them and based on physical arguments, to derive 
Switching functions. One such attempt was made formally by 
Schmid(27); it can be shown, however, that his "definition" 
is not unique, since it requires the evaluation of a line 
integral whose value depends on the path used to compute 
ite e i, s the quantity integrated is not an exact 
differential. Also, zeroes in the molecular state wave 
function appear to introduce poles and essential 
Singularities into Schmid's definition. 

A more successful attempt has been made by Ponce(31), 
who also applied his definition in a very approximate 
fashnuon to. calculations .on the 1s% ,2p0y and 2p7t states of 
H,*. He obtains "velocity field" functions which may be 
roughly compared with our switching functions, and in fact 
they do resemble our functions in a qualitative (but not 
quantitative) way. Difficulties with Ponce's method are that 
et uses Cartesian separations of electron coordinates (the 
molecular states are separable only in prolate spheroidal 
coordinates). 

Crothers and Hughes(30) have used the Euler-Lagrange 
scheme in a still more limited context, to define ETF's 
whichdanegequivaléntes tosolowesty order ain theetcollision 
velocity, to the use of a particular type of switching 


function, for the special case of a symmetric’ system 
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We can illustrate their method by its application to the IsZ, 
and 2p0) states of H.* : They define "atomic" fragments 
b+) given by 


_ o-!/2 
es tie ee ey 


] III-2 
+ ‘ 2p0.. 


which correspond, asymptotically, to 1S atomic orbitals on 
nucleus B and A, respectively. With each such d+ they 


associate a Bates-McCarroll type (single-centre) ETF, 


d= ¢ exp [+ (im/2h) A(R) Ver - (im/ah) vot] III-3 


where ~ (R) is a variable parameter. Finally, they 


reconstruct molecular basis states of g and u symmetry by 


and determine the parameter A (R) using the Euler-Lagrange 
method along the lines proposed by Riley and Green(21). It 
Gateebe, sshown(33), ‘that .(,to firstsorder in velocity) wthis 
procedure is equivalent to the use of the Switching 
functions 
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and 


Rana te =e = 
Fopo (F9R) = )(R) 20, FR 5, (r;R)] TIT-5 
u 


For the 1s@j state, this switching function is in quite 
Surprisingly good agreement with the one we use (see the 
paper by Thorson,Kimura,Choi and Knudson(33), Fig.(7)) (for 
R24,,0<c.8.U.),- but for-the 2p7j: state, it is.very different, 
and leads to some very different coupling matrix elements. 
We will discuss this point further when comparing our 
results for cross sections with those found by Crothers and 


Hughes(48) (Chapter V). 


Independently there is also a recent paper by Taulbjerg 
and Vaaben(28) which attempts to derive a common switching 
function for all states of a molecular electronic system by 
general arguments based on the electronic Hamiltonian. One 
topic to which these authors pay much attention is the 
united-atom limit R — > 0; they especially emphasize the 
idea that as R —> 0 the electronic density must be moving 
withtedtheg centretofecharge; rather than with either 
individual nucleus, or with the geometric centre. AS a 
matter of fact, we find that the switching functions we have 
obtained for HeH?* by the methods used here do seem to have 
this limiting property; however, we find it does not become 


important until very small internuclear distances, and we 
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have found no evidence that such limiting behaviour has a 
Significant effect on the cross sections calculated for the 
systems we study here. 

Much more important to our present work is the fact 
that different switching functions are required for each 
electronic state, and this possibility has not been allowed 
by Taulbjerg and Vaaben. These authors do not seem to. have 
done any extensive cross section calculations testing their 
proposed switching functions (see Section VI.A). 

In this work, I have followed the methods and ideas 
employed by Thorson and his coworkers (Levy & Thorson(8b) 
and Rankin & MThorson(29)) in this laboratory to determine 


Switching functions. I have used the form: 


PR) ee PORTH IRe aR) (h SUR 111-6 
n n nN 


where By (R) and Om (R) are variable parameters depending 
upon the state n and internuclear distance R, and 
We is the "angle" variable of prolate spheroidal 
electron coordinates (-14%7 &+1). This form is closely 
related to that first proposed by Thorson et _al(8) in 
conjunction with the calculation of coupling matrix elements 
for direct impact ionization in H*-H(1s) collisions, and was 
extensively studied by Rankin for ionization couplings in 
both H,* and HeH?* systems. Advantages of this form are 
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(1). Theoretically, as we will show below ( and is shown 
ijepdetaala sith recent work by Thorson, Kimura, Choi, and 
Knudson | (hereafter -’TKCK) (33) )\ this form ‘can . be 
justified by its derivation for H,* using analytical 
methods, and is not therefore a purely arbitrary choice, 
and 

(2). From a practical viewpoint, it is easy to compute 
the nonadiabatic coupling matrix elements for this form, 


compared with other proposed forms discussed earlier. 


In this Chapter, I will describe in detail the studies 
I have made of "optimization" calculations to obtain "best" 
Switching functions of the form (111-6), Similar to the 
Studies made by Rankin(29), but in the present case looking 
also at discrete-state couplings in addition to ionization 
couplings. These calculations give results qualitatively 
Similar to those found by Rankin(29), but slightly changed 
quantitatively, when couplings between discrete states are 
considered. The work described here was done before the 
analytical derivations done by Professor Thorson were 
obtained, and they served as an important confirmation of 
the correctness of the analysis he obtained because they are 
iniexcedlent agreement with hiswanalyticalerestits foreH:? 
In Section B. we discuss the methods and ideas to be used 
and also describe the ideas of the analytical method of TKCK 


briefly. The rest of the chapter is devoted to the detailed 
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optimization procedures used and their results. 


B. GENERAL METHOD 

In this laboratory, two entirely distinct approaches 
Nave sbeens, ideveloped. for the ~determination...of» switching 
functions for molecular bound states. The first approach is 
the so-called coupling "optimization" scheme used previously 
by Thorson and coworkers. The second approach is suggested 
by a. result obtained by Delos and Thorson(12) in their 
discussion of diabatic and adiabatic state descriptions of 
atomic collisions, and is based on an analytical 
decomposition of two-centre Ro yerueee Sac into "one-centre 
components". In my thesis work, I have concentrated mainly 


on the former method. 


1. Coupling "Optimization" Method 

Some time ago, it was found by Thorson and coworkers(8) 
that the corrected non-adiabatic coupling matrix elements 
connecting bound molecular states to the two-centre 
continuum states in H,* are hypersensitive to the choice of 
Switching function used. The problem has practical 
importance because the PSS couplings if not corrected are 
HenVielangeetG:O0eb «hOeks 140eneeUsded MeEYe—ghumerous 2430540 
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very long range (30-50 a.u.). 

Since it is believed that most of these couplings are a 
result of the fictitious "displacement" terms, an argument 
appealing to variational ideas suggests that a "good" choice 
for a switching function should have the result of greatly 
reducing or even minimizing the corrected couplings(Levy and 
Thorson(8b)). After empirical studies, it was found that 
when form (III-6) is used there are quite well-defined 
choices of the parameters B,(R) Chor’ sug? wth =0 by 
Symmetry) for each bound state, which produce very great 
reductions in magnitude of all continuum couplings 
Simultaneously, except those for the first few partial waves 
(see Rankin(29)). The "optimum parameters" obtained are 
independent of the continuum states, energies or quantum 
numbers, and also independent of the type of coupling 
studied (radial Or angular ). However they depend 
sensitively on each discrete state in which the electron is 
initially bound. This suggests they are actually properties 
of the bound states, ce sa timesliS © COGbOras | at hess) wale 
Rankin(29) formalized the optimization procedure more 
Systematically, and also extended iit to asymmetric 
one-electron systems, particularly HeH’* . He was able to 
determine parameters for BAR), Dn (R) for six or seven 
States) “of H,* and HeH?* for OSRS12 a.u., and his coupling 
matrix elements for ionization show very great reductions of 
magnitude relative to the uncorrected PSS values and in some 


cases are also much smaller (at finite internuclear 
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separations) than those obtained using Bates-McCarroll 


ETF’ s: 


2. Analytical Decomposition Method 

Using a classical trajectory formulation, Delos and 
Thorson presented~a unified description of low-velocity 
collisions, in which either one-centre (atomic) or 
two-centre (molecular) states, including the appropriate 
ETF's, could be used as basis states, and considered the 
effects of transformations between such basis descriptions. 
In particular, suppose a basis transformation is defined 


uniquely, 


Ys = my Oe III-7 


which links atomic. states Pp to molecular states Yu : 
Since an atomic state op has a uniquely defined 
Bates-McCarroll ETF whose "switching function" is just fy 
=+1, J=B,A, then, to maintain physical invariance of the 


description under transformation, the switching function Ep 


associated with tu is uniquely defined by the relation 


By thee) eras Xf III-8 


i Loin éluhoees 60 en i! 


sale a 7 ri | 
ra fat ip) te lgrne? eolipotaah’ ‘traivoste ua 
eitetev-vei> 9s webten 1 etek been ee 
teods) e#rtiea- one rechie Hen * 
oi te ee da OReBe ae ae Namba 
ste abel ae: Bea shy eye 
fyoeet elaad douse now tie se . 


btn te ied thet @) \sseni 


4 
* 
re 
a 
the 
be 
> 


x, Vee ae head Aa 
i 


go <etese talyneiém oF r #oy6iea™ places peers 
igat tet 4 fe pps J 6 * Bait ‘ah aese  /? ne 
es seer el aNIEN of asia ite’ “8 as 
a¢3 f -sOnmat seus! cero ei | qisvabem ‘az rr ited i 
mi rattan) geidorive\ ade Pper este ere Tbry, sqtto 
(olieéles ad -yd tiphl ded erry ae Pg dite belgian La 
os 2guee 


68 


Where. wJpiRBsOtyA nespectively, if bp is attached to nucleus 
B or A. This decomposition formula shows how a switching 
function represents, in.a local way, the degree OF 
attachment to each nucleus. Therefore, if a "two-centre 
decomposition” of molecular states into one-centre 
components could be made, it would be possible to define 
Switching functions from the properties of the bound states 
evone.7 FOr Hosa peslnOrsoni» Kimtira ,/iChol,. and Knudson(33) have 
shown that it 1S possible to achieve an analytical 
decomposition of the "angular" wavefunctions (in prolate 
spheroidal coordinates) which leads to switching functions 


of the form 


f = tanh[8 (R) Rn] III-6' 
u u 


The parameters Bu (R) may be calculated by a uniquely defined 
analytical formula. 

For asymmetric systems like HeH’?*, the same ideas may 
be extended in an approximate fashion to obtain the form 
(III-6) as a result. A discussion of the detailed analysis 
of TKCK(33) is outside the scope of my work; however, the 
extremely close agreement found between the parameters 
obtained by the “optimization” calculations I describe 
below, and the analytical values for H.*, confirms that the 


ideas are correct. 
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C. OPTIMIZATION CALCULATIONS 


1. Overview 

The principles and methods used in my calculations are 
basically the same as those employed by Rankin(29). 

We use a Switching function f of the form (III-6), with 
variable parameters lop lucay We (R) to be determined at 
each pRaw BOr.a thansition nus nte(wheresthe final eSstatecn' 
may be either a discrete state or a continuum state), we may 


define a "reduction parameter" ees ; 
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where Pain and Amn are radial or angular nonadiabatic 
coupling matrix elements defined previously 
(Eqs.(11-22b,c)). The basic hypothesis of the optimization 
scheme is that many of the uncorrected PSS couplings i 
are spuriously large due to "displacement" terms; we 
therefore look for a particular choice of the parameters 
Ao wR), Tn (R) which achieves large simultaneous 
Beductions, inipiasw manyee Of a:theadenyvor B62 Possible Whiis 
hypothesis cannot be rigorously derived from formal 
Variational. principles; .it ,aisvedustifaed byydtsjempirical 
Success, 1i1.e., such parameters can be found with 


well-defined values. 
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Rankin's study was confined to continuum final states 


/ 
n'. For a given continuum energy, € 


, ana a particular 
type of coupling (radial: 4 |m|=0; angular: A|m|=+1 or 


4|m|=-1) he defined sums of reduction parameters, 


= -10 
S ee Sa i 


the sum being performed over the continuum partial waves. 
The parameters eye were then varied to minimize S* . 
However, Rankin did not include all partial waves in 
these sums. He found that the first one or two partial waves 
have matrix elements which are relatively insensitive to 
changes in Ban and Ges (see Section C.2 below); since these 
couplings are also large, if they are included in the sum 
their presence tends to mask the effects of parameter 
variations on the couplings to higher partial waves. After 
making these exclusions, Rankin’ found that precisely 
determined values for Priel can be located which very 
greatly reduce for'in for all@rerms in” the? stums**{111=10)woat 
once. In some cases (for higher partial waves) Pas may be 
reduced by several orders of magnitude. The values of the 
Parameters obtained are independent of the continuum state 
energy €’ and of the type of coupling considered. By 
restricting attention to higher partial waves in the sums 
5 sh a 1 Ot Rankin was able in favorable cases to determine 
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somewhat lower precision (about two significant figures). 

After study of couplings to discrete states, I found 
that some of these couplings are also quite sensitive to the 
Switching function parameters. I also found that the 
precisely determined parameters of Rankin do not always give 
an “optimum” reduction of these discrete state couplings; 
nearby, but Slightly different, parameter values are 
indicated. As was found for the lower continuum partial 
waves, some groups of discrete state couplings appear 
insensitive to the switching function parameters. 

As a result of these investigations, summarized in more 
detail below, I have reported "optimum" values for’ the 
parameters Bewnn which differ somewhat from those given 
by Rankin; they also have a larger uncertainty (about 10 % 
in worst cases). (In one case (the 3d Oy State for Hee) 41 
have determined a quite different value from that reported 
by Rankin; he seems to have found a secondary minimum in 
that case). 

From the standpoint of overall method, there is no 
Guestion thaticouplingstto continuum ostateso areri thelemost 
sensitive to the switching function parameters, and provide 
the most definite general indication what values the 
"optimum" parameters should have. (I found that couplings to 
Rydberg states with AL +0 show the same behavior, and 
agree with continuum state couplings for determining the 
region of "optimum" values). However, after establishing a 


range (+8 %) for the "optimum" values from the continuum and 
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Rydberg Steie “couplings, .. Study ).0f.. thes discrete, ustate 
couplings (as well as continuum couplings to lower partial 
waves) suggestS a choice for parameters which can deviate 
Slightly from the precise values found by Rankin, who 
considered only the higher continuum partial wave 
reductions. The new choices still make large reductions in 
the continuum ‘couplings, but also make a_= significant 


reduction of discrete state couplings. 


The symmetric case of H,* ivi =0) is much simpler, 
and the parameters 4, may be located somewhat more 
precisely, than for the asymmetric case (HeH?*) studied. 
Therefore for clarity and convenience I have presented the 


two cases separately, with H,* first. 


2. Bn -Determination in H,* Couplings 
For most (but not all) coupling matrix elements n-> 

n', it is possible to select parameters for the switching 
function f£ which cause the corrected coupling matrix element 
( Print Ann) to vanish exactly. Some of these zeroes occur 
at isolated values of the parameter for the specific 
couplings (they could be related to nodal structure in _ the 
coupling), but others are systematic in the sense that many 
different matrix elements have zeroes at or near the same 
parameter value. This is the situation in the continuum 


couplings, which Rankin used. Of course, we are not looking 


ions 


| ' 
rh! aie Pax 
6h Oa ae sn 


ese ederie2s, wh feo 
‘ote «awe -a ie 


ieuivet tao. tobhAw nese 


fir aia oF 
GN re ; 
so2 20h 
. Wi i sy 
2 ‘wen Yue 
; e ; " bed 7 
‘ Rel) Bere Se 4 


| eli nial 

jthemele x22 ee ant aver (Tle sey sis 90m ant 
ealtisiwa, af? se aresdqs tay Faken og aidiaeeg ataly 
tiurel) xl 376m atid Cage o5 bardeatis edt veuan fotite t 
thase  2eer,es, weed? 36 Peers arerte deinae ef 4, 
shitionae adi wed A iemeitieg, “eNg) .29 aeuligy tan 
$62 .nt. esvawivie tebe) os staal ae Sony 
yaen jel? semee ors, vi ok tamed ara hie aIeAeB. ce 
was oft weer, Joo de: asates Svad astnemeie ingen’ 


Muvlerss ety al mc ited bee eta aE aan vache Saal 


ethos, 200 826 aM (denyoa-20" somaw nbangn okie 


dss 


for such zeroes exactly, but if a large number of them for 
different couplings are clustered in a narrow range of the 
parameter, it is reasonable to conclude that somewhere 
within that range there is an "optimum" choice which will 
greatly reduce the magnitudes of all of the couplings which 
have nearby zeroes, and that such a choice has physical 
meaning for "displacement" corrections in the state n. I 
have found that when discrete as well as continuum couplings 
are considered, such a clustering of zeroes is also. found, 
but that the range is a little bit broader and the "best" 
choice appears to be a little different than in continuum 
cases alone. 

For both discrete and continuum coupling matrix 
elements in H,* , the general pattern found is as follows 
(let 1 be the united-atom orbital angular momentum quantum 
number): 

a. «For couplings y"=l,  the..matrix elements have 
non-zero values fon abl (oho Vauues,) dO.) note haveg a 
minimum for any reasonable (Bis avalue*( 050 5£0-'170) 
and are relatively insensitive to (Bre (see 
Figure (Iti— fa) ) 

b. Couplings with 1'=1+2 have two zeroes vs. orn in the 
reasonable range (Figure(III-1b) 

©. Couplings with J,'=1+4 have three zeroesy, vs- Bn in 
the reasonable range (Figure(III-1c); for AL> 5 
we can expect this behavior to continue. 
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(0) 2po,- 3po, (b) 2pa,- 4fo, 
R=12au 
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0.2 04 


Figure(III-1) A-values vs. B for (a) 2pos - 3pa5 (1'=1) 
coupling (b) 2p0q - 4f£05 (1'=1+2) coupling, (Ce) "pay — 6ho 
(1'=1+3) coupling at R=12 a.u.. The corresponding P-values 


at R=12 a.u. are also shown with dot. 
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with one particular set of zeroes in each group of 
sensitive couplings, can be defined, in which a 
common "optimum" value can be found for both 
discrete and continuum transitions (see procedure 
below) } 

e. These features are independent of principal quantum 
number n' or continuum energy €’° . 

f. They are also independent of the type of coupling 
considered (radial or angular). 

g. In general, the effects are sharper and more 
definite at larger R-=values (28.0 a.u.), and the 
total reduction in couplings is also greater. This 
is to be expected from the theory of "displacement" 


effects. 


Using these facts we have developed the procedure to choose 
optimum A, 's for H.* states as follows: 

1) Perform a systematic calculation minimizing sums 
(III-10) for continuum states only, repeating 
for different continuum energies and different 
types of couplings. Omit the first two or three 
partial waves to obtain a "best" continuum (-n 
(Rankin's procedure). 

2) Define from this an "allowed range" for further 
search, at most + 8 % around Bn . einteiciis 
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reduced, but not absolutely minimized.) 

3) Considering discrete-state couplings for as many 
transitvons with | A Wie2 sas possible (both 
radial and angular), locate a range of "best" 
discrete parameter, A, : 

4) The two ranges will lie within one another or at 
least overlap each other. Choose the optimum 
fm from within the overlapped region. The 
uncertainty in the value corresponds to width of 


the overlap region. 


Values for at the resulting parameters are listed in 
PacrecrL ii fOr 16 States Of Ho fat R=2-4678 -T0-ande™ 12 
a Sica CSCC GSCUSS1.0NW~ObeRESUIES) Gisubelow.)..(A..Similar 
Procedure can be used to locate Be fOr Hen states; alter 


the determination of dhs ee 


3. Determination of Parameters for HeH’* 

In the asymmetric case HeH’** , behavior of couplings 
between molecular states is strongly influenced by their 
asymptotic connection to atomic states of He* or H, except 
perhaps at small R values. The parameter as gives clear 
indication of the limiting atomic character of a state, 


since it is the "crossover point" for which the switching 


function changes sign: for example, if the He nucleus is the 
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Table III.1 H,* Switching Function Parameters 


+ 
H, Switching Function Parameters £g a (u-state) 
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£ (r/R) = tanh[ g Ry) 
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Table III.2 H,* Switching Function Parameters 


R (a.u.) 
State 
2.0 10.0 210) 
0.065 0.065 0.066 0.067 0.068 0.068 
0.051 0.052 O505s 0.054 0.054 0.055 
4£6 
u 
5£6 
u 


ae 0.061 0.061 0.063 0.064 0.065 0.066 
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Table III.3 H.* Switching Function Parameters 


+ 
Ho Switching Function Parameters 8 (g-state) 
£ (rR) = tanh[s nkn ) 
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A nucleus (7 =-1), then a state for which 7), 20 (or even 
>+1) is evidently a "He*-like” state. The 1s@ and 2p7 
States of HeH?* are clear cases of such states. At large 
enough distances R, some HeH?* states also can be found with 
m,. <0, corresponding to "H-atom" states; the 2p@ state is a 
good example. However, as a rule the "H-atom" states become 
really molecular states at smaller R-values and have 
two-centre character; some "He*" states (e.g.3d”M ) have 
some mixing and two-centre character,but it is usually less 
important. 

a. For a completely atomic (one-centre) state, we know 
that Bates-McCarroll ETF's are correct. This means, 
foutstch wattistave, Sonr=anWe(He), Sorcri;e=+10° (nH) 
everywhere, which would correspond to >> +1 Otte Kc 
-1, respectively. For such a case, we should then 
find that couplings to other discrete states should 
be quite insensitive to the parameters of the 
switching function Brn, Wn . Any choice of these 
Parameters which gives f= -1 (or =+1) everywhere 
will give good behavior of the matrix elements for 
most transitions. This is what we find, especially 
at larger -R-values; for the "He*" type states like 
ie Cee} 26 0- Ov 2p ete. UFPor S2beistates (li kene2pge 2) 
this behavior is found only at very large R values. 

b. For the same reason, even at smaller R-values, 


coupling matrix elements from such an "atomic" He* 


State to other tightly bound discrete states are not 
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very sensitive to parameters Bn, Mn of the 
Switching function, because the switching function 
is close to -1 everywhere in the region of overlap 
of the wave functions. As a rule, we find that for 
States where mn 2+0.2, couplings to discrete states 
are not very sensitive for this reason. For. such 
states most information about the optimum parameters 
is) foundytrom ~the acontinuum couplings (and also 
Rydberg state couplings) because their wavefunctions 
can have overlap in regions where f is different 
fromaast. 

On theeronher whand aston thes molecular’ .type.of 
States, both the continuum and discrete state 
couplings are sensitive to the switching function 
parameters, and determination of optimum values can 
be made using both types of data. 

For those couplings which are sensitive to switching 
function parameters, qualitative features similar to 
those described in items (a)-(g) for H,* (C.2 above) 
are found, but most consistently for couplings to 
continuum and Rydberg states. To some extent, the 
location #of.,.“clustering” zeroes vs. Bin can be 
affected by the value taken for Nn anGsen~ltervs 
therefore important to locate the correct Wn value 


before fixing Bn more carefully. 
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The procedure used to find parameters for HeH?* is as 
follows: 

1) Perform systematic calculations minimizing sums 
(III-10) for continuum states only, repeating 
for different ae and coupling types. Omit the 
first two or three partial waves in such sums, 
to obtain "best" continuum parameters 
(Rankin's procedure). 

2) From these calculations an "allowed range" for 
further search may be defined (+8 % about BS 
nen Nsyee ny 4 about Jn ). Within this range, continuum 
couplings are greatly reduced, though not 
absolutely minimized. 

3) The range of acceptable gi, may be reduced to 
+10 % by consideration of couplings to Rydberg 
states. The “range sot Bathe obtained by the 
above steps iS narrow’ enough to identify 
"atomic" and "molecular" type states by their 
parameter values, and for "atomic He" states 
like is& OTra2P JI , the switching functions 
are adequately defined at this stage. 

A) Ol essenOLecllanm iat YPeesStates, couplings to 
discrete states are next examined to determine 
which couplings are sensitive to Switching 
function parameters--within the ranges indicated 
by the continuum coupling study. (Because of the 


complications introduced by atomic character and 
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the asymmetry of the system, the couplings among 
the tightly-bound discrete states do not fall 
into a simple classification scheme like that 
found for H,* ). For as many discrete couplings 
as possible, the locations of clustering zeroes 
vs. B. and mn are determined, and. the final 
"optimum" values for these parameters are 
obtained from the region of overlap between the 


discrete and continuum parameter ranges. 


Values for the resulting parameters are listed in 
Table(III-2), for 12 states of HeH?* at the same R-values as 
for H,*. In general the uncertainties in these values are 


Pangers than cNOSe LObe Hoa. 


4. Behavior of Optimum Parameters (H,* and HeH?* ) 
anaes 

Table(11I1-1) shows # values vs. R for a number 
of H,* states. The values reported for 18 Oe Pe2Dlus 1-3 
2Dlo. cs ane 250, are 5 % to 10 % larger than those 
previously reported by Rankin, and the estimated 
uncertainties are larger. For the 34 0g State, a 
quite different set of A, values is found. This 
value corresponds to a different minimum in the 
continuum couplings than that selected by Rankin, 


and its correctness is confirmed by the behavior of 
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Table II1.4HeH?* Switching Function Parameters 


f- 
Hen” Switching Function Parameters G. ne 


£ (7,R) = tanh[8 R(n-n°)] 


1so 
<r 0.554 0.541 0.539 0.537 0.533 
2s0 
ae few fo foe Tie fe 1.17 
ean og 0.732 0.278 0.076 | -0.145 | -0.226 | -0.332 
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2p 
ae ee 0.970 0.878 0.798 0.828 0.837 
ie 0.232 0.257 0.266 0.274 0.282 0.286 


3do 
0.724 0.513 0.484 0.485 0.506 
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Table II1.5 HeH?* Switching Function Parameters 
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the discrete couplings, which do not have clustering 
zeroes near Rankin's value for Brag weeContirmation 
for the values I have obtained here is obtained by 
their extraordinarily close agreement with the 
parameters Bp. determined by Thorson et al using 
analytical decomposition of the angular wave 
functions (see Figure(III-2)). 

An interesting regular behavior pattern is 
observed in the ratios of Brn for succesSively 


increasing principal quantum number: (n21+2) 


ero 070 fOr 05 III-lla 
Baim §31m =: 0) betta. 05 III-llb 
Bs im Pain = 0-8 +20.02 eae 


irrespective of the value of Rand of (1,m). These 
and certain other regularities in the behavior of 
the parameters om are confirmed and explained by 
the analysis of MTKCK. According to this analysis, 
the quantities bh», = Bb R can depend only on 
(1,m) and the energy parameter c, Caer Riz: 
Hence 
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Figure(I1I1-2) Parameter 8 vs. R for H,°* states. Solid 
curves show analytical results. All data points are values 
from “optimization” calculations. (@,@,4,8 , valves of 
Rankin and Thorson(29) for 180, ,280, .2PT, , 309% ; 
©,0,0,4,% 8S modified by inclusion of discrete-state 
coupling effects). Error bars indicate uncertainties in the 


latter. 


87 


| Te 
bites .eviane * 4M cod # 4) epapeaiel (S~ 
deutay ow aonieq aoaf (aA aan fastiylaes vale 
jy aeatay , &, das tal diaiibbuetes “opting aa T 
ee ee | 10t (9B nowaeitt bea sbtoet a 
siepeneretaai®, 36 nokia aa “ Bah Aggies) A 


i) ve eeeerberecay etapigeb dyed’ evr tarap tio gqwhtqves: 


es 


ae i 


=i") a 


88 


For the lower lying states (e.g. 250 / Bec, ,etc) the 
ratios deviate from the rule (n/n+1) because the 
electronic energies deviate from the Bohr formula 
~(1/n? ) due to bonding arrears however, the first 
equality smelileaii)), Stodtaholrds. Other regularities 
in the H,* parameters may also be deduced from the 


analytical theory. 


HeH? * 


Table(III-2) gives the parameters 6 my. vs. 
R for HeH?* states. 
In his study, Rankin proposed a scaling law 


relating GC tor ‘Hela to 6 Values tor, Hi-.,, 
we Oa 


and he argued on the basis of continuum matrix 
element behavior that this relation was confirmed. 
Comparison of the results in Tables (III-1, and 
III-2) shows that this rule is not very accurate 
when the discrete state couplings and a more 
complete study of continuum couplings is made. 
Further studies on this subject should be 
pursued to understand the characteristic features of 
the detailed switching function in the different 


states in the different systems. 
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Mésanalysisepresentea jin TKEK(33) for H,* may ‘be 
extended only in an approximate way to asymmetric systems 
mike HeH+*)-. For ~the more |"molecular" ‘types of states, 
estimates of the parameters Boe based on decomposition 
of the angular wave functions have been made. These agree 
fairly well with the values I have obtained here (see 
Figure(III-3,4)). 

If empirical regularities of the type found for H,* 
could be found here also, they could be useful to predict 

C values for upper states without performing optimization 
calculations. Some regularities of this kind are observed, 


but only between states which belong to the same nucleus as 


R—voo and also large R(28 a.u.) only. 


D. COUPLING MATRIX ELEMENTS 

Figures(III-5) to (III-15) show nonadiabatic couplings 
calculated using switching functions determined by the above 
procedures-for the H.° “and “HeH® systems. For practical 
applications in the collision studies done here, I have used 
linear approximations to describe the R-dependence of the 
parameters C. ; ns for each state. Such an approximation 
appears quite drastic in some cases (especially for the 
behaviour of aieat small R) but in fact it has only a small 
effecteHon: the*®coupling matrix elements and no effect on 


computed cross sections. 
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Figure(III-3) 
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More than 240 corrected nonadiabatic coupling matrix 
elements have been computed for the H,* system, and more 
than 130 matrix elements for HeH?* (this does not include 
the couplings to continuum states which were studied in the 
optimizing process). Figures(Il1-5) to (ItI-15) give only a 
representative sample of these matrix elements, but include 
especially the couplings which are most important in the 
collision processes studied in this work. I. shall not 
discuss individual couplings and excitation processes in 
detail here (see Chapter V). However, a rough classification 
of the different types of couplings shown in the figures is 
helpful to understand where ETF correction effects are most 
likely to be important: 

1. Classification of Couplings 

The efficiency of coupling between two states in 
quantum mechanics is affected not only by the size and 
range of the coupling matrix element which Baie ous them 
but also by the energy spacing between’ them. In 
particular, EOS low and even intermediate energy 
collisions the most efficient excitation processes 
always involve a resonance or near-resonance of two or 
more levels, that is, a "mediating degeneracy" which 
makes efficient coupling possible; such a degeneracy may 
be associated with the asymptotic limit R —»ya , with 
the united-atom limit R -—» 0, or with an isolated 
crossing or avoided crossing of levels at some finite 


R-value. For the one-electron systems’ studied here, 
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crossings and avoided crossings do occur, but they play 
only a minor part in excitations. Couplings ‘associated 
with orbital degeneracy of united atom levels are the 
most important excitation mechanisms in these systems 
(ror example, the .Strong»srotatronal 2pd; -2pit, 

coupling), but there are also some important couplings 


at very long range connected with degeneracy as R— > oo, 


Below collision energies of st KeV/amu, 
degeneracy-mediated couplings are really the only 
effective ones for excitation but in the intermediate 
range 1-10 KeV/amu this simple picture (based on _ the 
"Massey adidbatic “Criterion (2) ) "sis = complicated by 
increased efficiency of excitation processes which are 
not so clearly connected with a mediating degeneracy. 
This) is what makes collision processes “in this’ energy 
region so much more interesting, and harder to predict 
without actual calculations. As the energy increases 
Still more, above 10-15 KeV/amu, so many "non-resonant" 
excitation processes play a part that the idea of 
mediating degeneracy loses its significance. However for 
the intermediate energy region it is still very 
important and so energy gaps as well as coupling matrix 
elements must be considered to understand excitation 
processes. 

Looking at coupling matrix elements themselves, it 


is useful to group them into two main categories, those 


a ee 


a, i 
ik ee, Y i 


A eee i 


yeig: ysnt Jue S08, shin 
beet oomte: agri teed | ov 7 — 
sit oe elsvel mode bis ting! aa 
muta Je ‘* wets om ~hetnneiell | 7 

jwgs wat fucei reyes or) " - a fe 
eprrlqweg 2 est Oye S ~ ROA dis = Li 


pee T Be YoRCenepeb Lain, 


ee 
aT ae Caeiats 


une Vee , 2 ty es iar, s faa 
eine gra ef pert, ots’ scot: th 
a7 ga penta ou ony we one. net seatoue per ome 
afd no srad) Oth +E bathe. | fe 
vo Batenileeas vey Br: , 
ayhedw zor 2H ION Ad 7 | 
i teTsd0pah Powe mebaae, peer by 
‘7 ate seas a snag 


aon canes Yosene inthe ea | 
"tonceeet-mon” yew oe paar ovods waded tt ae 
to ambi Gay | Peed, Smeg yatq Hewessorg: ngtsgdioxd 
wo3 9OvSwOR sama ceute eat oaes ‘gost sanpeb enizetbem 
qi 806 Ered) St oetgat Geiens etelisarsynt ws 
kitesm onliquos 2é Lise aa ap Nersae ov bas. ‘head edeqel . A 
Koizaiion bneyerstw Biack —. od suum esnemsts. : A 

-** seenasgeg 
$k ,asvhonnsild e ‘etek abide pnt igen te purieal | 
meoes .eeitepagess em: owt athe med? quetp o2 tuigew €f Ae i 


94 


which involve short-range coupling and those which 
involve couplings with very long range. [In a few cases - 
a coupling may show both characteristics-for example, 
ces sou. 50%, coupling in Ho serig.( 111-4), but this is 
not usual]. 

For the systems studied here the short range 
couplings mainly occur at distances RS12 a.u.. They are 
associated with "quaSimolecule" formation or even with 
united-atom behaviour. It is useful to divide them into 
Subsets; first, the rotational couplings connected with 
Orbital degeneracy of the united atom (for example, 2p 

Ur a2 pi my sad (-sOlme topes 2p. syetc ies and ssecondly; 
the many short-range radial couplings. The primary 
excitation steps in the H.* and HeH’?* collision systems 
all involve these short-range couplings. 

Long-range couplings may involve Significant 
effects out to distances as large as 100 a.u. or more in 
some cases. They can be understood as perturbations of 
atomic states by the collision partner, which is an ion 
in these cases and causes Coulomb field perturbations. 
Here also we can distinguish two classes of long-range 
couplings relevant to the present case; ELES ts 
asymptotic couplings within a degenerate atomic level 
(both Coriolis and radial couplings), and secondly, some 
long-range radial couplings between non-degenerate 
levels caused by polarization in the ion field (for 


example, %wW-W couplings in H.,*). The long-range 


Gain awed 
go909,. 99) a 


eine, 199-24 


To. HOA RST +d 
oon: these vs SG NVR a | titi ky 4 
iw beeenee> ahad ined hae t 
ws , otmarmite, 1O8)), wad io: pretn, 
\VIPAOSCSS ONE) 4 ( ota ny igi ‘ i Nig. an 
cinmlag ith lame tga Bene spn ay! 
agetaya Te tie pL oro a ‘ine Ai + at ‘2 
00 Laue 2 


SOF PAA Bir IAQ ch onatystni | 
hod nat an 
IG Steet tah Sy hie! ome hy 
epee rel iy pee Rue? Awd wens Meo ahd er coset 
~aent! br el Set coal tiaveiar ape} Lawes, 
fove!l ainoze asyaerensoon 4 | | 
Sgoz vitnous? han 20nkk Snitar hoe ei falve® Aiod) 
nepeh-ndn esha ad age a Terbes.. » Sener yioT 
piebt. not .sneyve icy cs dll im dame ) 7 
es 9AT .(" Se WE SpReiniies- ar. agers 


7 wee) io . 
7, a ' ote e Ee 


35 


couplings play an important secondary part in producing 
further excitations following the primary (short-range) 
excitation steps. 

The distinction between long-range and short-range 
couplings is important when considering ETF corrections 
to coupling matrix elements. As pointed out in the 
discussion on PSS theory, the most obvious defects of 
PSS theory are connected with its failure to describe 
asymptotic (R — BD ) behaviour correctly; many 
couplings which actually have short-range character 
appear to have long-range behaviour due to the spurious 
displacement terms. However, these defects are also the 
easiest ones to remove, uSing almost any type of ETF 
description, because the states are almost atomic states 
(or linear combinations of them). At large R-values, 
therefore (say, R212 a.u.), there is no_ significant 
difference between corrected couplings computed using 
our switching functions and those computed with (say) 
appropriately constructed Bates-McCarroll-type ETF’s, or 
with some other switching function. Therefore although 
the effects of long-range couplings are interesting by 
themselves, they do not provide any test of switching 
functions in comparison to other ETF descriptions. 

On the other hand, the detailed form of the 
short-range couplings, which depend on molecular state 
behaviour, can be very sensitive to the choice used for 


ETF description. 
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Therefore we can expect the most’ significant 
effects resulting from differences in ETF descriptions 
to be seen in the intermediate energy region 1-10 Kev, 
where short-range non-resonant radial Or 
degeneracy-mediated angular couplings can play an 
important part in excitation. As I shall show in Chapter 
Vand Wi, this turns outlto be the case for both ‘of ‘the 
collision systems I have studied. 

In comparing the behaviour of radial and angular 
coupling matrix elements shown in the Figures, it should 
be remembered that the radial coupling matrix elements 
must be multiplied by the radial velocity R, but angular 
coupling matrix elements by 5 =Ves7R*, torobtain the 
actual couplings; this moderates the apparently long 
range of many angular coupling matrix elements shown in 
the Figures. 

Examples of Coupling Matrix Elements 
(a) H,* couplings 

Figures(111-5) to @i#I-9)e) show Wall the coupling 
matrix elements used in the Ht-H(1s) collision 
calculations, for the 5 u-states and 5 g-states basis 
sets listed in Table(V-1). Examples of most of the types 


of couplings discussed above appear in these figures. 


ne 


¥oa O¢ +? colget youene wind Pee 
(ahi. {Ph ACOECT-O | e 

o tole aon techlowe Jahre 
hd eit nd wore Liane 5 a x 
wis hp) wiped #03 4089 a> Pr ot 8 
be ing +“? 

raLuete bas LA ‘Da _ wa tvintiod 3 | 
fa Sao kohe a vance wh tbs ett agi 


Birks i a As Oh Be tape 7 he 
aan a 


re tunik twee: yh Cai secain feline 
ait Aiadee se at & a iby fad 


ed 


2 


pitinues . eee, 084 5 oa eons away 
oope pt Lao (it yuer8 ad ee ‘eengmete Pre 
piasd eetese-g 2 “ae eds 202 «etabteluaies, 
2euiis ale » J2op 20 Leight we bereil ages 

‘canes sgane of sg mnie TTT: 


ome 2 ay. 


oF” 


+ 
H urs tate 


ay SEE 
2P0., coupling 


Radial Coupling Angular Coupling 


Figures(III-5 through 9) ETF-corrected nonadiabatic coupling 
matrix elements for H,* system —— , n—Pn' coupling 
and @—@ , n'— »n coupling (ée.0.7, n=2pQy in 
Fig.(III-5)(shown on the top), and n' are shown in each 


figure.) 
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Figure II1.6 ETF-corrected nonadiabatic coupling matrix for 
H + 

; n corresponds to the state described on the left 
Side in the figure, and n' corresponds to the. state 


described on the right side in the figure. 
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Figure III.7 ETF-corrected nonadiabatic coupling matrix 
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Figure I11.8 ETF-corrected nonadiabatic coupling matrix for 
Hs 
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Figure II11.9 ETF-corrected nonadiabatic coupling matrix for 
Bow 
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(b) HeH?* couplings 
Bigumes(til—-10) to (111-15); \show “coupling matrix 


elemenes (for Heh”, 1n ‘particular forvcouplings linking 
2pq- sUinitialsstate for He**-H(Gis) collisions) and 1so- 

(initial, state for H*-He*(1s)) collisions) with other 
States? and yior couplings involving “other important 


States in the excitation processes which occur. 


(c) Asymptotic couplings in degenerate manifolds 
The atomic H(2s),H(2p) levels and He’ (2s) ,He* (2p) 


levels are degenerate asymptotically. The molecular 
States involved are (2pm ,3pG% ,4f0u ) and (3d, ,3d05 

12805 JernoHs anoseece i. oCcG: 6, 260ge" ) “in HeH?*. 4 In 
each set of three states the two @® states are 
correlated asymptotically to sp hybrid (Stark field) 
orbitals..-They-—“are™ both )»coupled..to the W level) by 
Coriolis coupling [within the atomic level, the matrix 
element of the atomic orbital angular momentum ‘Ly 

between 2p, and 2p, is a constant], and they are coupled 
to each other by a radial coupling which decreases as 
R-? (the effective Coriolis coupling also decreases as 
Ras when the angular velocity is included). The 
Solittings within each set of states ~also  sdecrease sro 
zero as RR“? asymptotically, <endP=chereresult rs” =a 
long-range coupling problem which we found most 


convenient to solve separately by a method described in 


Chapter IV.(Section D). 
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Figures(II1I-10 through +5) ETF-corrected nonadiabatic 
coupling matrix elements for HeH’* SyStem c= , ——> 
n' coupling and @=—=——@ , n'—»n coupling (e.g., n=2po 

in Fig.(III-10) (shown on the top), and n' are shown in each 


figure.) 
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Figure III.11 ETF-corrected nonadiabatic coupling matrix for 
He? * 
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Figure III.12 ETF-corrected nonadiabatic coupling matrix for 
Her*’ 
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Figure III.13 ETF-corrected nonadiabatic coupling matrix for 
He? 
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Figure II1.14 ETF-corrected nonadiabatic coupling matrix for 
He=< 
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Fi I1I1.15 ETF-corrected nonadiabatic coupling matrix 
‘Figure : 
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Figure(III-16) shows typical examples of these long 
range Coriolis and-'radial couplings | along with some 
examples of short range radial and angular couplings in 
HeH?* system. (Note that the couplings shown in the 
Figure are not multiplied by the radial and angular 


velocities, respectively.) 


(d)ebongarangetradial (iN -"\ ) couplings in H,? 
As ishown/}in Figure(11I-9), the 2piy's3pT, radial 


coupling has a long "tail" extending well beyond R220 
a.u.. This effect is due to polarization of the diffuse 
atomic orbitals by the ion field. The phenomenon is 
quite general and Figure(III-9) shows a whole series of 
these YW - W radial couplings for W states of the H,* 
System. Asymptotically these tails decrease as R™?, a 
result which can be predicted by considering the field 
of the distant proton as the R-dependent potential (seen 
from the atom to which an electron is bound) and _ using 
the Hellmann-Feynman theorem to compute the corrected 
nonadiabatic coupling. The same couplings can be_ seen 
among Ng states. 

We have found that’ these couplings play a 
Sugniticant pale mma -higlS JencOLusvOns as.Siicestney. 
provide the mechanism for sequential excitation of 
higher Rydberg levels by a "ladder climbing" process. A 
substantial portion of the flux involved may eventually 


appear as ionization. 
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Figure(III-16) long-range coupling matrix elements for both 


Coriolis and radial couplings along with short range radial 


and angular couplings in HeH** system. 
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IV. METHODS OF COMPUTATION 

In this section, we define some mathematical details of the 
wave functions and coupling matrix elements and discuss 
briefly the computational methods used in this study. The 
discussion can be divided into four topics: 

A. Calculation of electronic wavefunctions. 

B. Calculation of coupling matrix elements. 

C. Numerical integration of coupled equations. 

D 


. Long range coupling problem 


A. CALCULATION OF ELECTRONIC WAVEFUNCTIONS 

The one-electron two-centre problem has _ well known 
solutions. The electronic Schroedinger equation for fixed 
nuclei is 


2 >? 2 2 = : : ars, 
[-h/2u)V) - (Z,e /x, + Z,e°/r,) ce (R)1Y, (Bs) 


IV—1 


where E,.®) is the electronic binding energy. Introduction 


of prolate spheroidal coordinates 


E = (xr + r,)/R F ce auace 
n= (xr, - r,)/R j RG oe 2 IV-2 
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leads to separation of the eigenfunctions in the form 


> 
v, eR) = FL (ENS (n)® (6) TVe3 


The factors obey the eigenvalue equations: 


a. Radial equation 


22 2 2D 
a/aé (E -1)aF, /ae) + lece a eaere A om /(E “LIF, = 0 


IV-4a 


b. Angular equation 
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IV-4b 


c. Azimuthal equation 


2 2 2 
ado? /dd +m = 0 IV-4c 
u you 


In these equations, q=R(Z, + 2) , paw hE, Ee aa and 
c = ~e R’/2 rae & 3 Ep PS inra.U.. A, is the separation constant 


associated with the angular equation. 
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We represent the eigenfunctions in the forms 


2 = oo ‘ 
A? ae Oy (ele ian c,(E-1/6+1)" IV-5a 


ANE m 
= IV-5b 
2, So eee au 


® (9) = ue Gout TVeSc 


m is the azimuthal quantum number (=0,+1,+2,etc). In the 
united atom limit, R-—> 0, Anre il (+1), and the states 
may be labelled at allR by the limiting quantum numbers 


(nlm). 


To compute eigenfunctions and eigenvalues, the main 
task is the search procedure and calculation of separation 
constant Bene energy parameter C in a self-consistent way. 
The problem is a well-known one and has been adequately 
described elsewhere(35-37). We found it most convenient to 
compute eigenvalues and eigenfunctions by following a_ state 
from its known limiting properties at R=0 and continuing 
outward with a sufficiently small stepsize in R (70.25 a.u.) 
to guarantee convergence on the same state. Eigenvalue 
computation methods are based on that described by Bates and 
Carson(35) and some parts of the scheme developed by 


Rankin(29) are also used. 
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Precision errors in the computed eigenvalues and 


eigenfunctions are less than 10°* and 10°’ respectively. 


B. CALCULATION OF COUPLING MATRIX ELEMENTS 

In Eqs.(II-20) there are nonadiabatic coupling terms of 
first- and second-order in the collision velocity. The 
first-order terms are more important and we discuss their 
computation in detail. The second-order terms are small 
enough to be neglected over the whole energy range I studied 
for the HeH** collision system, and are non-negligible only 


at energies E25 KeV for the H,* collision system. 


1. Computation of First-Order Coupling Matrix Elements 


a. PSS Radial Coupling 


p= -ifi<xm|3/3R|n> Iy-6 
n 


A number of computational techniques may be used: 
(i)Direct Differentiation 

ExplMeit’™@application of the formula (ivV-6) 
requires the calculation (numerically) of 
derivatives of parameters and expansion coefficients 
in Eqs.(I1V-6) with respect to R and this is an 


inefficient procedure. 
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(ii)Usual Hellmann-Feynman Theorem 


Alternatively, we can use Eq.(IV-6) 


-ifi<m| (8/8R)2|n> = -ifi<m| (av/aR)>|n>/(e - ©.) Iv-7 


Unfortunately, the evaluation of the numerator in 
the RHS of Eq(IV-7) involves integrals with 
logarithmic singularities which are also tedious and 
relatively inefficient to compute. 
(iii)Modified Hellmann-Feynman Relation 

A third alternative involves the best of both 
of the above methods and requires neither numerical 
differentiation nor unpleasant integrals. Noting 
that 
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Since 


a 


G Zz os 2 
a a a = ~2/RH . 4 - (e /a,)V/R IvV-10 


it follows that the evaluation of all integrals in 
(IvV-9) is completely straightforward. To check 
numerical accuracy we have compared results obtained 
by all three of the above techniques, and agreement 
to five or six significant figures was obtained. In 
practice we have used the third method as it is 


about three times faster and much eaSier to program. 


PSS Angular Coupling 
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a is the electronic orbital angular momentum 
component on the y-axis ( i to R), measured from 
the geometric centre. In prolate spheroidal 


coordinates, it is given by 
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The resulting integrals may be done 


Straightforwardly. 


c. ETF Correction Matrix Elements Kee ne 


R 


oe = (im/2h) (e - e,)<m|z£_ (r;R) [n> JAILS 
i oi ae = 
rae okt eM ete €,)<m| xf (r7R) [n> Iv-14 
Here 
Z 
z= R/2En x = R/2[(E--1) (1-n*)]*/2cos¢ 


and £,(%;R) (defined in Eq(II1-6)) depends only ‘upon 


ae 


It is most convenient ih this case to do the 
integrals over “i by Gauss-Legendre quadrature, 
rather than employ analytical methods. Again 


evaluation presents no difficulties. 


2. Computation of Second-Order Coupling Matrix Elements 
The matrix elements which must be evaluated for _ the 


second-order couplings are catalogued explicitly in the 
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following equations: 


(im/f) {<k| [V- (5-5, [Vs (24. + (im/i) Pca | n> W-i5o 


). <k| [V+ (S.-8,)1|5><3 [¥> (-an¥, + (im/n) fH. ,8.]) [n> 
j j°°k R el’ °n WV -1Sb 


2 Se a > 2 
<k| (m/8) [(£-1)V + (Ver) (V efor) 


>> >> >> 
+ 2(Ver) (EVV £, + aveV £ )1|n> Rise 


To evaluate these very messy expressions requires some 
unavoidable algebra. In the case -of* terms (111-1542) 
mivolving, Ye ) we cannot make use. of the 
Hellmann-Feynman trick and it 1S necessary to compute 
Gerivatives of the wavefunction coefficients numerically. 
However, all the integrals which appear are completely 


straightforward to evaluate numerically. 
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The most important second-order terms are those of 
Eqs.(IV-15a). The complicated radial-angular cross terms 
(IV-15) are quite unimportant, since in nearly all couplings 
either the radial or the angular coupling is the dominant 


part. 


We found it most efficient to evaluate both integrals 
over § and over “| by numerical quadrature for all matrix 
elements, Gauss-Laguerre for 4 _ , and Gauss-Legendre for 

uy . We tested 10,28,64, and 98 point Gauss-Laguerre and 
10,32,80 point Gauss-Legendre quadratures and found that 28 
point Gauss-Laguerre and 32-point Gauss-Legendre quadrature 
give a precision in every case of at least five significant 


figures. 


C. NUMERICAL INTEGRATION OF COUPLED EQUATIONS 

The coupled equations (I11-19),(II-23) must be solved to 
determine the final state scattering amplitudes am(+;F ), 
given the initial conditions for each collision problem. I 
integrated the coupled equations numerically using time t 
(a.u.) aS an independent variable. For most calculations, 
straight-line constant velocity trajectories were used; for 
cases where Coulombic trajectories were employed, the 
configuration variables (R, @ ) must be related to the 


time t by the Coulomb trajectory formula. 
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I used the method of Bulirsch and Stoer(39) which is 
based © on athe principle of rational “extrapolation, for 
numerical integration of the coupled equations. This method 
seems to have found wide acceptance by workers in atomic 
eollision theory. It is faster and more accurate than 
multistep predictor-corrector methods such as the 
Adams-Moulton-Bashforth method and requires no_ special 
Start-up codes as they do. In my calculations, the relative 
truncation error (controlling the step-size) was 


automatically maintained between 1 X 10°‘and 1 X 10°. 


Gonservatiion tiof Siprobabils tyiwietm. &iChapter ect lb.o)y was 
checked and maintained between 1 X 10°* and 1 xX 10°75. I 
found that in much of the integration , Surprisingly large 
step-sizes could be used (3-12 a.u. in time t, depending on 


the collision energy). 


To compute cross sections for individual transitions, 
it is necessary to combine the coherent amplitudes (11-39) 
to obtain atomic state amplitudes according to formulas of 
thedy typera,(1I=4a),, “andeithenpyntortgntegrabe: tee? resulting 
probabilities for individual states over impact parameter, 
Eq.(II1-40). The integral was performed by the trapezoidal 
methods and multipoint Simpson's: rules <oPreciusron}.ofs ‘this 
last integration depends on the number of impact parameter 
grid points used. Depending upon the detailed cross section 


being evaluated, I always took enough grid points to ensure 
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D. ANALYTICAL FORMULATION OF COUPLED EQUATIONS FOR 
ASYMPTOTIC LONG RANGE COUPLINGS 

In these ion-atom collision systems, the molecular 
States which belong asymptotically to the same (degenerate) 
atomic level are nonadiabatically coupled over a very long 
range, due to effects of the Coulomb field of the ion on the 
atom levels. Both the splittings between the coupled states, 
and the couplings themselves, decay only as R~?, and the 
effects of these couplings must be taken into account if the 
tranSition amplitudes for individual molecular and atomic 
states are to be'‘calculated correctly. For the H,* system we 
required to solve this problem for the H(n=2) manifold 
C2s"200 -2p.), “and similarly “fer ithe He *(i=2) imani fold **rn 
HeH?*. The molecular states involved are (3d&0 ,2sm ,2pTr ) 
in HeH*™ ®Wand “in HP®, “(3panreretos: P2pie*) “for “u-symmetry “and 
(3d0q 1280, , 30M, )~ for g-symmetry. Beyond some "cutoff" 
Gistance R,, couplings to stateS outside one of these 
manifolds are negligible, so (in the cases treated here) we 
have to solve separately a three-state coupling problem from 
Re to Re Byoo. The effects of this coupling may be 
represented in terms of a unitary propagator U(@ :R,) which 


converts molecular state amplitudes at R, to those at R-y 
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(or vice-versa). In our present problems, the incident 
channel states do not belong to these degenerate manifolds, 
hence we really only have to consider these couplings in the 
Outgoing part of the collision trajectory. Because of the 
long range of the coupling, neither R nor the time t is a 
Suitable progress variable for the integration of coupled 
equations. 

To illustrate the methods used we will consider the 
HeH** case in detail; exactly the same principles may be 
extended to the cases occuring in H,*, with at most minor 
modifications. 

The 3d*% and 2se~ states of HeH?* correspond 
asymptotically to the bonding and antibonding sp O -hybrid 
atomic orbitals of He’(n=2). These are coupled to the 2pT 
State by long-range Coriolis coupling, and to each other by 
an equally long-range radial coupling. For R2R,, we _ found 
that both the couplings of these states and the splittings 
between them may be accurately described by analytical 
expressions obtained from atomic perturbation theory. 

To obtain these expressions, we use the asymptotic wave 


functions for the three states, 
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B37 2 e, -p_/2 ns 

Peete 1/22 /4V2n) [op ,cos6,+ (2 p,) le A IV-16b 
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and represent the molecular electronic Hamiltonian as 


A 
la = - = 
Lowel ier on As NEG 


where "A" is the He?* nucleus and B the proton. If the 
perturbing potential is expressed as a multipole expansion 


(assuming ry << R), then we have 


2 0 3 
-Z7 EE i -_ _ = eee -_ 
ay *s i/R r,cos®,/R r,P,, (cos8,)/R IV-18 


1. Asymptotic energies for the three states are found by 


evaluating the diagonal matrix elements of he for the 


three states: 


(1) 2 3 
; = «0,5 = 27R’-/ G72R! = 3/2R +" eter 
«340 / / / 
33 
-) = -0,5 - 1/R + 3/2R° - 3/2R + eee IV-19 
2s0 
3. 
2pT 


2. Asymptotic couplings. The Coriolis couplings between the 


2pm state and 2s~ , 3d are evaluated using the 


operator 
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the radial couplings between 3df and 2s@ are evaluated 


using the Hellman-Feynman theorem, 


3 - ) 
<2s0| (8/8R)> |3do> = <280| [9 (-25/rq)/9R1Z |3407/ (34g ea 
A 


Iv-21 
and expanding the perturbation via (IV-18) (in both 
Cases;7m tne reference origin, for: the electron is thus 
tanenayco be the He” nucleus). The cresulting couplings 


are 


° 2 3 4 * 
= -j - = + QR = 
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V = ipV JV2R" | 1+ 1/2R + CPR oe 8) are IV-22 
2p, 3do0 aa 3dc,2pz 
V = -ipV /V2R7 | i= 1/2R + C/R? ] = v 
Dp, eso 2s0,2pt 


where Pf is impact parameter, V, the asymptotic collision 


velocity and the radial velocity R is 
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for straight-line trajectories. For both HeH?* and H,”’, 


keeping terms up through order R~* in Eqs.(IV-19) and 
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(IV-22) gives splittings and couplings accurate to three 
Signiticance figures for RS tSS(Giwatuin andme"cutof£" 
distances R, of this size were used. 


Coupled equations to be solved have the form 


ia, = Lact Ven ene en! IV-24a 
where 
ee ee oe IV-24b 
Phesspiactings Oe = § ae ale 
“2s0,3do 3/R° 
“2so,2pT 3/2R° v 3/R° Iv-25 
ED ado = 3/2R° + 3/R° 


and Vak are given by Eqs.(IV-22). 
Progress variable Considering explicitly the outgoing 
trajectory case, we set u= f /R (taking positive sign in 
(IV-23)), and define a progress variable 
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oh _ ol 
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so 6 ranges from 0 to "/2 as R goes fromato § . The 


coupled equations are now expressed in terms of § as 


independent variable, 


i (da, /ds) = bu Vien en ere [7144 IV-24c 
with 
980,340 mea 
A -a[6 + 4/ ey 
2s0,2pr “ rVe27 
Z 2 
- 6 - 4 
Apt, 340 a [ JOU Sin 107 2) 
where a= 3/20V,_ , while 
V = - 1/2 cosé[{1 + 2/psiné - Weceins) 
3do,2so0 
V = 1/v2{1 + 1/2psin6éd + Cf ence) IV-28 
2pt , 3do 
V =- 1/v2[{1 = 1/2psiné6 + C/ ant 1 
2p7,2S0 : c y 
For the collisions of interest, Ro. >» f and § thus 


corresponds to values «1/2. The coupled equations are 
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easily and efficiently solved by the numerical method of 
Bulirsh and Stoer described in SectionIII.C., and the 
unitary propagator U(0; 8,) is constructed as the required 
connection between outgoing amplitudes at R, and those as R 


—7 OO. 


roy 


te 86>bOddew 


os Po 6 


beriupes add 96 be2 ouudeney! at tevin se 


a 2G owen? DS 4S Fe esbu dit sem, 


V. H*-H(1s) COLLISIONS 


In this chapter we present and discuss our calculated 


results for direct excitation, 


H*+ H(1s) — H*+ H(nl) 


and charge exchange, 


H*+cH(1s) —=—e)H(nl)y+.H? 


in proton-hydrogen atom collisions at projectile energies 
je7eckeVese eSwitching aiunchronsesobtained) by) the methods 
described in Chapter III were used to construct ETF's and 
compute ETF,corrections) (n.b: for H2%, switching, functions 
found by the "optimization" scheme are in quantitative 
agreement with the analytical results of MTKCK(33)). Basis 
sets with up to 10 molecular states have been uSed in 
systematic cross section calculations, and good convergence 
of cross section values is found for excitation to all 
atomic states with ns<2. 

Our results can be compared with those of other 
molecular calculations, especially the recent study by 
Crothers and Hughes (48), and with a variety of experimental 
measurements of detailed atomic state excitation cross 


sections, polarization of Lyman- QQ radiation excited by 
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collision, and total and differential charge exchange cross 
sections. Special emphasis is placed on the excitation cross 
sections for atomic H(2s) and H(2p) states, for which there 
is at present a very chaotic situation both experimentally 
and theoretically. Later in the chapter, we present results 
of extended-basis-set calculations which bear on excitation 
to higher levels and a new mechanism for ionization in 
We-H(ie)acéllisions, 

Early molecular state studies of H*-H(1s) collisions 
were mostly concerned with resonant charge transfer 
oscillations or the efficient Ly-@ excitation process (see 
Knudson and Thorson), and employed bases of two or at most 
three states (1s 0g -200G4,2pwlu Yeebatenrecalculationsenusing 
larger bases were done by Rosenthal(45), Chidichimo-Frank 
and Piacentini(46), Schinke and Kruger(47) and most recently 
and completely, by Crothers and Hughes(1979) (work published 
after the present study was undertaken). With the exception 
of the work by Crothers and Hughes, these later calculations 
all employed uncorrected PSS theory and completely neglected 
ETF corrections. The calculations by Chidichimo-Frank and 
Piacentini and by Schinke and Kruger also entirely neglected 
radial couplings. 

When compared with the PSS calculations, our results, 
in agreement with the findings of Crothers and Hughes, show 
that both ETF corrections and radial couplings have 
significant effects and must be taken into account to obtain 


a correct description of excitations in this system. 
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Comparison between our results and those reported by 
Crothers and Hughes is also interesting, because their ETF 
description and resulting coupling matrix elements are quite 
different at some points from ours. In a few cases there are 
substantial disagreements between their cross sections for 
individual atomic state excitations and those we have 
obtained. Discrepancies in the total H(2p) excitation cross 
section and in the 2p./2p, sublevel distribution are in fact 
traceable to the differences in ETF's used; however, an even 
more marked disagreement on the cross sections for charge 
exchange and direct H(2s) excitation does not seem to arise 
from this, and itsS origin is so far not explained. 

In all cases there iS good agreement between our 
results and recent experimental measurements, where these 
are available. In particular, in those cases where our 
results ‘disagree ‘with “those "ore *Crothers*and* Hughes, our 
values are in much better agreement with the most recent 
experiments. 

Finally, we have also carried out selected studies with 
larger basis sets (up to 16 ungerade states) to examine the 
behavior and convergence of cross sections for excitation to 
atomic levels with n23. These studies strongly suggest that 
a "ladder-climbing" sequence of excitations via these upper 
levels may be the dominant mechanism for ionization in 


H*-H(1s) collisions at these energies. 
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A. COUPLINGS AND EXCITATION PATHS 
1. Basis States 

We used sets of 4,8 and 10 molecular states for 
Systematic calculations, as listed in Table(v-1). These 
contain equal numbers of gerade and ungerade basis_ states, 
but as we shall see most of the excitation occurs in the 
u-manifold. Table(V-2) lists the basis sets of up to 16 
u-states which we used for selected studies of convergence 
and excitation to upper levels. States were chosen for 
inclusion in these sets if preliminary study showed that 
they are Significantly linked to relevant lower levels; the 
high preponderance of 7 states, for example, is explained on 
this basis, because strong 7] - 71 couplings have an important 
role in higher excitations. 

Figures(V-1) and (V-2) show electronic binding energy 
Gurves vs (R forpthemW@-state baovorot Tabie(vV-") and for 
the 16 u-state basis of Table(V-2), respectively. 

2. Couplings and Excitation Paths 

Figures(V-3) and (V-4) show some of the more prominent 
couplings linking lower levels and associated energy gaps. 
For H*-H(1s) collisions in the energy range considered, the 
strong 2pey-2pTlu rotational coupling associated with united 
atom orbital degeneracy is the main source of excitation and 
excitations from the entire g-manifold , or from 2pd@ via 
radial couplings, are less important. Thus for example, a 
good qualitative picture for the 3pa excitation is the two 


step process 2pé> -2p7m -3pé> . However the radial 2pq- ~3pd> 
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Table V.1 4,8,10 Molecular basis sets 


Molecular Basis Sets 


for systematic close-coupling calculations 
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states basis states basis states 
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Table V.2 16 u-State basis sets 


Molecular Basis Sets (u-states only) 
for upper state excitation/convergence studies 
Atomic 


Basis States principal 
quantum # 


# of u 
states 


|2po,>s 
|2pn >, |3po>, 
[3pm >, [4fn >, 


[4pm >; 


8 above, plus 


[4f5 >, |5f0 > 
u u 


[Stn >, |5po,>, [5£6 > 


[5p > 
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u 


|6pn > 


a 


| 
-— = eager - 


a 


——<— vee 


[ (lee avte?e i) ‘aa 1? eine person 


orhuja sateg? varins\aehtan tains. sexe coger wot if 
1 Hm wichita sais i 


1 3a98 adend / 


a 


rey ae 
Leet mii ow 


“ety ney, i me i; 
“aia hall at nel 
uf al evan 7 eee 


:% 


a a enn ll ieresiiies 
4 

ae 

a 


ake ies haat 


“tn ase: . 


= 

. 
ee ee 
; 

@ 


, nant 
we we ea 


— ee a 
a 
i. 


134 


Electronic Energy (Ry) 
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Figure(V-1) Electronic potential energy curves, €(R) vs. R, 
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Figure(V-2) Electronic potential energy curves, €(R) vs. R, 
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coupling is not negligible and at higher energies it 

plays an important secondary part in the 3pf% excitation and 
resulting atomic (2s,2p,) cross sections. 

In the g-manifold, the most significant couplings from 
Is dp are radial 1sé@> -3dd> and angular Is@> -3dilg , but due 
to the large energy gaps, neither is very effective except 
at the higher energies studied. At small impact parameters, 
radial 185 —2s0, coupling plays} Some part;} but it is 
unimportant at the intermediate and larger impact parameters 
which make most contribution to the cross sections. Which 
initial process is more important in the g-manifold depends 
sensitively on both energy and impact parameter, but in any 
case the strong rotational 3d Uy ~30 Mg coupling populates 
both states in a second step. Because the contribution of 
the g-manifold to overall jexditation is smallJ most/ of our 
discussion (and study) focusses on the u-components. 

Figure(V-5) depicts in a qualitative way the main 
pathways for excitation of the 5 lowest molecular states of 
each symmetry: thicknesses of the arrows indicate the 
effective coupling strengths. Note that both angular and 
radial couplings to 4f% are comparable and relatively weak 
( a different result than is found by Crothers and Hughes 
who have a stronger 2p%y -4£% coupling). In the g-manifold 
we can see that couplings among excited levels (e.g., 3d oq 
-3dig) are as strong as those in the »u-manitold, “but tie 
difference which makes g-manifold excitation unimportant is 


the very weak initial links from 1s@¢ . 
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Figure(V-3) Prominent radial and angular coupling matrix 


elements and corresponding energy separations, vs. R_ for 


u-manifold. 
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Figure(V-4) Prominent radial and angular coupling matrix 


elements and corresponding energy separations, vs. R for 


g-manifold. 
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As noted in Chapter III, there are many strong 
couplings among the more highly excited levels; the states 
involved are diffuse and easily polarized, so that couplings 
Mave quite Tong = range, ‘aSm@main the radial 7-71 couplings 
shown in Figures(III-9). Our study of upper (ungerade) state 
excitation paths has not been as systematic as for the lower 
States, but Figure(V-6) indicates in an approximate way some 
of the important flux pathways leading to upper-state 
excitation in the u-manifold. The 7 - 7 couplings lead to 
upper-state excitation through the sequential neighbouring 
np my ~(nt+1)pNy radial coupling. Some flux is promoted to the 
np@s states from the initial 2p@>\state through the radial 
Gourling, but »thel strong nog Snpyivpmangular “couplings 
populate the flux back into 7Jy states. The role of nf Tu 
States through radial and angular couplings is relatively 
weak and hence no significant contributions have been 
observed for upper-state excitations by this route. 

lt may be of interests ito” pointiout again that the 
upper-state g-manifold couplings (e.g., 3d ~3d a ,etc.) are 
just as strong as’ those “form, the gu-States, but play no 


important role only because of the weak initial step from 1s 
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3. Asymptotic Couplings within the H(n=2) manifold 
Within the asymptotically degenerate (2s,2p) manifold, 


long-range radial and Coriolis couplings have a_ significant 
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Figure(V-6) Important flux pathways leading to upper-state 
excitation in the u-manifold; thickness of connecting arrow 
indicates qualitative importance and dot-line indicates weak 


couplings. 
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effect on the detailed atomic state populations. Both the 
couplings and the splittings between coupled levels decrease 
as R°’*. Beyond R 220 a.u., however, the n=2 sublevels are 
effectively decoupled from levels not degenerate with them, 
and for each symmetry there is therefore a_ three-state 
asymptotic coupling problem, which we have found convenient 
to treat separately. The three states involved for the 
u-manifold are (3pé@> ,2pmu ,4£%> ), and for the g-manifold 
(300g 301g , 2805 Wednmeachsicase the two 0 =<states are 
sp-hybrid atomic states and are linked to the 7 -state by 
Coriolis coupling and to each other by radial coupling. For 
R 2 15-16 a.u. we found we can describe both couplings and 
splittings accurately by analytical models using 
atomic-state perturbation theory. By introducing an 
appropriate action variable aS a progress variable (instead 
of R, or time t), the resulting asymptotic coupling problem 
can be solved efficiently. In effect, its solution may be 
represented as a three-state propagator U(oo;R ) which 
converts molecular state amplitudes at R, to those at R00; 
This is attached to the output from the integration over the 
interior "real collisions" region. We used this device to 
generate amplitudes for all three final molecular states, 
even in those cases where the basis used in the interior did 
not span all those states (e.g. in the 8-state basis 4£%7 

and 2804 do not appear in the close-coupling problem, but 
amplitudes for these states at R-oo are fed from amplitudes 


in the other components (2pfly ,3pfu ) and (3dM%g ,3d% ) at 


sat hae ie 
Wie! ce yaya 
% | 


wz &4ot .enokgeluged wel 


javideh alent Galquen’ peswiod 


spaiteilqe is ae 
oom  elowedciae. Sieg, ade vevesweut: er Oss a bre: . 
marc ot diiw etaremwgeb tom, — oh ‘Be lqueseb ievida 
atayereesdy ss stelswedd we, seta; Yatemmye foee a 
iheiwevien Srwet vad as doidw iaeidorg pntiquos had 
sii? 20a few Lovad | vias ‘ait  dednasqna’ 

fnmrg ate a02 bag, | aa fi Qh, eet sah las 

eovere> 7% wail ait oa. dose or ft pas, gine 

gin KH wil 49 haanad ond Bi esvede, ahmed ‘pia - 
ist .patinwos bethany ws rete. Lila od ‘ore gaitquess LOkT 
bie agedilqwos Aged sd xoaat bu ow H bine Lud ate ater. 
Suen aba ai Coney | teas taunt gd Roe: igi 
& el Mines ai 4 “ynostts 3 ; 
Sie they d | dadvavr aediparg Sj an eiatvad nobipe saat! 


ae 


a 


‘ip 
an ia 


alee oni Lawes ke ee eatstume, oda, 144 wats 10 al | 
oa Ven nol Hikoe ars 2 catia oh ehangiaiire hevios od 
estiw 4. Be oe7g a epee 2g’ sina -weds 6 @s hetoent 
1 Od Je otoade a3 oh $4. rabad tiga s2e32 eda Siam aa4 | 
gid 3ev6. sol+arpomini wd neal Mugave add of baioesse a 
o2 ppeveh adds iia one tok pe "anoteilioo Ise”. rolseind : 
8osnt2 yvalvoedion bends sont de iin nod aebus i lqus orers09g | - / 
bth yoiuisnt ade oi nees abel efi drei peres saody nt neve ie 


| ah 


td ziend otasar) wd WP 4 Bea) aeteva etod? Ifs fncre ton: y 


re 


wd «6 waldo nw, dnékquestweds’, Say ai tesqas son ob pes bas 
estuttique wer! bel esa aeRO @eveve seedd yo? eebus! fqms 
28) 4. Geet RE) ban 4 THE. Gites) asAsnogmes Yadi¢ ad? wh, 
Can ea ; ; F 
: vt | ae 


143 


R=R, ). Effects of the asymptotic coupling are small (<10 % 
changes) but are certainly not negligible. This should be 
taken into account to give predictions of atomic 2s,2p, and 
2p, excitation probabilities. 

Section D. in Chapter IV gives a detailed account of 


our procedures for treating this asymptotic coupling. 


4. ETF Effects on the Couplings 

In this section we compare our coupling matrix elements 
(with ETF corrections based on switching functions, as 
described in Chapter III) with those employed by Crothers 
and Hughes who used a different treatment of EFAs: 
(Comparison with uncorrected PSS matrix elements is not 
considered here since it is now generally well known that 
there are very subStantial corrections in most cases; the 
papers by Crothers and Hughes give specific examples of such 
comparisons with their corrected couplings for this system). 
As pointed out in ChapterIII(see also TKCK(33)) the ETF 
description used by Crothers and Hughes is equivalent, to 
lowest order in the collision velocity, to the use of a 
particular type of switching function. For the isd state, 
their effective switching function closely resembles the one 
used. there extR2=2c0 eracuals sbuticior . the. 2p75  state.1t ws 
qualitatively very different (see TKCK(33)) and leads to 


quite different coupling matrix elements. 
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As relevant examples, Figures(V-7,8 & 9) show our 
Souplingsmatrix elements for 2p¢; -2p7). » 2D e4io” ,2D Tu 
=4f¢@y , and 1s oq -3d0, , in comparison with those employed 
by Crothers and Hughes (Hermitian averages). The asymptotic 
behavior of the couplings is the same, but in most cases 
there are substantial differences at finite internuclear 
separations(s10-12 a.u.). 

As we have found by actual calculations (see Section 
V.C below) these differences in coupling matrix elements can 
lead to significant differences in excitation probabilities 
and cross sections. For the 2pé0) -2p7y coupling, the node 
and steeper slope in the matrix elements of Crothers and 
Hughes leads to a considerable increase in the 2p7y 
excitation cross#section ( \tog first order inv’ velocity), 
relative to our values. The much larger 2p@y -4f% radial 
Coupling of Crothers ang’ HugWwes 1s "mainly presponsible for 
discrepancies between our results and theirs on the detailed 
magnitudes of 2p, and 2p, cross sections, the resulting 
polarization of Lyman- @ radiation, and (to a_ smaller 
extent) the total 2p excitation cross section (see Section 
C.2. below). Similar effects might be expected to result 
from differences in the Is0g ~3d0q couplings, but of course 
these are much less important because of the much smaller 


g-manifold excitation cross sections. 
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9) Coupling matrix elements linking 2p 


O> .2pt 4f0u states, and 1s0g ,3d0q states, as computed 


with analytically determined switching functions (solid 


curves--note non-Hermitian character), and as computed by 


Crothers and Hughes(48) (dashed curves--Hermitian average). 
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B. EXCITATION PROBABILITIES 


1. Collision Histories 

Fagures(Ve1i0°& 11) show “two “aypical “histories” ‘of 
collisions (molecular state probability vs. time t) for the 
5 u-state basis at E=5 KeV and impact parameters f =0.5 and 
5.0 a.u., respectively. Early in the collisions some effects 
Seene Vongerange. tail, in the 2pdg" -4f@, coupling’ can be 
seen, followed by 2p¢g> -3p@ radial coupling and (at small 
impact parameters) the dominant 2p@> -2p7iy coupling. Later 
execitavicn jOft 3pi/lu Vviev-frotational coupling), fromespd,. can 
also be seen at smaller R. Finally during the outgoing part 
of the collisions there are effects of 3p7fy -4f£4 and 3p7T7vu 
-2pTlu couplings. The sudden dip of 4£%) probability around 
t=0 for Ff =5 a.u. is probably ‘connected with the 3pm, -4f0y 
curve crossing at R=6.5 a.u. and a smaller effect at both 
earlier and later times may be similarly related to the 3p 


Usecetc, crossing at 'osavu.- 


2. Excitation Probabilities vs. Impact Parameter 
Figure(V-12) shows molecular state excitation 
Probabilities vS. Ampact. ~parameter Gat \) 5, ‘Kev * for the 5 
u-state basis. The importance of radial couplings is 
illustrated clearly by comparison with Figure(V-13), which 
shows the same data when only angular couplings are 


retained. The most obvious effect of radial coupling is the 
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Figure(v-10) Collision 


probabilities vs. time (a.u.)) for E=5 KeV, impact parameter 
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Figure(V-11) Same as Fig.(V-10), but 
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Figure(vV-12) Probabilities for excitation of molecular 


States vs. impact parameter at 5 KeV. Angular and radial 


couplings included. 
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Figure(V-13) Same as Fig.(V-12), but only angular couplings 


included. 
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change in the 3p@ excitation probability, especially at 
small impact parameters, arising from direct excitation 2p 

OS -3p0y> . A secondary effect is the 5-fold decrease in the 
3p My probability which occurs when radial coupling is 
included: 3p)[y amplitude, produced at small R by angular 
S0upriIng, Leturns tO 2pyp dace in the*cotlision due ‘to™*the 
long-range 7] - 7i radial coupling. The effect of the direct 
radial 2p05 -4£ 07 coupling on the 4£07 probability can also 
be seen, though this is much smaller and mainly influences 
the shape rather than magnitude of the 4fO0> curve. 

The effects of radial coupling are also very evident in 
Figure(V-14), which compares atomic state excitation 
probabilities vs. impact parameter at 1 KeV for the states 
2p,,28s+2p,., and 3p.,, with the results of Schinke and Kruger 
who neglected radial couplings (and also ETF corrections). 
The dominant 2p, (2p7ly ) probability is not much affected 
and hardly even changes its shape when radial coupling is 
included, but the (2s+2p,) probability is very different, 
Gue-"to “our inclusion “ol -™ene 2D7> spe, Girect couplings 
Also, our 3p7v probability is significantly increased at 
larger f , relative to that found by Schinke and Kruger, 
(though the absolute magnitude is in any case much smaller). 
From these comparisons, we can see that there is no 
theoretical justification for neglecting the radial 
couplings; although their matrix elements are smaller, they 
act effectively over a wider range of impact parameters and 


become particularly important for increasing collision 
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energies. 


C. CROSS SECTIONS 

1. Total and Differential Charge Exchange Cross 
Sections 

Charge exchange in H‘*-H(1s) collisions at these 
energies is dominated by the resonant charge transfer 
process. Figure(V-15) shows the H(1s) charge transfer cross 
section vs. energy together with the H(2p) and H(2s) charge 
transfer cross sections; the contribution from H(n=2) charge 
exchange is only 3.5 % even at the highest energy studied 
here. Because of this dominance of the resonant’ process, 
nearly all theoretical results, both from atomic-state and 
from molecular-state basis calculations, are in _ reasonably 
good agreement with each other and also with the 
experimental values reported by McClure(63) and by Fite, 
Smith and Stebbings(80) (below 3 KeV, the results of Fite et 
al are somewhat higher than those of McClure). The flat 
portion of the cross section slowly increases toward a 
maximum in the neighburhood of 1 KeV. Uncorrected PSS 
calculations are quite sufficient to obtain a good account 
of the total charge transfer cross section in this’ region, 


as has been shown by F.J.Smith(41). This is easy to explain, 


since 
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Figure(v-14) Atomic state excitation probabilities vs. 


impact parameter at E=1 KeV. Solid curves, present work; 


dashed curves, results of Schinke and Kruger(47-a) (angular 


couplings only, no ETF corrections). 
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(Cm) 


CHARGE EXCHANGE CROSS SECTION 
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Figure(V-15) Charge exchange cross sections for H(1s),H(2p) 


and H(2s) vs. E: ,present work. Experimental data: 


@ ,McClure(63); @,Fite et al(80) for H(1s) charge exchange 
cross section. For H(2p) and H(2s) charge exchange cross 


sections, data are the same as in Figs.(V-17) & (V-20). 
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(1) the PSS results are exact within the two-state (1s 
Og ,2pG>- ) approximation where only elastic scattering 
occurs, and 

(2) below 1-2 KeV the dominant excitation occurs via the 
Stronge2pdy, -2p No rotational coupling, and) its not? ‘much 


affected by ETF corrections below 1 KeV. 


Our results for the total charge exchange cross section 
are therefore in good agreement with other molecular state 
calculations and comparison with these other results is not 
Particularly significant. (All our “numérical results for 
direct and charge exchange cross sections are tabulated in 
Tables(V-3) and (V-4)). 

We have also computed the differential cross_ section 
for charge transfer into all states, using the small-angle 
formula given in Section(II.D) (Eq.11-42). This is ncun vs. 
scattering angle 9 for E=1 KeV, in Figure(V-16). Here also 
the dominance of resonant charge transfer and the 2p(@> -2p 
Tl, excitation process at low energies means that the 
calculated results from nearly all molecular state studies 
should be in good agreement, especially at small angles, and 
this is confirmed by the very good agreement between our 
values and those found in the three-state calculation by 
McCarroll and Piacentini(44), for @ $2°. At larger angles, 
however, our results show improved agreement with 


experiments by Helbig and Everhart(60) and by Houver, 
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Table V.3 Charge exchange cross sections 


Atomic Level Charge Exchange Cross Sections 


H* + H(ls) (10-State Basis) 
(cm ) 


en (KeV) 
aan 0.637 x 10 + 0.761 x 1077 2420-X 410" 
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Table V.4 Direct excitation cross sections 


Atomic Level Direct Excitation Cross Sections 


H’ + H(1s) (10-State Basis) os 


uae H (2s) H(2p)) 
Energy (KeV) 
ar: Je | Te eee 
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Figure(V-16) Charge exchange differential cross section (all 


States) VS. lab. scattering angle at E=1 KeV. cammmmmum 
,present work; —<—|$——= Schinke et al(47-a)3; ececccces 
McCarroll et al(44). Results of Crothers et al (48), not 


shown, nearly coincide with ours. Experimental data: o 


Houver et al(61); A Everhart et al(60). 


; 7 - 7 7 oe Renny La; cat) 


ie 


gar 


toad eae 


a vi i" 


- 
eS 


7 a ~ 
- a vy ine : : f sn 
ims : al a _ ; a7) 
F iv : z = a 


<< 


(ig) mobedea sted bisnaeii i ; senitons nse 01-5 wn ups 


| eel ‘VoR P92 20 ofens ‘eetvegsage a his: 
; | hesvene (este) ie ga sanition “a poy 1O¥ 
“ton ae) boas viadson a9 eniueat ) iim, go bk 
| “© peers sarqemi tees _ ate 


aan sy: ; ij 

) i t G 
Sabi 7 
ars i 


161 


Fayeton and Barat(61), relative to the three-state results; 
this was also found by Crothers and Hughes, whose results 
(not shown) essentially coincide with ours. The improvement 
is mainly due to the inclusion of the radial coupling 
ettects: {including ETF ‘corrections). ‘Ber still larger 
angles, the approximations involved in the small-angle 
formula (and the use of a straight-line trajectory) are no 
longer valid and Eq.(I1-42) should not be expected to give 


results in agreement with experiment. 


2. H(2p) Cross Sections 

In) Figure(V-7) our |} 10-state results for the H(2p) 
charge exchange cross section (2p,+2p,) are compared with 
results of Schinke and Kruger and Crothers and Hughes, and 
with the experimental values reported by Morgan, Geddes and 
Gilbody(66) in the energy range 1-7 KeV. For 2p excitation, 
the contribution of the g-manifold is relatively small, 
hence direct and charge exchange cross sections for H(2p) 
differ by less than 8 % over the energy range studied. The 
(included) corrections due to second-order terms in velocity 
(see Section(V.D) below) are at most 16 %, at the highest 
energy studied. The agreement of our results with experiment 
is generally very good, and is somewhat better than that of 
the other theoretical results, though these also are in 
general accord with experiment. [The success of molecular 


State calculations in predicting correctly at least the 
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Figure(V-17) Charge exchange cross section for H(2p) (all 
components), vS. E: mJ: present work; —==j—m_seme: Crothers 


et_al(48); ----— Schinke et al(47-a). Experimental data: 


O Morgan et al(66). 
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qualitative behavior vs. energy for this cross’ section 
confirms that they are appropriate at these energies; 
atomic-state basis calculations (not shown) give much larger 
(more than a factor of 3) cross sections which vary rapidly 
at low energies and predict the location of the maximum 
cross section inconrestlyga adjthoughsisthest inclusion prof 
pseudostates(Cheshire,Gallaher and Taylor(57)) causes marked 
improvement in the results]. 

The discrepancy between our results and those of 
Schinke and Kruger can be shown to arise mainly from the 
neglect of ETF corrections, and (to lesser extent) from the 
neglect of radial couplings. (The curve shown for their 
results was constructed by addition of their published 
curves fotfae2py:r ands fork Gsm Gegquads —toiel2pw. , inv iitheiz 
calculation)). At 5 KeV, for example,’ more than 90 % of 
their total! |H(2p)recross section qarase’s p from). the 2p, 
excitation via angular coupling from the initial 2p@> state. 
As Crothers and Hughes also point out, the uncorrected PSS 
2p0> -2p7y coupling matrix element behaves in a completely 
incorrect).way _ for. R230a.u., due to the, linearly increasing 
"moment-arm" term in the uncorrected operator ag ; 
Calculations made by Dr. V. SethuRaman and Professor Thorson 
(unpublished) in this laboratory indicate that above 1-2 Kev 
the resulting spurious coupling leads to _ significantly 
larger 2p7ly excitation probabilities, as were found by 


Schinke and Kruger. 


EBt: 


nobavee enor eis 


jastptans sents of 


wep iel anu 44 co iineg 


vibiges yey doisiw enolsoeg swore, te te pe ® sad 8 | 
eiinixen oft t6 °° sofsacol sae: soiheang bie aslex 9 4 
> dolevwfoat ef? Rt te qtsoerienni sneha am 
bad am seannn (£72) oad te etal dehenedoh ona 
Soe alah a oh 

sodA ine ailsers te nedwise edie 
mort gli aw Batam, Ber pilik ae 1 soot ml 


fis. mo73 (7nees Treane) oy ae aon r e 
fords m4 aweste ev are ont) a" 


ry 


SKtieilidaa teed he ietsibke 9 . a 
tisf3 i «yrs O29 Cage). a8 coe ake. 
2 6flCU OP? CU? 6 aoe Dalat daa ws an 
ot od2 rrr Mag iis awtengs deeb; nee 
Sigta Guys Leadial ets aon Bit) Seren: wake « 
a2" ied pet4oonu 98% |, tuo snpog nots panded bre co 
vivtwignion « Ab gevedad dremel kegs Paitquon 
pRiaegs s04 tines d sacl a 4 abe ad 104 uw 


; <a seg eaete) Gatcer’ 5 Sani add ni pants + 
noaios? rogervter? Soe Viemiodek ov cae va obese anglss 
Ys) S-] qvode tadd sdeathnd adtevodal abs at pane 


yiweeeitaneta of ebay) gatiquan pore entstinty oat 
‘@ Boucl arad 24 epbdhibdadony nohtas Lae oh a3 


164 


The results of Crothers and Hughes are in _ good 
agreement with ours up to about 6 KeV and then show a_ rapid 
increase relative to ours and to the experimental values. 
Calculations we have made (see sub-Section(D.4(b)) below) 
indicate that this increase is due to the much larger 4f% 
populations they computed on the basis of their ETF 
description and the resulting coupling matrix elements 
(Figures(V-7)-(V-9)). Our own calculations (Table(vV-3)) show 
that the 2p, and 2p, charge exchange cross sections are 
approaching their maximum values near 6-8 KeV and that the 
H(2p,) charge exchange contribution is at most 25 % (at the 
highest energy shown). For the direct excitation cross 
sections, the 2p, cross section (as noted earlier) is only 
Slightly different from that for charge exchange, with a 
maximum near 6-8 KeV, while the 2p, direct excitation cross 
section is still increasing (maximum near 20 KeV ?); at the 
highest energy studied here, the 2p, contribution to direct 
He2p) rekqitation tsestwklconlySee%lcof sthe)letotal:.) These 
findings are in contrast with those of Crothers and Hughes 
who obtained 2p, cross sections larger than those for 2p, at 
these energies, again directly related to larger 4£% 


excitation probabilities. 


3. Polarization 
The polarization of the Lyman- A radiation emitted by 


the inelastically scattered atom can be calculated by using 


Syap ni SIS Ste 


bigea & soe Wes 


(eotad (va) d.ai16 
Hit see vet Agem ais oy wih Br * semen ai die 
ETF iohe 26 steed did ae _ bepuqns Lanes 
at¥uewala -eit2ed ers Lae “og | int themay, aha line, ; 
votes Uls-V)elidat) ane 5 etn 
gue lsoée 22e%0 Perasdiowe att ae ie 4 
eis fade Bite Vea Beal tebe outa aed na i wt 
oft 4a) ¥ BE doom Fe ei ross Anthea veiadoxe oan 
Betis Not eSs (Ike $9078 td ht Aneta ‘agen i 
Gf) Oe) dort bse avons, = cml : oil 
n “ake spnaiaxe soteds. a8 ted ibe 4) i} ; ” g fs 
zeuts Hbicds tons 4ed2he ae eae paint (en 09 008 
et2 #8 (6b "eet OS Réee SER homies fide od 
soeu lh os weitadi Avs oge os \ acti deh base ieasne | 
saent iveteyn |g ee ome xine fee ‘ya - 
vatoh< bya etied sna = events ives abd toma wie: ate apn 
te.QS 309 asiegtongds pas gab lane whet as bentaadg 6 
Mat weprel ot bharalet rere. wings auigrens 
| i" aap Heimer ne | 
coral | 
ne Paes 
aosgactaten 2k 3 
qa betzins noizeiia: A -neeyD ede Yo cotdagientoe ‘edt | 
enienw 70 Seveluaie>s sf as mote Baned sane gilestseeteni ei? 


iu 
uP 
= 
—~ 
o4 
ZA 
e 
~r 
th 
—. 


165 
the formula derived by Percival and Seaton(49) 


D 
Q°(2p9) - 9 (2P)) 
I = ine te Dosh 
aQ (2p,) + 8Q (2p,) 


where QO =2.375, B =3.749. Results are shown in 
Figure(V-18) together with some other theoretical 
calculations and one set of experimental values (Kauppila 
et _al(69)). Since our excitation cross sections for 2p, are 
much larger than those for 2p, over the entire energy range 
considered, we obtain consistently negative values. These 
disagree with those of Crothers & Hughes, which become 
positive between 4 and 5 Kev. The sign change in their 
results is directly attributable to the much larger 4f% 
excitation probability they predicted Gue to their 
description of ETF effects. 

Actually, all sorts of numbers between +1 and -1 have 
been reported in the various theoretical studies, using many 
different approximations. Only one set of experimental 
values iS available to test these predictions and 
unfortunately its reliability is not really adequate even 
for a qualitative conclusion which might reduce the 
theoretical chaos. It would clearly be very useful to have 


more reliable data on the Ly-QA@ polarization in this energy 


range. 
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Figure(V-18)Polarization of Lyman- @ radiation vs. E: 


present work; =—_—_/=—_—=— _,Crothers et al(48); 


—— ,Rapp et al(58); —~— ,Gallaher et al(57). 


Experimental data: O,Kauppila et al(69). 
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4. H(2s) Cross Section 

(a)Cross Section 

Cross sections for direct and charge exchange 
excitation to the 2S atomic state are shown in 
Figures(V-19 & 20). Effects of Isdy -3d@%- and other 
g-manifold couplings play some role at higher energies 
so we find that One and Or are Slightly different. 
For charge exchange, our results agree rather well with 
the recent data of Morgan, Stone and Mayo(68) and of 
Hill, Geddes and Gilbody(67) from Hoo) ton6 Kev (though 
Hill et al give somewhat larger values than those of 
Morgan et al below 3 KeV); Our values seem to agree more 
closely with those of Hill et al in this region, and 
fall somewhat below both experimental values at 7 Kev. 

On the other hand our results differ very much from 
those reported by Crothers and Hughes, which show a 
pronounced minimum near 3-4 KeV in both direct and 
charge exchange cross sections. This behavior is similar 
to results found in some atomic basis set calculations’ 
(not shown) and appears to follow closely the older 
experimental values reported by Bayfield(64). 
(b)Calculations using Crothers and Hughes matrix 
elements 

In an effort to determine the origin of this and 
other discrepancies between our results and those of 


2The atomic calculations (Cheshire, Gallaher and Taylor(57) ) 
show a change of slope, but not such a deep minimum, as that 
found by Crothers and Hughes, 
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Figure(V-19) Direct excitation cross section for H(2s) vs. 


—m— Crothers et al. 
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Experimental data: O Morgan et al(66); 4 Chong et al(65). 
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Figure(V-20) Charge exchange cross section for H(2s) vs. E: 


present work; ——_———,Crothers et al(48); 


——=—— Schinke et al(47-a). Experimental data: @,Morgan 


t al(68); M,Hill et _al(67); a,Bayfield(64). 
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Crothers and Hughes, we have carried out close-coupling 
calculations using the same 10-state basis, but 
employing the matrix elements shown in their papers, for 
E=2,3,4 and 5 KeV. Straight-line trajectories were used, 
and we retained only terms of first-order in velocity. 
The essential features of the results we obtained are 
Shown in Table(v-5). For all except the 2s excitation 
cross sections, we have obtained reasonable qualitative 
agreement with their published curves. For example, the 
2p, excitation cross sections we obtained in these 
calculations differ from their reported values by at 
most 20 % (at 5 KeV) and follow the same uniform trend; 
the discrepancy in this case is probably attributable to 
the effects of higher-order terms in the velocity, which 
are included in their reported values, as well as in our 
own calculations. However, in the case of the 2s cross 
sections, we have not been able to reproduce Crothers 
and Hughes' results even qualitatively; the results we 
obtained from these calculations lie about 30-40 % above 
our curves in Figure(V-19) and show no sign of any dip 
or even a marked change in slope. We therefore cannot 
prove that this large discrepancy necessarily results 
from the different treatment of ETF effects by Crothers 
and Hughes. 

From remarks made by these authors and comparison 
with our own experience, (Section(V.E.1) below) it 


appears that their second-order velocity corrections are 
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Table V.5 ETF effects on cross sections 


ETF Effects on Individual Cross Section 


First-Order in V Only (units of nD 
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rather larger than those found in our approach. As shown 
“PraSectron( Wir) etsde ‘jalso! irKeK(33)i):, the ETF 
description employed by Crothers and Hughes is 
equivalent to the use of switching functions only up to 
first-order terms in V. Therefore it is possible that 
their second- and higher-order corrections behave quite 
differently from ours, and we cannot exclude the 
possibility that such second-order terms are the origin 
of the 2S cross section discrepancy. In our own case, 
however, we find that the second-order corrections are 
much smaller than would be needed to do this, and they 
are also steadily increasing functions of the energy 


(ctacSection (Veh below). 


(c)Comparison with Charge Exchange Results of Schinke 


and Kruger 
Figure(V-20) also shows the results for the 2s 


charge exchange cross section reported by Schinke and 
Kruger; the results found by Chidichimo-Frank and 
Piacentini (not shown) are very Similar. There is a 
substantial discrepancy between our results and these, 
which is masked to some extent by the linear plot used 
for Figure(vV-20); at 1 KeV our value is at least 4 or 5 
times larger than that of Schinke and Kruger and at 5 
Kevyit @isestillvaboutees0 % larger. As €use"shown thy 
¢ompakisonwertofesrigures(V-12) aandi(Ved3) (5eKeviymour 
larger cross sections result from “Inclusion of (| sthe 
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higher population of the 3p%> level at smaller impact 
parameters. Figure(V-14) (1 KeV) shows that this effect 
is even larger at lower energies. Given the limits of 
error on the experimental values reported in Hill et al 
and Morgan et al, we cannot’ say that Piece clearly 
confirm our results as correct, but they seem to agree 
somewhat more closely with our predictions than with 


those made by Schinke and Kruger. 


D. VELOCITY AND TRAJECTORY EFFECTS 
1. Effects of Second-Order Terms in Velocity 
In this study we have included the effects of 
second-order terms in the collision velocity where these 
could be shown to have a significant effect on the cross 
sections. In this section I present data showing the 
size of these corrections for various cross sections as 
a function of the collision energy. Table(V-6) shows the 
comparison of the cross sections calculated using the 
first-order terms only and by the inclusion of up to 
second-order terms. The effects of the second-order 
terms can not be seen clearly below 3 KevV; however , 
they become significant as energy increases. At the 
highest energies studied, the second-order contribution 
rises up to  1/° % for some individual cross sections. 


Therefore in H*-H(1s) collisions, the second-order terms 
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Table V.6 Effects of second-order terms in velocity 


Effects of Second-Order Terms in Velocity on 
Charge Exchange Cross Sections 


(units of — 


ei 2s 2Po 
[ci ges En ie 
Cea a 


ove 


(77P30L4y HE eMIeg 


2. 


175 


play an important role, (as I already pointed out the 
Sizes of the matrix elements are relatively larger at 
large distances than those for the asSymmetric case. (see 
Section(I111.D))). 
Effect of Coulombic vs. Straight-line trajectories 

I have also made some studies at low energy of the 
effects of using a straight-line trajectory, by 
comparing with calculations made using Coulombic 
trajectories, and these studies are summarized here. 
Table(V-7) compares the effect of the two trajectories 
at 1 KeV. The discrepancy between the two different 
trajectories employed is at most 4 % in the H(2p,) cross 
section, due to the creation of the new peaks at small 
impact Parameter in the Coulombic trajectory 
calculation, and in H(2s) and H(2p.), the effect is much 
smaller, about 2 %. Therefore we can conclude that a 
Straight-line trajectory 1S adequate to the present 


purposes. 
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Table V.7 Trajectory effects on cross sections 


Effect of Trajectories on Atomic Level 


Charge Exchange Cross Sections 


+ 
H energy E=l Kev (units of cma) 


Trajectory 


Straigh-line 


| 1.98 X 107. "| 0.637 x.107-° | 0.761 x 1007° 
| 2.06 0.625 0.744 
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E. CONVERGENCE STUDIES 


1. Convergence of n=2 Levels 

Tables(V-8) and (V-9) shows excitation cross sections 
for H(n=2) levels as a function of basis size, for the basis 
sets of Table(V-1) >, and for energies 1,3,5 and 7 Kev. In 
all cases the H(2p,)cross section appears to have converged; 
on the other hand, for the-H(2s)cross~sectionm there is a 15 
% change at 7 KeV when the basis is augmented from 8 to 10 
States, and similarly a 20 % change in the H(2p,) cross 
section, so in these cases the effects of further increases 
in basis size are worth further study. Unfortunately no data 
on convergence vs. basis size were provided in the study by 
Crothers and Hughes or by earlier calculations on this 


system. 


2. Convergence and Flux Loss to Ionization 

To explore further the convergence properties of a 
close-coupling calculation on this system, we increased the 
size of the u-basis by adding more highly excited 
molecular-states strongly coupled to those already 
considered, and performed calculations of molecular state 
excitation probabilities at selected energies and impact 
parameters. In “additions to the. -S<state § basis “aiready 


considered, we used basis with 8,14,16 u-state. Table(v-2) 
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Table V.8 Convergence of cross sections 
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Table V.9 Convergence of cross sections 
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lists these states and the corresponding atomic. state 
manifolds to which they dissociate. As these Rydberg states 
are very loosely bound and closely spaced in energy, they 
can hardly be expected to behave even approximately in 
adiabatic fashion. Moreover the coupling matrix elements 
linking them, both radial and angular, are of very long 
range and reflect the very large polarizabilities of these 
states. Even at the low end of the energy range we 
considered, it has proved impossible to obtain reasonable 
convergence of excitation probabilities for atomic levels 
n=3. Moreover, as the baSisS Size 1S increased a small but 
non-negligible portion of the flux persistently moves up to 
the highest levels accessible; a collision history diagram 
in Figure(V-21) shows that this is occurring in the outgoing 
portion of theg collision “trayectorys and wise due to the 
Jong-range 7 =- “Tl! couplings £ane particular; “Table(v-10) 
indicates the size of the effects in question for E=0.7 Kev 
and three impact parameters. It can be seen that a small 
portion of the flux always escapes to the highest Rydberg 
states as the basis size increases, while the probability 
for n=2 levels is relatively stable for all impact 
parameters. However, aS Nn increases, the probabilities 
change by more than 30 % for each impact parameters and this 
is particularly so at large impact parameter when the basis 
is expanded from 14 to 16 states. 

Since close-coupling molecular state calculations, as 


currently performed, include only bound states which follow 
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Table V.10 Excitation probabilities for higher n levels 
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the nuclei, they make no allowance for flux loss due to 
ionization. However, aS waS pointed out by Thorson & 
Levy(8), it is physically unreasonable to regard a state as 
"bound" in a collision, if the transport kinetic energy of 
the electron, relative to the other collision partner, is 
comparable to or greater than the static binding energy; 
events corresponding ina tightly bound state to charge 
exchange will lead in such a case to ionization. At E=0.7 
KeV, levels with n25 fall in this category. This suggests 
tietiwiee sagnaficant “portion ‘“ofi® the™probabildaty. l4a'sted™in 
Table(V-6) for the high levels corresponds in reality to 
ionizing events. SethuRaman, Thorson and Lebeda(8) and 
recently Choi and Thorson(8) have calculated the direct 
ionization from the close-coupled 1S 0q ,2pG> and 2p7u 
molecular states of H* for projectile energies E= 1 KeV, 
assuming that ionization dominantly occurs by a Single 
impulsive excitation from a tightly bound level to the 
continuum. The results obtained show that at 0.7 Kev and f 
Sa eO aa s:lie70) Ce resulting u-state impact ionization 
probabilities are 50-100 times smaller than the aggregates 
for n=5 and n=4 listed in Table(V-10) and drop off much more 
rapidly with increasing impact parameter “than, those in 
Table(v-10). 

From reasonable extrapolation of the experimental data 
of Fite, Stebbings, Hummer and Brackman(62) on ionization, 
the experimental values should be about 50-100 times larger 


than those theoretical results and give nearly the same 
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order of magnitude as are found here due to the "ladder 
climbing" process. 

In conclusion, we may expect that a small but 
non-negligible portion of the flux leaving the ground state 
in H*-H(1s) collisions, even at energies ~1 KeV or less, 
will be promoted to higher (Rydberg-type) states by the 
multistep "ladder-climbing" process seen above, and that a 
large portion of that will actually be ionized in the end. 
We may also expect that close-coupling methods which do not 
mciude flux loss” to the continuum as ‘part. of their 
formulation cannot represent the excitation to upper levels 
correctly or convergently. Further theoretical study of the 
"ladder-climbing" ionization mechanism in this prototype 
system is obviously highly desirable to further understand 


Slow collision processes. 
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VI. He?*-H(1s) AND H*-He*(1s) COLLISIONS 
Pn tnis chapter we present and discuss our calculated cross 


sections for the processes 


(a) Het > + Hi 1s) ——— > He**+ H(nl) (DIRECT) 


—-- He*(nl)+ H* (CHARGE TRANSFER) 


at He** projectile energies 1-20 KeV (c.m. energies 0.25-4 


Kev), and 


(b) H*+ He*(1s) —— > H*+ He*(nl) (DIRECT) 
—-»> H(nl) + He** (CHARGE TRANSFER) 

ftecom. energies 146-8. “Kev. sbasise sets with (op to ‘12 
molecular states (Table(VI-1)) have been used, and good 
convergence of results as a function of basis size has been 
found. These collision processes have earlier been studied 
extensively by Winter and Lane(71), Hatton, Lane and Winter 
(HLW)(74), Winter and Hatton (WH)(73), and Winter, Hatton 
and Lane (WHL)(75), using both the uncorrected PSS method 
and PSS method with Bates-McCarroll ETF corrections. 
Comparison with their results provides information about the 
influence of ETF descriptions on cross sections. Our results 
may also be compared with the available experimental 
measurements on these processes. 

For process (a), our total cross sections for charge 
transfer are in generally good agreement with those found by 


WH & HLW using Bates-McCarroll ETF's, except perhaps at the 
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Table VI.1 Molecular state basis sets 


HeH?* Molecular State Basis Sets 


for He-*-H(1s) and iene Ce) Collisions 


. of States Limiting 
2pm, 3do 


as above, PLUS 


2s0 


as above, PLUS 


1so 


as above, PLUS 
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highest energy studied. There is less quantitative but still 
reasonably good agreement with Winter and Hatton's cross 
sections for charge transfer to individual states of 
He’(2s,2p.,2p;,). Both our results and those of WH & HLW are 
in good agreement with experimental measurements of total 
cross section for charge transfer and the individual state 
cross section for He*(2s)(Gilbody's group) (see 
Figure(VI-18)). 

For process (b), on the other hand, some discrepancies 
between our values and those of Winter, Hatton and Lane(WHL) 
are found, using basis sets of comparable size, for both 
charge transfer and direct excitation cross sections; the 
charge transfer cross sections of WHL are 14 % to 32 % 
larger than ours, and their individual state direct 
excitation cross sections (He* 2Ee 20a 20 are 
systematically larger as well, in some cases by as much as a 
factor of two. These differences appear to be traceable to 
differences in ETF description used and the resulting 
coupling matrix elements; given the convergence behavior 
found in augmenting the basis from 5 to 10 states, it is our 
view that the "converged" values obtained by the two methods 
will in fact be significantly different. Comparison of the 
charge transfer cross sections with measurements by Peart, 
Grey and Dolder(85) suggests that our values are in slightly 
better accord with these than the values found by WHL, but a 
more definite conclusion on this question will probably 


require further experimental and theoretical investigations. 
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We will discuss the results for process (a) 
calculations and then those for process (b). However, we 
will give first a discussion of the basis sets used, and the 
coupling matrix elements, since these are common to. both 
calculations. 

Figure(VI-1) depicts electronic binding energies vs. 
internuclear distance R for 22 molecular states of HeH?*; 
this diagram is useful in understanding coupling and 


excitation paths in this system. 


A. BASIS SETS AND COUPLING MATRIX ELEMENTS 


1. Basis Set 

The basis sets used were selected as follows: 

(ay 2pq@- 2p 7 ;3d@ and 2sCVr states form “an essentyel 
first block of basis states for He?*-H(1s) collisions, 
and (although it plays a negligible role in process (a)) 
the 1s@ state must be added to these for a treatment of 
H*-He*(1s) collisions, where it is the initial state. 
(b) In their calculation, WH included as the next block 
all remaining states associated with principal quantum 
number n=3 (united atom), that is, the states 3po ,3p 
Tw ,3am ,3d§ , and 3s@ . The first three of these 


States couple Strongly with the Gd@- state. On the other 
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handsnt Winter) sandshHatitontfound thats the 3s@ and 3d$ 
States make a very small contribution to the total 
charge transfer cross section in process (a), and I also 
found after comparing coupling matrix elements computed 
here that coupling to these states from the first block 
of states should be quite small. Therefore I excluded 
these states from the 10-state basis, and added instead 
the states 4d@ ,4fN , which correspond asymptotically 
to He?*+H(n=2) levels. Inclusion of these states gives a 
rough indication of the role of H(n=2) production in 
processes (a) and (b). 

(c) Selected calculations of molecular state 
probabilities at various impact parameters and energies 
have been done using a 12-state basis consisting of the 
ten’ states! above; plus 4i@ytand 5q0rm. 1 4fi7~sisy strongly 
couplede-tonostatesolin the 10-state block, and its 
inclusion amplifies excitation to the 4f7] state (via 
rotational coupling) by factor of as much as 4 or 5. .5g 
O~- was eneluded because it plays an important role as an 
intermediate state: when 5g@® is in the basis, there is 
a significant increase in H(n=2) populations, even 
though the 5g@® population itself is negligible at the 
end of the collision. A similar effect or effects was 


noted by WH & HLW. 
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2. Dominant Couplings and Excitation Paths 

The two factors which determine effective coupling 
strength, i.e. coupling matrix elements and the associated 
energy gaps, are shown for the dominant couplings linking 2p 
G- (the initial channel®for process (a)), 2pMihand 3ace 
States* in ,Figure(VI—-2). As for H*-H{1s) colUisiogs, there 
is very strong rotational coupling between 2p@ and 2p7t 
States! in’ the uniteds, atom limit “due to their orbital 
degeneracy and this coupling is the primary one for process 
(avn A) Second’ ‘important,’ Link’ is «the padialy 2pG -3dc~ 
coupling which is strong around R=6-7 a.u. because there is 
a fairly small energy gap between these states at larger R, 
and this connection may be expected to be a dominant one at 
larger impact parameters. Also given in the Figure is the 3d 
G@ -2pm ‘angular ‘coupling matmix “element, showing the 
typical long-range Coriolis coupling associated with orbital 
degeneracy of the separated atom; the problem of asymptotic 
couplings is most conveniently handled separately (see 
Section(IV.D)). The molecular state 2s@ is not coupled 
strongly to the above three states at small R values, but te 
must be considered in the asymptotic coupling problem. 
Coupling to 1s@ is very weak due to the large energy gap to 
all other states, hence charge transfer to He*(1s) is almost 
negligible in process (a) and the total cross sections for 
process (b) charge transfer are also very small in 


comparison to those for process (a). 
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inmeprocess (b), the most important coupling out of is 
G- is the is@ -2p@, coupling, secondarily the Is@ -2p 7 
angular coupling, and (much less important) 1s@® -2s@ 
coupling (see Figures(III-10-15)). After these weak first 
steps, the strong 200-20 N |) 20G woe a 7 207-30 oO 
couplings play an important role in determining the further 
path of excitation, similar to their role for process (a), 
and most of the important transitions in process (b) involve 
two-step mechanisms in which 2p@ is the intermediate or 
"gateway" state. Figure(VI-3) depicts in a qualitative way 
the main coupling pathway for excitation in both processes 
(a) and (b); thickness of the arrows indicates qualitatively 
the coupling strength. 

For excitations beyond these five molecular states, the 
3dQ® state plays a most important role as a "gateway", since 
it is correlated with united atom manifold n=3 and is 
strongly linked to 3d 7) ,3pC@ 7, and opie (e@yese “States are 
of course also coupled strongly among themselves). As was 
found by Winter and Hatton, 3d6 is less important (it can 
only be reached through 3d7 or 3p7 ) and 3s® still less 
important. 

AS one “goes still) furthers up ankenergy,7the smaller 
energy gaps and greater number of states leads to a very 
rapid increase in the number of couplings which can play a 


(minor) part. 
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Figure (VI-3) Effective couplings and excitation paths for 
He?*-H(1s) collsions. ——» radial coupling; -<-—-=»> 
angular coupling; thichness of connecting arrow indicates 


qualitative importance. 
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3. Asymptotic Coupling 

Within an asymptotically degenerate atomic manifold---- 
such has she ae’ (n=2)..manhholdy,of, .25.00y ,34d 0a.uands 2p JT 
States---- a long range coupling problem arises which we 
found convenient to treat separately since these states are 
completely decoupled from other states for R215 a.u.. The 
two hybrid) (sp) @snstates -are.|linked «sto, .2pis by, Corioli 
couplings and to each other by long range radial coupling, 
SnegebotnecoupLingSsand Splittings decrease only as eRi7.. Hor 
ReiO med «ap aeDOth couplings and splittings can be accurately 
modelled by an analytical perturbation expression and, using 
a new progress variable, the numerical problem may be solved 
Vetyeset ficiently. ‘The,setfects uot this coupling are 
Suanuticant; for example, “at Ss Kev, the probabilities for 
He*(2s) excitation as computed by this procedure differ by 
about 25 % from those found by direct numerical integration 
of Eq.(II-20) even as far as R =25 a.u.. The details of the 


method are given in Sec.D. of Chap.IV.. 


4. ETF Effects on Couplings 

In this section I compare effects of different ETF 
descriptions on the resulting coupling matrix elements. for 
the HeH?* system (comparison with uncorrected PSS couplings 
is not so much of interest, since the work of WH & HLW shows 
that ETF corrections of some kind are necessary). We compare 
our coupling matrix elements (based on Switching functions 


"optimized" for each molecular State) with those of WH & HLW 


-veo"Bigiiaa® oinoj¢ oe pnt 
ne Beis “Bht. Wat a6 axotiti: i | 
at ide eeetusn aeldo3q gnliquas ‘gen’. ood - 
“14 #6fn9e eeed? a2kle eR cee baal an: aie 
s @0s8 191 sedeia vealzo ‘inci “betquened 

> yd Wes of eaeeply ae agteae atu 
.woliques feihesopae? eaet’ eg yeu ag doaw et ve 
4o8 8.'" eo vino seeensei eontayitge pene ) 
eistvens #d. nes fpmi tts.ige ihe apaitques: idee. 1% ie 
piiev .bas sOreanigad i Seder ee depiaytane ts wb 


3 i 


bavige ed yeu aefdos Land 1éitan Ssiz elidel F) 2xPOTR a 


Le 


: ; _ 
sia paifeues eid? fe rivdthe wnt ilar " 


e, a 


102 aardllidadgag ond vam a as \ 
4d oyettib Mwbhsgoiw’ wide a pesuanod. Ra wiaiasions | 

of seseasat teed tetuw saedib yd bauwd spas 038, 2a od 
7 YS athedet om ..toeS= © am el ae oeve (OST tame 
| ME Gat Ae coal aalecr be 

 egailgied a0 yarn 

WH srs79hbih i¢ evadedis otegees Tt nokzoes aide Ma 
rol adviemete xitgew Amiiques pristowed gai 9 anotsgty3 
sunt iguen 229 beiveedozn dete neni oomene? maseye nan ite 
svete Wi 2 WW te Saoe ody oonde ,deexagad- pcr 
ev addine > ow ,lviseeereh sie BAIR ence te meee v 
bed trond” ‘Cnidetivg no beendd ita 
as," 2. HW 90 sg0da dotw (sede Gelgowion- 


ms 


seat 


" 
ae 


196 


who used the simpler Bates-McCarroll (atomic state) ETF's, 
and with some matrix elements we have computed using 
Taulbjerg and Vaaben's(28) "universal" switching function 
(the same for all states). The couplings shown are for the 
dominant couplings in processes (a) and (b) (or for the 
primary excitation steps even though the coupling is_ weak, 
as in (b)). Even more pronounced effects of different ETF 
descriptions may be seen in couplings of secondary or more 
highly excited states, but these are not so significant for 
the collision processes considered. 

Figures(VI-4) & (VI-5) show coupling matrix elements 
Sinking 2pVwwith 2po7] wand) Saw gtheses are the ‘dominant, 
primary couplings in process (a) and are also important in 
process (b). 

For thesstrongu2p¢o- -Zp yi Peouplings (Figure(V1-4)) our 
matrix elements and those of WH & HLW, or Taulbjerg and 
Vaaben, are different mainly at distances between 2 and 10 
as. .wrhis coupling. is strong because of the united aivoem 
Orbital degeneracy of 2pVameand 2Zpjl) =; the ditferences 
between the coupling matrix elements beyond R=2 a.u. should 
therefore have little effect, except at higher energies, and 
this is the main reason why the total charge transfer cross 
section for process (a) is relatively insensitive to the ETF 
used, especially at lower energies. 

For the radial 2pa7 -3d@ coupling (Figure(VI-5)) our 
matrix elements are generally quite different from those of 
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Figures(VI-4) & (VI-5) Coupling matrix elements foOr=2pC> —Zp 
Tr angular coupling and 2po™ =300) radial coupling, 
respectively. Solid curves, matrix elements used in the 
present work, based on the different switching functions for 
each molucular state; dashed curves, matrix elements used by 


Winter et al, based on Bates-McCarroll ETE’ Ss. dot curves, 


used by Taulbjerg et al, based on the common switching 


function for each molecular state. 
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the important role at larger impact parameters and can be 
expected to have some effect on detailed excitation 
probabilities. Our couplings are generally smaller than 
those found using Bates-McCarroll ETF's, but the _ gross 
differences found at smaller R-values really have less 
importance than may appear because of the dominance of 2pa@ 
e2py\ coupling at small’ R. 

Figures(VI-6) & (VI-7) show the primary excitation 
Pouplings »tor sprocess.-(b), linking lsQ-with 2pc- 52pm; 
and 2s@ states; although these couplings are much_ weaker 
than those for process (a) (because of large energy gaps), 
they have a controlling influence onjall excitations in this 
collision process. 

The 1s0: °728sC@ jandi (2p Tl esrares Ware “all. -"He’"-like 
states; HO this reason, the couplings found using 
Bates-McCarroll ETF's are Hermitian (non-Hermitian couplings 
only occur for states with different ETF's), and also the 
couplings linking them are relatively insensitive to ETF 
description (i.e., our "optimized" ETF's are nearly the same 
as the Bates-McCarroll ones in the region of wavefunction 
overlaps). For the 1s@ -2pT couplings our matrix elements 
have maxima at a somewhat larger distance than that of WH & 
HLW but the shapes are generally similar. 

For the 1s@ -2p@ couplings the situation is quite 
different. 2pd> Vs'a “truly "molecular ™*state torr s*4 ta ru; 
with substantial two-centre character, and the 2p@ 


switching function used in our calculations reflects this 
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Figure(VI-6) Coupling matrix elements linking 1sO and 2po~ 


HeH?* states. Solid curves, matrix elements used in present 


work, based on switching functions; dashed curves,matrix 


elements of Winter et al, based on Bates-McCarroll ETF's. 
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Figure(VI-7) Other coupling matrix elements linking 1sO@ to 
excited states in HeH?*. Solid curves, present work;dashed 


curves, matrix elements based on Bates-McCarroll ETF's. 
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fact; the Bates-McCarroll ETF corresponds to setting fapq =+1 
everywhere. AS a result, our matrix element for Lise 
|RAD| 2pa> (2p~ 1s 0 ) is much smaller than that calculated 
by WH & HLW. Even more important however is the complete 
difference in shape and character of the element <2pa 
|RAD| 1sA> (1sA—2pA ) resulting from the use of Switching 
functions. AS can be seen from the parameters given in 
Table (Lir=2) 5 our VstO> = “switching function differs 
Significantly from the Bates-McCarroll value fjsa =—1 at 
small R values, even for this "He*-like" state, and this is 
especially true in regions of important overlap with the 2p 
(\ state. 

It appears that the difference between our couplings 
and those of WH & HLW for 1sH@e 2p@ (especially 1s@®-% 2p 
o~ ) is the main reason for differences in computed cross 
sections for process (b). This is at first Surprising since 
the sizes of the two matrix elements are Similar; however, 
we have found in the course of studies on couplings in 
H*-H(1s) collisions that shape has an effect which 1s as 
important, or more important, than integrated area, 
especially for non-resonant couplings: steeper slopes, or 
the occurrence of sign changes, have a tendency to increase 
transition amplitudes. (Such an effect can be understood 
from a Fourier analysis point of view: the more efficient 
couplings are those with components at frequencies closer to 


the frequency difference of the nonresonant coupled states). 
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B. CQ - PARTICLE - H(1s) COLLISIONS (PROCESS(a) ) 


1. Excitation Paths 

Typical collision histories are shown in Figures(VI-8) 
through (VI-11) for process (a): probabilities P(f:E) for 
molecular states (in a four-state basis calculation) are 
shown vs. time t, at impact parameters , fe 0.5 and 5.0, and 
0.3 and 3.0 a.u. and energies 8, 3 KeV, respectively. For 
all impact parameters the effect of 2ph -3de radial 
coupling is evident during the early part of them collision, 
but for small impact parameters this is later dominated by 
the strong 2p@ -2p/{ rotational coupling at small 
internuclear distances (Figures(VI-8) & (VI-10)). For large 
f , on the other hand, the effects’ of fthe 2p? —-30dl? 
coupling are dominant; for example, here even the 2p7\ 
excitation results in most part from the two-step process 2p 
Oo -3d® -2pm . (The curve crossing of 2pjt and 3d@ at 4.5 
a.u. may be important for the efficiency of this mechanism). 
The small 2s@ probability for both large and small impact 
parameters is attributable to the processes 2p7 -3d? -2s- 
Q and 2po -2p7i -2s@ , rather than direct 2pg -2sa” 


excitation. 


2. Probabilities vs. Impact Parameter 


a. Molecular states 
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Figure(VI-8) Collision "history" (molecular state 
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probabilities vs. time (a.u.)) for E=8 KeV at impact 


parameter 0.5 a4.u. 
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Figure (VvI-9) Collision "history" (molecular state 


probabilities vs. time (a.u.)) for E=-8 Kev at impact 


parameter 5.0 a.u.. 
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Figure(VI-10) The same data as Bigss(Vieo) &s(Vi—-9), but for 


E=3 KeV at impact parameter 0.3 a.u.. 
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Figure(VI-11) The same data as Figs.(VI-8) & (VI-9), but for 


E=3 KeV at impact parameter 3.0 a.u.. 
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Excitation probability vs. impact parameter 
at E=6 Kev is “shown “in ,Figure(Vl-12pm The large 
peaks at small impact parameter is due to 2p -2p 
MT angular coupling . The 3d@ _ probability shows 
Rosen-Zener-Demkov oscillatory peak structure at 
larger impact parameters resulting from coherence 
between the 2p@ and 3d@ states in the extended 
region over which they are coupled. There are 
secondary 2p7’ peaks in phase with these, which are 
due to the two-step process 2p@ -3d@ -2pt . The 
small 2s@ probability shows its two-step origin in 
the same way. There is no coherent interference 
between the amplitudes produced by 2p@ -2pn and 2p 


O- -3d@ couplings, as was also noted by WH & HLW. 


Total charde transicex probability 


The total transfer probabilities times impact 
parameter have been plotted versus f in 
Figures(VI-13) through (VI-15) for 1,8 and 20 Kev, 
respectively, and for both 4-state and 10-state 
basis sets. These figures may be directly compared 
with figures of the paper by Winter and Hatton(73). 
For 20 Kev (Figure(VI-15)) we have included Winter 
and Hatton's 10-state results for direct comparison; 
at lower energies, the agreement between their 
10-state results (not shown) and ours is even better 


than at 20 Kev. 
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Impact Parameter ( a9) 


Ayqeqosd 


Figure(VI-12) Molecular state excitation probabilities vs. 


impact parameter for He?*-H(1s) collisions at 8 KeV. 
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Figures(VI-13 through 15) Probability times impact parameter 
vs. impact parameter at 1, 8, and 20 Kev, respectively. 
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Figure VI.14 Probability xX impact parameter vs. impact 
parameter 
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Some general trends may be seen in the Figures: 

(i) fxP(P:E) for the 4-state and 10-state results agree 
closely at low energy and disagree at higher energy, 
showing that excitation to higher levels becomes more 
important at higher energies. However, even at 20 Kev 
the disagreement of 4-state and 10-state results is only 
12 % at the highest peak maximum. 

(ii) As noted in (i), above, the peak at small impact 
parameter (ficcae Saku, iis: “Drimarr ry: due YC angular 
coupling of 2p@ -2pim states. 

(iii) The larger impact parameter peaks (Piet ar Mare 
primarily due to radial Coupling Of -2p0 =" -3d7~ ‘states. 
and show typical Rosen-Zener-Demkov oscillations 
associated with coherence of the 2pm and 3d states 
which have an almost constant energy gap over = an 
extended region. 

(iv) These two different peaks caused by angular~ and 
radial couplings do not interfere coherently. The two 
major couplings occur in different regions and the 2p¢ 
-3dQ@ coupling is weak at the smaller R-values where 
rotational 2p@ -2p71 coupling plays a part. 

(v) The phase difference in the oscillatory peaks at 
smaller if between the 4-state and 10-state basis results 
is mainly the result of strong coupling from 3d@ to the 
upper levels 3p S2onee san © 43d0 esta "gateway" to 
these higher states, and as a result of this coupling 


the effective phase of the peak amplitude is shifted. As 
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waS noted some years ago by Knudson and Thorson(42), the 
influence of couplings gtot\lothers statestais toften 
detectable most easily in their effects on _ phases, 


rather than the magnitudes of excitation probabilities. 


3. Cross Sections 

In this section, we present our results for the total 
as well as the individual (n,1l) cross sections calculated by 
using 3,4,10,;12 states as in Table(VI-1). Emphasis is placed 
on the comparisons of the cross sections and the convergence 
of the cross sections in terms of the basis sets with other 
theoretical results particularly those of HLW & WH and 


Taulbjerg et al(72). 


(a)Total Cross Section 

Table(VI-2) shows our total charge transfer cross 
sections for 3-, 4-, and 10-state basis, together with 
results of previous calculations. There is generally good 
agreement at most energies between our results and those of 
HLW; like them we find that convergence of the results as 
function of basis size is much improved (relative to the PSS 
results of Winter and Lane), when ETF effects are included. 
Setting aside for the moment the results reported by Vaaben 
and Taulbjerg, it also appears that the convergence limits 
off ¢ athe ETF-corrected calculations are Significantly 


different from those of the uncorrected calculations, at 
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least at the higher energies. At 8 and 20 KeV the 
convergence of our total cross sections is somewhat faster 
than that of WH & HLW (changes from 4-state to 10-state 
basis cross sections are 64 % of those of WH & HLW in both 
cases); however, the comparison is not clear cut, since the 
10-state bases differ slightly. Additional information about 
convergence is provided by results of calculations using 12 
basis states at 8-10 values of impact parameter C120 a < Ua) 
Lor both energies. Table(VI-3) shows a comparison of 
molecular state excitation probabilities for fhe, tZ2-state 
and 10-state basis. It can be seen that in every case the 
total charge transfer probability has changed by less than 1 
%, relative to the 10-state result. 

At 20 KeV, it appears that the convergence imiutss 1.0L 
total charge transfer cross section resulting from the 
differnt ETF treatments by ourselves and by WH & HLW may be 
diverging slightly also. As shown in Figures(VI-13) through 
(vI-15), the main differences result from collisions at 
small impact parameter and are probably associated with 2p 
OQ. -2p7m coupling. However, calculations using still larger 
basis sets would be needed to settle the question clearly. 

The results reported by Vaaben and Taulbjerg(72), using 
a common switching function for all molecular states, are in 
marked disagreement with our values and those of Hatton,Lane 
and Winter, especially at the lowest energies, where the 
effects of ETF corrections appear to be least significant 


for this system and where all other calculations reported in 
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Table VI.2 Total charge transfer cross section 


Total Charge Transfer Cross Sections: 


He +H(le) Heth (Cunitesl0ce” em-} 


Basis fre Energy, keV 


Piacentini and 
Salin (1977) rare aamale we SareiErs 


0.994 
Winter and ; 1.03 
Lane (1978) 

22 


Calculation 


Vaaben and 10.36 


Taulbjerg (1979) 


- 5.19 
5.93 
Hatton, Winter 6.07 
and Lane (1980) 
6.30 
5.90 
Kimura and 6.01 
Thorson (1981) 
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Table VI.3 Convergence of excitation probabilities 


Convergence of Excitation Probabilities for the 


10- and 12-State Basis. 


of 
He? energy E=20 Kev 


Impact Parameter (a.u.) 


2.0 3.0 4.0 


0.765 0.801 0.799 
0.806 0.805 


] 
| 120 
= = 
i} 


0.573 0.768 


S50 


12-state 


\f 7 


reeat 


ia | 
tittle oy » Nt 
vos ; 
e.¢ 
—- -~ 
t j fy 
te 


idetdag occas aml 


Mee 


1a; euhsellenMeR nob nd hGH 


jak eit, goto Yeah 
| ek CRAG varies et 


a 


218 


Table(VI-2) agree most closely. These results were 
sufficiently disturbing that we performed a calculation at 1 
KeV uSing their switching function. For the 3-state basis at 
jitey we find the “value 0.247 X° 10°$°% yemz, ini general 
agreement with other values in Table(V1-2) gOur? result was 
obtained using straight-line trajectories, while Vaaben and 
Taulbjerg used the HeH?* 2p potential curve to determine 
trajectories.* However, the extensive studies by Winter and 
Lane on the effects of Coulombic vs. straight-line 
trajectories, show that these are relatively small, and we 
think it unlikely that the discrepancy 18 a trajectory 
effect. 

Figure(VI-16) compares our 10-state results and those 
of HLW with experimental measurements of the total charge 
exchange cross section. Except at the lowest energies (<4 
KeV), there is good agreement of the theoretical results 
with the values reported by Shah and Gilbody(83) and by 
Nutt, McCullough ,Shah and Gilbody(84). The measurements of 
Bayfield and Khayrallah(82) lie above the other experimental 
values as well as above our result and those of HLW & WH, 
for energies above 10 KeV. At lower energies, the older data 
of Fite Smith and Stebbings(80) lie considerably above the 
values given by Gilbody et al, who have attributed the 
discrepancy in part to _ the nonnegligible formation of 
molecular species, not taken into account by Fite et al. 


2Dr.J. Vaaben (private communication) has confirmed that a 
result close to the above value is also obtained by him 
using straight-line trajectories. 
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Figure(VI-16) Total charge transfer cross section for 
He?*-H(1s) collisions, vs. projectile energy E. Theoretical 
TeSU1tS > cums , DOSEN work; memeeeec= ,Winter et 
al(73,74); m———, Piacentini et al(70). Experimental 


data: 01,Fite et al(80); @,Bayfield et al(82); 0, Gilbody 


et al(83,84); M ,Olson et al(81). 
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The success of molecular state calculations in 
accounting for the marked decrease in charge transfer cross 
section seen experimentally below 10-20 KeV confirms’ their 
appropriateness for collisions at low energies. Atomic state 
basis calculationss by Rapp et al and by Gallaher et al (not 
depicted in Figure(VI-16)) show reasonably good agreement 
with experiment above 15 KeV, but below this energy give 
much larger cross sections (about one order of magnitude at 
1 KeV), agreeing only with the older data of Fite et al. Our 
results practically coincide with the measured values 
reported by Gilbody and coworkers in the region 5-20 Kev. 

Below ~4 KeV, the experimental values given by Gilbody 
étralefall off |istilis” imore, rapidly) than’ “the* theoretical 
values both of ourselves and of Hatton,Lane and Winter, 
while the two theoretical results are in excellent agreement 
at lower energies. Given the difficulties of experimental 
measurements at such low energies, we believe the 
theoretical results are more reliable in this region. 

As waS pointed out by Winter and Hatton, when ETF 
effects are taken into account the four states 2p® ; 2p7 
,3d@ ,2sG are almost completely sufficient to study charge 
exchange and to determine the total cross section 
quantitatively, since the levels He’*+H(1s) and He*(n=2)+H’ 
are near-resonant. Our calculations and those of HLW & WH 
for total cross sections establish that inclusion of ETF 
corrections has an important practical effect on 


convergence, and that uncorrected calculations probably do 
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not converge on the same result as corrected ones (except 
perhaps at the lowest energies considered, where only 2p% 

-2p N coupling at small R (insensitive to ETF corrections) 
has much effect). Moreover, there are small but not 
negligible differences in the apparent convergence limit 
values obtained using different ETF descriptions, at least 
at the higher energies, which suggests that in general a 
more elaborate description of ETF's for molecular states may 


be needed than is provided by Bates-McCarroll ETF's. 


(b)Detailed Cross Sections 

Cross sections for electron transfer to individual 
levels He*(n=1), He*(n=2), He*(n=3), and direct excitation 
to H(n=2) are shown for the 10-state basis in Table(vI-4) 
and Figure(VI-17). By far the largest contribution to charge 
transfer comes in the He*(n=2) mee eld due to the 
near-resonance with H(1s)--99.2 %, 95.6 % and 93.9 % at 1, 
8, and 20 KeV, respectively. A secondary contribution is 
made to. He (n=3))1 000 eo mance «buy. respectively. 
Capture into the He*(1s) ground state is unimportant at 
these energies. Direct excitation of H(n=2) amounts to about 
1.5 % of the charge transfer cross section at E=8 and 20 
KeV; however, these higher level cross sections may be 
expected to change significantly if a larger basis set is 
used. We find, for example, using the 12-state basis set at 
8 KeV, that the population of H(n=2) is increased by ~40 i%, 


mainly because the inclusion of 4£@% greatly increases the 
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Table V1.4 Atomic level excitation cross sections 


Atomic Level Excitation Cross Sections 
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Figure(VI-17) Detailed level cross sections vS. E.  __ 
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4£ 77 population via the rotational coupling between them; a 
secondary effect is introduced by inclusion of Sga -, (even 
though the final 5g@ population is itself negligible. 
(Similar effects are mentioned by WH & HLW). Expansion of 
the” basis “to include such ‘states as 4d? ,3d6 ,4dil , etc 
would be likely to shift upper level populations back (from 
H(n=2)) toward He*(n=3,4) levels again. However, these 
results do confirm that (when ETF effects are included) the 
4-state basis is sufficient for a qualitative account of 


He?*+H(1s) collisions at these energies. 


(c)He(n=2) Manifold Cross sections 

Table(VI-5) presents an analysis of individual atomic 
state cross sections within the He*(n=2) level. 

At low energies the principal contribution comes from 
2pyn excitation associated with the 2p” -2p7N rotational 
coupling, but as the energy increases this is overshadowed 
by the radial 2p® -3d@ coupling and , to a lesser extent, 
the 2p@ -2s@ coupling, which populate the 2s and 2p, 
atomic states. He’*(2s) and He*(2p,) cross sections are 
sensitive to the basis size and energy since as already 
discussed 3d@ state plays the role of a gateway to higher 
Rydberg states. Increase in the basis from 4- to 10-states 
affects the 2p, cross section much less than those for 2s 
and 2p, . This happens because a number of higher’ states 
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Table VI.5 He*(n=2) atomic state cross sections 


+ 
He (n=2) Atomic State Cross Sections 


He-*+H (1s) Collisions (units one cat) 
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States are relatively isolated from the higher levels 
(ote Figure (vi- 1) )-, 

Figure(VI-18) shows a comparison of our and other 
theoretical results for the He*(2s) cross section with 
experimental values obtained by Bayfield and Khayrallah, and 
Shah and Gilbody. Our results and those of Winter and Hatton 
are in good agreement with the experiments, especially those 
of Gilbody et al above 5 Kev, although our values are 
somewhat larger than the experimental value: as already 
noted by Winter and Hatton, however, the energy dependence 
at ‘lowrenergies isf fin better accord with the jtrend of 
Bayfield et al's measurements. The 3-molecular state results 
by Piacentini and Salin(70) are also shown. Also as noted by 
WH ,the better agreement of their values with the experiment 
results should be considered fortuitous, since their total 
Capture cross section does not agree with the experimental 
values and they did not test the convergence of the cross 
section as a function of basis set. 

Figure(VI-19) shows the individual 2s,2p.,2p, cross 
sections along with He*(2s) excitation cross sections 


measurements. 


(d)Effect of Higher Order Terms in Collision Velocity 
In the study of He?*-H(1s) collisions and H*-He*(1s) 
collisions, we used the coupled equations(Eq.(I1-20)), which 
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Figure(VI-18) He*(2s) capture cross section vs. E for 
He?*-H(1s) collisions. Theoretical resultS: —__ 
,»present WOTK; m_oemeneaenm=,Winter et al(73); —_—_——— 


,Piacentini et _al(70). Experimental values: @,Bayfield et 


al(82); O,Gilbody et al(83,84). 
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Figure(VI-19) Detailed He*(n=2) level capture cross sections 


vs. E. Experimental data are the same in Fig.(VI-18). 
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Section(II.C.1b), the effect of second-order and higher 
terms in Vis much less important for the HeH?* collision 
system than it is for the H,* system, primarily because of 
the shorter range of couplings between states correlating to 
the size of the higher-order terms. We did some calculations 
of molecular state probabilities vs. impact parameter at 20 
Kev (4 KeV in c.m. frame), with the 4-state basis, using 
coupled equations accurate to second-order in V. The results 
are tabulated in Table(VI-6) along with the 4-state basis 
calculated values which include only the first-order in V 
term only for eight different impact parameters. 

Tor can be Seen trom fable (Vico) that win all cases the 
individual state probabilities change by at most 1 &. 

We have not made tests for process (b) at the higher 
energy of 8 KeV (c.m. frame); however, since our studies in 
H,* suggested that second-order effects increase in roughly 
linear fashion with increasing energy (see Section(V.D.1)), 
we believe that errors of about 2 % (at most) in individual 
state probabilities in process (b) might result due to 


neglect of second-order terms at 8 KeV (c.m. frame). 
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Table V1.6 Effects of second-order terms in velocity 


Effects of Second-Order Terms in Velocity on 


x 
Total Electron Transfer Probability (4-state basis set) 


Here energy E=20 KeV 


O25 1.0 


Abas) 255 3.0 4.0 
0.261 0.556 0.625 e 
0.261 0.555 
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2 
0.281 0.571 0.928 0.742 
0.625 0.279 0.570 0.927 0.742 
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The four states 2po, 2pm, 3do, 2so0. 
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C. He*(1s)-PROTON COLLISIONS (PROCESS(b) ) 

Since He*(1s) formation plays a negligible part in 
He?*-H(1s) collisions, it is not surprising that the total 
charge transfer and direct excitation cross sections for 
H*-He*(1s) collisions are very much smaller than those for 
process (a) --- as has already been shown by Winter, Hatton 
and Lane(WHL)(75). The main reason for this is the highly 
nonresonant energy transfer required to go from is@ (the 
incident channel) to all upper states. Moreover, the 
couplings responsible for such excitation (primarily 1s@ 
-2p@ and to a secondary extent, 1sS@ -2pM™ ) are effective 
at relatively small internuclear separations (< 5 a.u.), and 
one might expect molecular state ETF'’S to give some 
different effects than a Bates-McCarroll description does. 
In fact, we find in this case that the treatment given to 
ETF corrections has a more significant effect on the values 
found for both charge transfer and direct excitation cross 
sections, than is the case for process (a). Our 5-Sstate 
basis and that of WHL are identical; as noted previously, 
the two 10-state basis sets differ slightly, but they do so 
only in respect to states which are relatively weakly and 
indirectly linked to the more important lower levels. In 
both cases we find that the cross sections computed by WHL 
(using Bates-McCarroll ETF's) are all systematically larger 
than those obtained here. 

We have already described the coupling matrix elements 


linking 1s@® to the manifold (2pA ,2pm , 2sf ,3ad?.) and 


(ialeononbe- . § ic 
1; dyey sidber tee ~s agente alma junk an: ink i : 
eo ENEC On {ark 

"73 enrertese et0e> nolsesiona Jeger th Seem iabatiess | 
seodt? madt 1effeme sou YIsv Sas enobethtes. ” ‘(sy hy 
nfiisH testo vc avedea nade giaesie anal ae nee Ke), a 
vidgin eto ef eid vel sosaee whem eet Sedan — 
zi Set mort dp d+ oettivge va Fanagdi Nog ti og 
lavesTomM veecede aeday 22k oF ied S — 


Gah ylisient 4) detguptoxe gina Lamas sceneiagoan 
Vijnertte asus | Pate Aal .. snes yisbnogen s of 


feces si3 feds puigetieie. Joa at 34 


ia «ted 2H) (Je rages waa Lounramta bhai: eave: 
moe evip of --' TRS onsite sehusetom, sree no on 
aooh Motigiuoash Licvadie ede 2 neds eioetie 4 oe 22 Eb 
wip Joentes?? aff sada sane etda nt boit ow To) 
eliaed ont ra ioette igead E! 15 e108 a eed anotsrex0a % 
anois nolisdinza- doaceh-tne welanend waxy dicg 10d 
ssada-G ‘wo he) aesuo7g wi eteb ane a nmet3 venels: 
‘Ciawoiverg beten es glesficeth) eve’ Si Ba sey Bos ela 
Oe Ob watt duc ere ae eth Stem whecd ogege~or ows wt 
Ges ylteew itevenepss ove Cotie eGinve of szseames at waa 
1 6 .Sievei «6©sewol jnastoomt “sian sd4 04 betore L ulsoestiad 
HA vd Betuqwios anoidose teats edo tars bit sv ReseD. djod 
ere: yilsobtenpedave {fe-s18 (2 ota i Lov ahegeces puieu), | 
| ‘o1ad Benietde gaody ~ 7 
Riesteis ‘itdam poilquoes 643 tediysest ybenstg-onalt ow : 


Baa { S6E, BES, WES. MAE) Bloltinee otgot Bae | ‘pnzande ; 


eo 


shown that the different ETF description resulting from our 
molecular state switching functions leads to markedly 
different coupling matrix elements from those of WHL (See 
Section VI.A.4). Furthermore, the most important differences 
appear form-the ='l'sf “-2p0 poupling.) which is the primary 


mechanism for excitation in these collisions. 


1. Excitation Paths 

Figures(VI-20) & )(VIF2i)f showa.coliision,) histories 
(rolecular State probability pvs. time (a.u.)0, using the 
5-state basis, at E=4 Kev (c.m. system) and impact 
parameters QD and 2.0 a.u., respectively. Dominant 
interactions, occur for |tisi2@ga.u. (or R <4°7 asu.), since 
the couplings (Figures(VI-6)) have short ranges. While some 
very small 3dO and 2s@ amplitudes can be seen early in the 
Collision aue . to’ direae Beouplamdgeetmon. is? {, both are 
eclipsed by the primary excitation process 1s@ -2p@ . At 
higher energy (8 Kev), direct 1s@ -2pyH excitation plays 
Some part but ( for example) the peak in 2pm probability in 
Figure(VI-20) (0.5 a.u.) results mainly from the two-step 
process 1SQ@ -2p@A -2pmM rather than a direct: “coupling. 
Similarly, population of 3d@ , and the (small) 2s@ 
population, come from two-step processes KGa 0;) SiGe = 2p Uae, 
followed by 2p "sac ‘radial’ coupling laters’ an stehe 
Collision). Tnus’ “a ‘Yarge’~ portivon eof the flux which 
eventually appears in He*(n=2) levels passes through the 
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10°? 


H* + He* (1s) 
p =0.5 au. 


10°? 


1G 


Probability ——~ 


Figures(VI-20) Collision "history" (probability vs. time) 
for 5-state basis, H*-He’*(1s) collision at c.m. energy 4 Kev 


and impact parameter 0.5 a.u. . 
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10°3 


H* + He* (1s) 
Ecu. = 4 keV 
p =2.0 au. 
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Figures(VI-21) Collision "history" (probability vs. time) 


for 5-state basis, H*-He*(1s) collision at c.m. energy 4 Kev 


and impact parameter 2.0 a.u. 
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intermediate; in this sense, the term "direct excitation", 
though it is the correct technical name for the excitation 
processes leading to He*(n=2), is perhaps a little 
misleading. This picture is not changed by increasing the 
basis to ten states, and it explains why both the charge 
transfer and "direct excitation" cross sections in process 
(db) are sensitive.to the 1s@ -2p(Q coupling. 

Figure(VI-22) shows impact parameter times probability 
Of charge ‘transfer tofPH(1s) (2p state) vs: impact 
Parameteru, at) Ec m.=47 Kev. for both 5-state and 10-state 
basis. All sthertransirions occur for le Seeorae 4g Whiech 
demonstrates the short range of the primary coupling, and no 
Significant range golf the primary coupling, and no 
significant change occurs when the basis is increased. Thus 
the five states’ (1sA  ;2p@a0 ,2p"e,cde ,2s7 ) are 
sufficient both qualitatively and quantitatively to evaluate 
the cross section for charge transfer ( to H(1s)). 

Probabilities for excitation to He*(n=3) and charge 
transfer to H(n=2) (found using our 10-state basis) are 
relatively small (less than 10 % of the total excitation 
probability in all cases). These more highly excited levels 
are accessed mainly via the 3d@ state as "gateway", just as 
in process (a), while the 2s7 ,2pf\ levels are relatively 


"Isolated" from these upper states. 
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Figure(VI-22) Probability times impact parameter vs. impact 


parameter at c.m. energy 4 KeV. 
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2. Cross Sections 
Table(VI-7) presents our calculated charge transfer and 
"direct excitation" cross sections for the 5-state (KT-5) 
and 10-state (KT-10) basis. Also shown for comparison are 
the values reported for their 5-state and 10-state basis by 
WHL (WHL-5,WHL-10). Note that for WHL-10 no basis state 
correlating to H(m=2) occurs, and also that cross sections 
for He*(n=3) were not reported in their paper. 
atm Convergence 
For the H(1s) charge transfer cross _ section, 
the largest change we find on augmenting the basis 
from five to ten states is -4.3.% (8 KeV), compared 
top £27" 99 “% 2.04" MeV) Pfor® -WHLG! “For | excitation of 
He*(n=2) levels, the stability of the individual 
26°2057 wand, 2b uno Osc sections wars «much poorer, 
particularly at 8 KeV (changes of +14 %, -26 %, ~-12 
vi) Tespectively ,coMlpared tO t13y%, -55 %, ~19- 1% for 
WHL).. The reason fornithismiostability Jot. cour Seiiu1s 
the substantial redistribution of flux into He’ (n=3) 
and H(n=2) levels due to coupling from 3d@ to 3d7 
,3p0, 3p7 (and 4d@ in our case), when these states 
are included; that it is 3d#& which serves as_ the 
connection is reflected in the somewhat better 
Stability of 2pm cross sections. A more significant 
measure of "overall convergence" is given by the 
changes: of valay tee; cand) 410% (at. obA6iy 4 and a6 


KeV, respectively) in the total cross section for 
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Table V1.7 Atomic level excitation cross sections 


+ 
Atomic Level Excitation Cross Sections for H He (is) Collisions 


6 maa) [Present Results: KT5, KT10; Results of Ref. 75: WHL5, WHL10) 


(units 1072 


EoM. Basis & 
(keV) Calc. 
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both charge transfer plus excitations of all kinds; 
these values make it clear that augmentation of the 
basis from 5 to 10-states has not caused massive 
changes in the fx drained from ist 3 
Unfortunately the corresponding data is not 
available from WHL. 

While it appears that the convergence of our 
resultssisVslightlyo fastereethan *thaet Sof**wWHL; = ea 
comparison is weakened by the differences in the two 
ten-state bases. In both cases, further 
augmentations of the basis may be expected to cause 
significant changes in the He*(n=2) individual state 
cross Séctidns® wOnestheseother.thand* it. omay~ be 
reasonable to expect that the charge transfer to 
H(1s), and especially the total flux loss from is@, 


will be much more stable (< 10 % change at 8 KeV?). 


Cross Section Values 


In our opinion, a more significant difference 
is! *found”™ in “comparisons =‘of” the’ ‘cross section 
magnitudes found here with those reported in WHL 
(see Table(VI-7)), and comparisons of both 5-state 
and 10-state results are relevant for this purpose. 
In every case, the cross section values reported by 
WHL are appreciably larger than those we obtain for 
the comparable basis. For the more stable H(1s) 
charge transfer cross section, the cross sections of 


WHL are 18 ‘eto we2,-% larger than / ours» for’ the 
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five-state basis, and 14 % to 26 % larger than ours 
for the 10-state basis. These differences are at 
least twice as large as changes occurring due to 
basis augmentation in WHL, and at least four times 
as large as those due to augmentation in our 
results. While it is conceivable that further, large 
augmentation of the basis (discrete states only) 
mi ghtuwieadsntowreconciliation ,ofe these differences, 
iia 1SPOULLODI DAOnRa tNatwethicmuwill (not, occur. In 
principle, there is no reason to expect that such 
convergence on a common value must occur, unless 
completeness of the two basis sets is explicitl 

ensured by inclusion of both continuum and discrete 
states in each. 

An alternative explanation (ope these 
differences, i.e., that they result from effects of 
higher order terms in the velocity, (we retained 
only first-order terms in our calculation), seems 
improbable to us because (a) the differences are 18 
%¥-14 % even at 1.6 KeV, (b) as stated previously, we 
Studied effects of the second-order terms MONS 
process (a) at EB. =4 KeV (20 Kev projectile energy) 
and found them less than 1% for individual _ state 
probabilities, and (c) the sign of these 
second-order corrections is always negative-- which 


would increase the observed discrepancies. 


Broo. ood taptal. £45 af or oi itu | 
+e ere senndiottPh oeed? ane | 
ot wh .onfiesse cepnats 
mid 205% 26a! 28 Bs zea “i note 
wo in] Ol sonraemmcs - ats yaad 
spisi ,rehtue? cade eidevteoton as ae ares, webu 
ig @3eeate sisvoeih) fea ait bias ed eaainRINpET ES ae i 
S4aNetSIIio wba} .%¢ Ros tat lPeteass me bees, apes 
tog. gon [Siw ids Semy me wo ak di) 


ivie FeaqJ 19G"e Qt oOfeeet 6@,,828 waeds aaah) 


iolinnaigite svi sonaeaie vs 
veer yey sef eothebs epagnayst ib 
Sewieds1 of) ,yaisctew ate af nated ‘sapn0 aedpid 4 
amyee. | no} Haleciee, aaard aut rebae aa ine | 
Si $26 Roane telIs6 age (al aduseet Bb. of aidadowut 


o4 .¥lavolvang Fogete ee Me) Aa ey) ta stove : Br-x 


= 
— 


te *«ioasiis mos? 3 


' 


tO; 0 ais iahiosbaadse hee * werselta Seibuse 
i(Npcans siisvepend v5) OSY Volt’ dim: ssa 26. (6) steceag | 
avace fevbivibat. sc? % 1 me Coat sey haved. ins 
feet; | 1D) pte. 6atle to) os eedabiddeder 


Gsidv ~~avizepan cyawds ci Bnet ioens0> wehtorhnones 
O29 Fonage toe te bayesian s2 sanmsant ati 


wee 2 Th 


241 


Similarly large or even larger discrepancies 
(Hsiimiich ‘asita- sfactor tof 2iifori2smand2payiiand any &% 
tick. 2288 'Y, Maton Gi2p A)woccin kt om the Hindi viduadedaiy ect 
excitation | cross sections. Given the poorer 
convergence in both sets of calculations for these 
detailed cross sections, it is not so clear that the 
values would converge on different limits; however, 
the consistency of the discrepancies found suggests 
this is a reasonable conclusion. 

Cumulative level cross sections for He*(n=3) 
and H(n=2) are also computed in Table(VI-7) for the 
10-state calculations. We do not attribute firm 
significance to these numbers, as they may be 
expected to change appreciably when additional 
States belonging to the same manifold are included 
in the basis. Our estimates of the contribution of 
H(n=2) are probably a little high; we have included 
H(n=2) states, but not the comparable He*(n=3,4) 
States which are their competitors for coupling to 
3a@- as excitation source, and inclusion of these 
States might be expected to shift flux back toward 
He* formation. We do agree with the conclusion of 
WHL (based on PSS_ studies) that the H(n=2) cross 
section is "small" (4 %-10 % of H(1s): at, 8: Kev). 

In summary, the comparisons made suggest that 
the different treatments of ETF corrections made by 


ourselves and by WHL may lead, for process (b), to 
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Srontiricantiy (different. >cross sections “both for 
Gharge transter and for direct excitations, “and “it 
iSeecur )Opiniony yehat thisiis imstact’ che case... (It 
should be possible to test this by calculations 
using a basis perhaps twice as large as those 
employed here.) 
c. Comparison with Experiment 

Experimental values for the total charge 
transfer cross section, aS reported by Peart,Grey 
and Dolder(85) are shown in Figure(VI-23) together 
with our 10-state results (including H(n=2)), and 
those of WHL. The experimental values have been 
determined from the measured He?* formation; hence 
they include both charge transfer and ionization 
events, though it may be reasonable to assume that 
the contribution of terization is less than 1%. In 
spite of differences between them, both our. results 
and those of WHL appear to be in good agreement with 
the experiments. Our lower values at 4 and 8 KeV_ do 
lie slightly closer to the data values given by 
Peart et al than do those of WHL. The same data for 
total charge transfer cross section for. both 
theoretical and experimental values are plotted on a 
linear graph in Figure(VI-24). It clearly shows that 
discrepancies between ours and WHL are’ significant, 
especially at 8 KeV. A conservative opinion on the 


error limits for the experimental data does not 
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Figure(VI-23) Total charge transfer cross section for 
H*-He* (1s) collisions, vS. c.m. energy E. Theoretical 
results: @,present work; O,Winter et al(75); MA Rapp et 


al(76): Experimental result: & Peart et al(85). 
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(10 cri) 


CHARGE TRANSFER CROSS SECTION 


Ecom, (KeV) 


Figure(VI-24) The same data as Fig.(VI-23), but in linear 


scale. 
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permit a strong conclusion as to which set of 
results iS more accurate. 

More accurate measurements of this cross 
section, or accurate measurements of the direct 
excitation cross sections, would be of porate ea nie 


interest for the questions discussed here. 
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